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Part I. InrRopucTION 


DISLOCATIONS were first introduced into physical theory as part of the 
study of a mechanical aether. Burton (1892) postulated that ‘a given 
portion of matter consists, not of any individual portion of aethereal 
or other substance, but of modifications in the structure... . of the 
aether, and when matter moves it is merely these modifications of 
structure . ... which are transferred from one portion of the aether to 
another’. As a typical modification he considered a ‘ strain-figure ’. 
He realized that these figures could only persist in a medium with 
non-linear elastic properties: ‘Let the medium ... . be held in its 
strained condition by some compelling agency: .... provided the 
strain has at no point too great a value, the original condition will be 
completely regained after the compelling agency has been removed. 
But suppose ... . that the restoring stresses do not always increase 
with the strain, but that- beyond a certain point in the process they 
begin to fall off in value, until at last a point is reached at which the 
general tendency of the stresses is to further increase the strain. If the 
compelling agency is now withdrawn, the medium will subside into a 
new condition of stable equilibrium, involving stress and strain at every ° 
point. The state of things ... . will be referred to as a strain-figure, 
and we may now proceed to examine its dynamical properties’. These 
include the ability to travel freely with any velocity small compared 
with the speed at which disturbances of the same type as the strains 
in the strain-figure are propagated through the medium, the possession 
of an ‘ effective mass ’, and the existence of forces between strain-figures. 
‘The coarse-grainedness of the medium has an influence in determining 
the size of possible strain-figures.’ Larmor (1897) used a similar model, 
considering as examples a wire ‘ welded into a ring after twist has been 
put into it’, and a ‘breach of continuity ’ produced ‘ by the creation 
of a lens-shaped cavity, and the material on one side of the breach is 
twisted round in its plane, and continuity is then restored by cementing 
the two sides together’. He pointed out later (1900) that ‘It is the 
intrinsic strain-form alone that constitutes the electron; and it is a 
fundamental postulate that the form can move from one portion to 
another of the stagnant aether somewhat after the manner that a knot 
can slip along a cord. If this form is taken as an entity, so that its 
position is part of the specification of the system, then the dynamical 
analysis introduces forces acting on it.... When an intrinsic strain-form 
is implanted in an elastic solid, . . . . it cannot in general be removed 
to another position even by imposing additional strains in the medium, 
because the discontinuities involved in its creation cannot be entirely 
annulled by imposing any continuous strain: thus there cannot be 
phenomena of freely mobile strain-forms connected with a solid medium = 

Larmor’s theory of strain-figures in a medium showing strictly linear 
elasticity, except at some time in the past when violent cutting and 
welding operations were permitted, has developed into, the theory of 
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elastic dislocations described in Part II of this article; Burton’s theory 
of strain-figures in a granular medium which does not obey Hooke’s 
law for large strains corresponds to the theory of dislocations in a 
erystal lattice described in Part III. In this article we deal only with 
stationary dislocations, and we do not discuss the influence of dislocations 
on the scattering of x-rays (Wilson 1949, Frank 1949 a) or on the 
electrical (Koehler 1949, Mackenzie and Sondheimer 1950, Dexter 1952) 
and magnetic (Brown 1941) properties of the material. Nor do we 
describe the important role of dislocations in crystal growth, of which 
accounts have recently been given by Burton, Cabrera and Frank (1951), 
Burton (1951) and Frank (1952 a). 

This article was originally written during my tenure of a Royal 
Society Warren Research Fellowship at the University of Bristol. It 
has been considerably improved by the criticisms of my colleagues in 
Bristol and afterwards in Birmingham: in particular, I should like to 
mention the continual help of Dr. J. D. Eshelby. 


Part II. DistocaTIONS IN AN ExLastic CONTINUUM 
§1. Tom DrvELOPMENT OF DisLocATION THEORY 


The ordinary theory of elasticity requires that the displacement at 
any point of a body should be a single valued function of the coordinates 
of that point. This restriction disappears if the body is multiply 


Fig. 1 


(a) (b) 

(a) Dislocation of a hollow ring. The ring is cut along the radial lines AA’ 
and BB’, the wedge AA’B’B is removed, and AA’ is joined to BB’. The 
same dislocation could be made by cutting along CC’ and DD’. (6) The 
ramen ring. The neighbouring points Q, and Q, lie on opposite sides 
of the cut. ; 


connected. For example, the ring of fig. 1 (a) may be cut along the 
radial lines OAA’, OBB’, and the wedge AA’B’B of angle ¢ removed. 


en 
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If the cut edges are now joined, the ring is in a state of self-strain. The 
coordinates of a point P in the strained ring may be expressed in polar 
coordinates r=OP, 6=angle AOP (fig. 1(b)). The displacement wu, of 
the point P is given by up=r6¢/27. For the point Q, the displacement is 
zero; for the neighbouring point Q, the displacement is r¢. We may 
either put a mathematical barrier along AA’, and say that the displace- 
ment is discontinuous across the barrier, or else say that the displacement 
is continuous, but not single valued. The possibility of doing ‘this 
rests on the fact that the physical state of any part of the body 
does not depend directly on the components of displacement, but only 
on their derivatives. The addition of a displacement corresponding to 
that of a rigid body does not alter the strains, and the difference of any 
two of the multiple values of the strain is a rigid-body displacement. 
If we use continuous multiple valued displacements, the component w, 
at Q, has the values 0, rd, 2rd, ...., —r¢d, —2rd, ..., and the second 
of these is continuous with the displacement of Q,. Hither of these 
methods of approach is mathematically satisfactory, but the continuous 
displacement gives a better physical picture, for there is nothing in the 
ring of fig. 1 (b) to distinguish the line AA’ from any other radial line. 
Instead of cutting out the wedge AA’B’B, we might have cut out the 
equal wedge CC’D’D, and the resulting self-strained ring would have 
been identical with that of fig. 1 (6). 

Displacements of this kind in isotropic elastic bodies were studied by 
the Italian school of mathematicians under the name of ‘ distorsioni ’. 
The principal results were summarized by Volterra (1907), and more 
recent accounts have been given by Love (1927), Neményi (1931) and 
Goodier (1939). Love ‘ ventured to call them “ dislocations ” ’, and the 
word has been generally adopted. 

The first application of dislocations to the theory of plastic deformation 
was made by Prandtl in a lecture course in 1921-2, in which he developed 
a model to explain the plastic properties of crystals by assuming that 
in regions of imperfection there were molecules which could change their 
allegiance from one relatively perfect lattice to another. Prandtl (1928) 
did not indicate the connection of these imperfections with the dislocations 
of the theory of elasticity, although the mechanical model he devised to 
represent his ideas contains one idea that was essential to the further 
development of a theory of crystal slip. The solid is imagined divided 
into upper and lower parts. The periodic field of the lower part is 
represented by a corrugated surface, while a particular atom of the upper 
part is represented by a heavy roller which can move across these 
corrugations. The roller is connected by a system of springs and weights 
to a rigid framework representing the rest of the upper part. A shear 
of the solid is represented by moving the rigid frame in a direction 
perpendicular to the corrugations. For small strains the roller moves 
reversibly in one trough; for larger strains the roller jumps irreversibly 
into another trough. If the shift of the rigid framework is equal to the 
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wavelength of the corrugations, and the roller jumps to the next trough, 
the forces acting on the roller return to their original values. The 
process may be repeated from trough to trough. 

In considering the plastic deformation of single crystals, Masing and 
Polanyi (1923) were led to consider the imperfections which must os 
in a bent crystal. They suggested a type of deformation (‘ Biegegleitung °) 
in which sheets of crystal bent elastically were joined through surfaces 
of misfit as in fig. 2. The lattices on either side of each surface of misfit 
are uniformly expanded or compressed, and atoms above and below the 


‘ Biegegleitung ’ as envisaged by Masing and Polanyi. The crystal breaks up 
into a pile of elastically bent lamellae. 


Fig. 3 


Dehlinger’s picture of a ‘ Verhakung’. Six rows of atoms in a simple square 
lattice. The atoms in the three upper rows are displaced : the sinusoidal 
curve shows the potential in which the atoms in the third row move on 
account of their interaction with the atoms in the fourth row. 


surface are out of register over almost the whole surface. It follows 
from the work of van der Merwe (§ 2.31 of Part III) that a structure of 
this kind would not be stable. A closer approach to the modern idea of 
a dislocation was provided by the imperfection which Dehlinger (1929) 
introduced under the name of ‘ Verhakung’ (fig. 3). In this, the 
disturbance of the lattice is confined to a small region of the surface of 
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misfit, and outside this region the atoms above and below the surface 
are in almost perfect register. On Dehlinger’s rough mathematical 
approximations the localization of the disturbance is even more marked 
than in a more exact treatment. His ‘ Verhakung’ in fact represents 
a pair of dislocations of opposite sign, separated by only one lattice 
spacing. Dehlinger’s idea was important in directing attention to 
localized imperfections, and in indicating that the lattice planes 
perpendicular to the glide plane must become curved in the neighbourhood 
of the dislocation. 
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Orowan’s model of a local gliding. (a) Section in the glide direction perpendi- 
cular to the glide plane AB. The lattice planes were originally straight 
‘and orthogonal. The rings indicate the positions of the dislocations. 
(b) Schematic view in the glide plane. The circles were originally 
concentric ; the dislocation lies between the circles A and B. 


Fig. 5 


Polanyi’s picture of a dislocation. The glide plane is AB. 


The foundations of the modern dislocation theory of slip were laid by 
Orowan, Polanyi and Taylor in 1934. Orowan (1934 b) showed (fig. 4) 
that if a limited area of one plane slips by one atomic distance over the 
neighbouring plane, the boundary of this area is a closed loop of 
dislocation. Once this loop is formed, glide can propagate across the 
plane by the spreading of the area, which is a motion of the dislocation 
line across its glide plane. Polanyi (1934) sketched the form (fig. 5) of 
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a single dislocation at rest in an unstrained lattice, and showed that it 
could move across the glide plane under stresses much less than that 
required to slide one lattice plane rigidly across the other. Polanyi’s 
dislocation is confined to a small region of the glide plane, but the strains 
within the region of the dislocation are still uniform, as they are in fig. 2. 
The planes perpendicular to the glide plane are not bent, and the 
separation of pairs of planes passing near the dislocation does not tend 
to the unstrained value at large distances. Polanyi estimates that if the 
dislocation is n atoms wide, the shear stress required to move it is 1/n of 
the shear strength of a perfect lattice. This would require ” to be at 
least 1000 in a soft single crystal, which seems improbably high. The 
estimate of §4.1 of Part III leads to a much smaller stress. If the 
dislocation moves across the crystal from one side to the other, the 
material above the glide plane slips one atomic distance over that below. 


Fig. 6 


Negative Dislocation. 


: : ne 
Taylor’s model of the atomic positions before, during and after the passage of 
a dislocation across the crystal. 


Taylor (1934), like Polanyi, considered dislocations lying along straight 
lines in the glide plane perpendicular to the glide direction. He pointed 
out the connection of these dislocations with the elastic dislocations of 
Volterra, calculated (though with an error) the stress field surrounding 
such a dislocation, and studied the interaction between dislocations. 
His sketch (fig. 6) of the positions of atoms surrounding a dislocation 
showed for the first time that nearly all the lattice planes perpendicular 
to the glide plane remain continuous, though they are curved. 

Burgers (1939), considering the connection between Taylor’s two 
dimensional lattice dislocation and the elastic dislocations of Volterra, 
was led to consider ‘ screw’ dislocations, which lie along lines in the 
glide plane parallel to the direction of slip, and Frank (1948, 1949 b) has 
emphasized the importance of this type of dislocation. Since then the 
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theory has been extended by a number of workers, but it is only recently 
that attempts have been made to relate dislocation theory to the actual 
crystal structures of individual metals (see Mott 1951, 1952a). Surveys 
of the dislocation theory of slip have been given by Cottrell (1949) and 
Seitz (1950 a). 


§2. Some ELEMENTARY PRoPERTIES OF DisLocATIONS 
2.1, 7HE STRESS FIELD OF AN EDGE DISLOCATION 


Consider an isotropic medium having rigidity » and Poisson’s ratio o. 
In this medium lies a dislocation of the kind sketched in figs. 5 and 6, 
which we shall call an ‘edge’ dislocation. Conflicting formulae have 
been given for the stress distribution around such a dislocation: the 
present analysis is based on that of Koehler (1941). We take the Airy 
stress function 


X= —$ Dy In.b (2? +-7")+ Ba tan-| (y/~), . . . . (1) 
in terms of which the components of stress are given by 
07x 
07x 
Pry — Oxoy 2 Se Tas ome ies (2) 
0"y 
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To obtain Taylor’s expressions for the stresses we take 


DE alte A ee ee eer (9) 
while for J. M. Burgers’s expressions we take 
D=bp/2r(1—oc), B=0. eee Merc ety sat) 


Tn these formulae 6 is the strength of the dislocation, equal to the relative 
displacement of the two parts of a crystal when the dislocation crosses 
the glide plane. If the stresses are calculated from (1), (2) and (3), it is 
easy to see that the y component of the resultant force transmitted across 
a circular boundary centred at the origin does not vanish, but has the 
value 27B=—2by for each unit length of the dislocation. This stress 
system can only be maintained by a force of — 2b on the boundary at 
infinity, and an equal and opposite force applied at the origin. Instead 
of calculating the stresses and integrating, we may use the theorem of 
Michell (1899) that the resultant forces acting on the boundary are 
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where the integrals are taken round any closed circuit containing the 
dislocation. Single valued terms in dx/dy and 0x/dx do not contribute 
to these integrals: only the term Btan~! (y/x) in 0x/dx contributes a 
residue of 27B. 

Both these forces at the origin and the corresponding couple, which 
depends (Filon 1921) on an integral of the type of (5) involving x, 
vanish for Burgers’s stress function, which thus represents an 
unconstrained dislocation in an isotropic medium. 

The displacements are given (Love, p. 205) by 


0 0 
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and the function 4=— D[x In b-?(2?+-y?)—2a-+ 2y tan-! (x/y)] satisfies (7). 
The corresponding displacements are 


ee b xy b a 
«~ 4r(1—o) zi Latta, coe y (8) 
(1—20)b b mj 
ma aided ay eg MY Wht ioe Re 
IT Be (hae) es eI delet ee 


The stresses are given by 


Pog = — Dy (3x? +y") (a? +y?)? 
Pry= Dx? —y’)|(eP+y?)? cath ea eet 
Pyy=Dy(P—/P)|(e+y*). 


The equivalent expressions in polar coordinates are given in § 4.41. 
These formulae, which are given by Burgers, are equivalent to those 
of Koehler, since the origin of uw, and u, is arbitrary. It appears that wu 
is logarithmically divergent at large distances from the dislocation. 
The strain and stress correspondingly tend to zero only as fast as 7} 
at large distances 7, and the energy of the strained medium is not localized 
near the origin, but diverges logarithmically at infinity. Because of 


this divergence, the mathematical theory of dislocations presents 
unexpected difficulties. 
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The displacement w, is single valued, but w, has a residue of b for every 
circuit enclosing the origin, corresponding to the fact that every circuit 
of this kind cuts once the surface across which the lattice has slipped by 
a distance b in the x direction. 

There are no other two dimensional displacements which have the 
essential properties of a dislocation. These properties are (a) that one 
component of displacement (which we take to be w,) should have the 
same residue for any large circuit enclosing the dislocation, because such 
a circuit cuts the line over which slip has occurred once, (b) that the 
other component u,, should be single valued, (c) that no resultant force 
or moment acts across a circuit surrounding the dislocation, and (d) that 
the energy should be as small as possible. All of these conditions, except 
possibly (d), are satisfied by the Burgers dislocation. If we work in polar 
coordinates, the stress function is of the form y=—DrInrsiné@, and 
the stresses are everywhere of order 1/7. Suppose another dislocation 
with the same strength b, but less energy, can be made by adding to y 
another stress function y1. Then stresses corresponding to y! must also 
be of order 1/r at large distances. For if the stresses tended to zero 
more slowly, the self energy of y! would diverge at infinity more rapidly 
than the interaction energy of y and y1, and the addition of y! would 
increase the total energy. If the stresses tended to zero more rapidly, 
the interaction energy would be negligible in comparison with the 
self energy of y. We therefore require that each of the expressions 


Beil dx 1 @y1 
Pre 7 Or 1 OE 
a /1 ax! 
Vi = ay (; x) . . . . . (10) 
C2 1 : 
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should be of the form r—'a,,(0) when r is large. This requires that x! 
should be of the form 
yi=Fr In r cos 0+Gr In r sin 0-+-rH (6), 
with stresses 
Dr»=Fr- cos 0+Gr-t sin Or [i (0) 4-2 (6)] 
P= sin 0—Gr— cos 0 
Poo=Fr—! cos 0-+Gr~ sin 0. 

The term in G is of the same form as x, while that in F is the same 
deformation turned through a right angle, and so violates condition (b). 
To determine the form of H(0), we make a Fourier analysis of H’(0)+H(@). 
All terms in this analysis contribute to the self-energy of y1, while the 
only term which can reduce the energy by its interaction with G is 
the term in sin@. Thus H’(6)+H(@)=sin 0, which gives 

H(6)=—30 cos 0+ F"! cos 6+G" sin 6. 
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The terms in F! and G@! lead to zero stresses, and there remains the term 
y!=—4r0 cos@. This is the term with coefficient B in eqn. (1), which 
has already been rejected because it is not in equilibrium unless a force 
is applied at the origin. These three terms, together with H(#)=30 sin 0, 
are the only biharmonic functions of the form rH(6). 

Laszlé (1950) has suggested that the two dimensional dislocation 
discussed here may not be the form of lowest energy, since it contains 
a long cylinder of compressed region which might reduce its energy by 
buckling, while the energy of the surrounding material increases. The 
problem has not been studied quantitatively. 


2.2. THE SCREW DISLOCATION 


/ 
Figure 7 represents schematically the positive and negative dislocations 
of fig. 6. In fig. 7 (a) glide has started at T on one side of the crystal, and 
has not yet reached the other end T’ of the slip plane. There is a 


Fig. 7 


a 
— 


(a) (0) 
(a) Passage of a positive dislocation across the glide plane from T to T’. 
(b) A negative dislocation. 


region AA’ of the slip plane which contains an edge dislocation in which 
n atoms above the plane are associated with n—1 atoms below. In the 
ideal two dimensional edge dislocation of § 2.1, the same disturbance 
appears in all sections parallel to the plane of fig. 7. We must expect 
that in a real crystal this ideal two dimensional picture will not be valid, 
and that the region of slip TA will have progressed further across the 
slip plane in some sections than in others. Looking down on the glide 
plane of fig. 7, we shall see the dislocation schematically as in fig. 8. 
The parts of the dislocation represented by single lines are edge 
dislocations such as we have already considered. The parts represented 
by broken lines are screw dislocations of the type introduced into 
plasticity theory by Burgers. Whereas the line of an edge dislocation 
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is perpendicular to the slip direction, a screw dislocation is a line parallel 
to the slip direction which separates regions of the glide plane over 
which slip has already occurred from regions which have not yet slipped. 


Fig. 8 
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The glide plane of the crystal of fig. 7. The portion of crystal to the left of the 


boundary has slipped, while that to the right has not. Continuous lines 


in the boundary are segments of edge type, broken lines are segments of 
screw type dislocation. 


Fig. 9 
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A screw dislocation along the x axis (after J. M. Burgers). Continuous lines 
indicate rows of atoms above the zy plane, dotted lines indicate rows below 
this plane. 


Figure 9 shows the displacement of atoms above and below the slip 
plane in the neighbourhood of a screw dislocation. Although we have 
associated a screw dislocation with an edge dislocation lying in a particular 
glide plane, the screw dislocation can move in any plane passing through 
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its axis, and its sections by all of these planes are similar. For example, 
figs. 8 and 9 show the section of a screw dislocation by a plane passing 
through TT’ of figs. 7 and 8, and perpendicular to the plane of fig. 7. 
Suppose we now take a section parallel to that of fig. 7, and passing, 
not through TT’ of fig. 8, but through UOU’. Then (fig. 10) the section OU’ 
of the screw dislocation appears as in fig. 9. At O we have slight 
distortion of the lattice without any dislocation below the plane of the 
figure, and an edge dislocation in planes above that of the figure. This 
is an example of the general rule that a dislocation cannot end in the 
crystal. Its segments of edge, screw or mixed type must extend 
continuously from one boundary to the other, or else form closed loops 


Fig. 10 


The same dislocation as in fig. 9 lies along the positive x axis. At the origin 
the dislocation line bends into the direction Oz; the plane of atoms above 
the ay plane contains an edge dislocation (after J. M. Burgers). 


(cf. § 2.3). By a segment of mixed type is meant one which lies neither 
perpendicular nor parallel to the slip direction. Introducing these, we 
see that a dislocation need not form a rectangular zig-zag as in fig. 8, 
but may follow any smooth curve in the glide plane. The loop dislocation 
of fig. 4 is an example of a dislocation composed of segments of mixed type. 

The screw dislocation is less easy to visualize than the edge dislocation. 
It is perhaps most easily pictured by considering the ring of fig. 1 (a) to 
be a section of a long hollow cylinder. If this cylinder is cut by a plane 
such as C’C passing through its axis, and the surface above C’C slides 
over the surface below through one atomic distance in a direction parallel 
to the axis of the cylinder, the two surfaces may again be joined with 
almost perfect matching of the atomic lattices. If the same process is 
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tried with a solid cylinder, there is no difficulty in the neighbourhood 
of C’, but there is a singularity with infinite elastic strain at the centre O. 
This is a screw dislocation. 

Partly because of the difficulty of visualizing a three dimensional 
deformation, and partly because the original papers of Taylor and 
Polanyi considered only edge dislocations, little consideration was given 
to the effects of screw dislocations until about 1948. In fact, the screw 
dislocation, which can travel in any plane in which it lies, plays at least 
as important a part in plasticity as the edge dislocation, which can travel 
only in the direction of the displacement it produces. 

The mathematical theory of screw dislocations is simpler than that of 
edge dislocations, for only pure shear strains are involved, without 
dilatation. For a screw dislocation lying along the Oz axis, the only 
component of displacement is u,, and the displacement depends only 
on x and y, being given by 


b 
Ue ope PO Ne ein Bete oe a eng) (LT) 


We use tensorial components of strain e;, defined by 
_ Ou; , Ou; 
es Ox, dx; ’ 


and only components of strain which do not vanish are 


2e (12) 
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(13) 
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The corresponding stresses are 


(14) 


Tn cylindrical coordinates these become 
U,=b0/2a0 
€9,=b/4ar tt ad RSS Oy Bcak Ae (15) 


: Pose] 2nr, 
and all other components vanish. It will be seen that in this elastic 
solution the displacement is independent of r. 

As in the case of edge dislocations, we may show that of all displacements 
which have only u, components, and which depend only on @ and y, 
only that given by (11) has the properties of a dislocation. For suppose 
that 


Up=U,=9 
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Then the only stresses are 
Pne=HOU/@x 
and Pye LOU Oy. 


The equations of equilibrium in the x and y directions are automatically 
satisfied, and the condition for equilibrium in the z direction is 


OD ey Og 

Sr Nee 

aU UAE 
or Oe Oy =. 


In cylindrical coordinates, 


Pre=PO U/or (16) 
Po,=p0U|/rd8, 
Fig. 11 
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A dislocation line AB which ends in the medium at B may be surrounded by 
a cap PQR. Outside this cap the strain is continuous and can be 
differentiated at least twice. The circuit MjM encircles the cap. 


and, as in § 2.1, we need only consider harmonic functions U for which 
the stresses are proportional to 1/r at large distances. The only function 
of this kind is H+-F Inr+Gé@. The constant H leads to zero stress, and 
U=Gé@ is the screw dislocation just discussed. FInr gives a single 
valued displacement which does not interact with G@, and a combination 
of these terms has a higher energy than the term G@ alone. 


2.3. A DISLOCATION CANNOT END INSIDE THE MEDIUM 
We shall first give a formal proof of this theorem by means of the 


theory of elasticity, and then illustrate its meaning by considering a 
dislocation in a crystal lattice. Suppose AB in fig. 11 represents a 
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dislocation which ends in the medium at B. We take (cf. Love, p. 222) 
a point M, close to the line of the dislocation, write down the 
components u; of the displacement, and evaluate the integral 


a ae eo (17) 


for a closed circuit starting at M, and encircling AB. Then if AB is a 
dislocation, the integral (17) does not vanish, and its value is the strength 
of the dislocation, a vector b. Since the line of the dislocation is likely 
to be a singularity of the strain, we surround it with a surface PQR and 
take the circuit outside this surface, where the strain is continuous and 
can be differentiated at least twice. Though the displacement may not 
be single valued, the strain, which determines the physical state of the 
medium, must be so. 
Then 
Ou; 


du; 
as ds= Be, dx, 
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The repeated suffix 7 in these formulae implies a summation over all 
values of 7. Integrating from M, to an arbitrary point M on the circuit, 
we obtain 


ei (ec ean) 


M du; Ou; 
= [eis dey eM) 2 (M)} |e! Me) — 5! (Me) 


Ou, du; 
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du; du; a 1 
if ie — aa has more than one value, any one of these values may 
J 


be selected, and inserted in both terms of (19), so that we follow the same 
branch of the many valued displacement functions. 


Now 
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and so, changing the dummy index in the first integral on the right of 
eqn. (19) from j to k, we have 


eM 


a [exact fey 0M} (2 — Fit) | ae OB 


If the circuit is closed, so that M coincides with Mp, the first term on 
the right of (21) vanishes, while the second may be written [é;, dv, 


Oe, Olp; we 
where E n=l in t+ {x ;—2 (Mp) } (Za - ) iy, fe pied eee 
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Since the dislocation ends within the medium, we may replace the 
integral of €;, round the circuit by the integral of its curl X;, over a 
surface which has the circuit for an edge, and lies entirely within the 
medium and outside the surface PQR. ‘The value of X,, is given by 


OE sin 
Xin=€him Be? i; 8 Ce go Gg hee eee 
where €;,,, is +1 if (klm) is an even permutation of (123), —1 if (kim) is 
an odd permutation of (123), and zero otherwise. If this expression is 
evaluated, it is found that X,, vanishes on account of the conditions of 
compatibility of elastic strains. Thus (21) vanishes for a closed circuit, 
and the strength of the dislocation is zero. 

To visualize what happens if we try to end a dislocation in a crystal 
lattice, consider the dislocated plane of fig. 6(b). Take any closed 
circuit surrounding the dislocation once anticlockwise, and count the 
number of lattice spacings up, down, left and right in the circuit. The 
number of steps up is equal to the number of steps down, but the number 
left always exceeds the number right by unity. This residue, which is 
analogous to the integral (17), is a measure of the strength of the 
dislocation. Now consider a large circuit of this kind, far from 
the dislocation, so that the atomic spacings on this circuit are nearly 
uniform. Build up the next sheet of atoms by putting in the atoms 
which lie immediately above this circuit and then filling in the closed 
ring of atoms as well as possible. Usually, the atoms in this closed 
ring will lie in one-to-one correspondence with the atoms of the circuit 
below. In this case the excess of steps left over steps right is unity for 
the new circuit also, and the new plane must contain a dislocation within 
the circuit. If there is not a one-to-one correspondence, a line of 
dislocation must cross from inside the boundary to outside at the point 


where this one-to-one correspondence breaks down. In either case the 
theorem is verified. 


Mathematical Theory of Stationary Dislocations 287 


§3. Tor Format Exastic Tarory or Distocations 


b] 
3.1. WEINGARTEN’S THEOREM ON THE POSSIBLE TYPES OF DISLOCATION 


We consider multiply connected bodies of two types. The first is 
exemplified by a torus, such as the ring of fig. 1; the second may be 
represented by a sphere with a toroidal hole inside. Each of these is 
doubly connected, and has the property that a closed circuit can be 
described in the medium, and this circuit cannot shrink to a point while 
always remaining inside the medium. In either case a single cut removes 
this property. For the torus, a suitable cut is that made by an axial 
plane, such as AA’ in fig. 1 and the resulting body is topologically a 

, Sphere. For the toroidal hole, the cut is made outwards from the hole 
to the outer surface (fig. 12) and the result is again a topological sphere. 


Fig. 12 
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A body with a toroidal hole of small bore is made singly connected by a cut. 
M, and M, are neighbouring points on either side of the cut, N, and N, 
are a second pair of points separated by the same cut. 


In either case the body may be dislocated by making a cut (which need 
not be plane), displacing the two cut surfaces and possibly adding or 
removing thin layers of the medium, and finally joining the surfaces 
again in their disturbed positions. 

Weingarten (1901) has shown that if the strains in the deformed body 
are everywhere finite and twice differentiable, the disturbance of the cut 
faces cannot be arbitrary, but their relative displacement must be one 
which is a possible displacement for a rigid body. The proof follows 
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(Love, p. 223) from the discussion of § 2.3. We consider the body of 
fig. 12. The surface PQR of fig. 11 is now replaced by the toroidal hole, 
and the circuit is one which cannot shrink to a point without leaving 
the medium. We take M, to be a point lying close to the cut surface, 
and M, the corresponding point on the other side of the cut. Similarly 
N, is a different point on the same side of the cut as Mo, and N, is the 
neighbouring point on the other side of the cut. We join M, to No 
by a path lying in the medium and close to the cut, and M, to N, by a 
neighbouring path on the other side of the cut. From eqn. (21) 


du, du; 
w (Ny)—u(My)—=—B (Ny) —2()} {52 CM) — 5 M)} 
4. ( £e de, | 04s 
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Now in (24) and (25) 


w (M,)=2 (Mp) 


x (N,)=x,(No), 


the paths of the two integrals coincide, and €, has the same values at 
corresponding points of the two paths because the strain components 
are continuous. Subtracting (24) from (25), we obtain 


uN 1)—u(No)=[u,(My)—u (Mp) ]—3 {@(N)—ax,(M) } 
du; Ou; Ou; Ou; 

[geld — Fe ary} — {Fay Say} |. 26) 

The expressions in square brackets in (26) are independent of the 

coordinates of N, and we see that the relative displacement 5u,(N) of 

neighbouring points on either side of the cut at N is given by an expression 
of the form 

5u(N)=6,+4,,0,(N); 9°. . 9) ow we dee 


in which d;;=—d,,._ This represents an arbitrary rigid-body displacement. 
The theorem, then, shows that in a multiply connected body we may 
make a cut, not necessarily plane, which renders the body singly connected. 
If the body is free from strain, the cut surfaces will remain in equilibrium 
in their original positions. Keeping one surface fixed, we may give the 
other any small rigid displacement, remove any surplus material where 
the surfaces overlap and fill in any gaps, seal the cut again, and release 
the applied stresses. Then the strains in the resulting body are continuous. 
If the second surface is deformed in any way which does not reduce to a 
rigid displacement, the strains are not finite and twice differentiable. 
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This is a very remarkable result, and appears to contradict intuitive 
ideas that a lenticular hole could be excavated in an elastic solid and the 
faces then welded together without difficulty. The theorem holds only 
on the assumption that the strain is continuous and twice differentiable. 
As an example of an apparent exception to the theorems of §§ 2.3 and 3.1, 
we may consider again the ring of fig. 1. If the ring is thin, so that AA’ 
is much less than OA, the radial displacement is clearly given to a good 
approximation by w,—rd/27, and it is not immediately obvious why 
the dislocation shown can only be effected in a multiply connected ring 
and not in a solid cylinder. The exact solution (Love, p. 227), shows 
that w, contains a term proportional to rlnr. The strain is therefore 
infinite at the origin. 

The result has been proved for the doubly connected body and the 
type of cut shown in fig. 12, but it is equally valid for the torus of fig. 1, 
or for the body of fig. 12 if a cut is made inwards from the toroidal hole, 
so that the cut body has the connectivity of a hollow sphere. If the 
toroidal hole is of very small bore, the cut is plane, and the relative 
displacement of the cut surfaces is one atomic distance in a direction lying 
in their plane, the resulting strain is that produced by Orowan’s loop 
dislocation (fig. 4). 

In the present mode of description, the dislocation of fig. 1 was made 
by cutting along AA’, keeping the left face of the cut fixed, and rotating 
the right face rigidly about O until the line BB’ fell on AA’. The wedge 
AA’ B’B was then removed. It was clear that the ring could be brought 
into the same state of strain by making any radial cut such as CC’. 
Volterra (1907) has proved that any dislocation of Weingarten’s type is 
equivalent to the dislocation produced by applying the same translation 
and rotation to the surfaces of any cut which can be continuously 
deformed into the original cut. 


3.2, THE SIX SIMPLE TYPES OF DISLOCATION 

According to eqn. (27), the most general dislocation of a doubly 
connected body is described by three constants b; specifying the relative 
translation of the two cut surfaces, and three constants d;,; specifying 
their relative rotation. Each of these six may be regarded as an independent 
dislocation, and the corresponding constant is the strength of this elementary 
dislocation. For a closed ring (not necessarily a body of rotation) encircling 
the z axis, displacements which will maintain equilibrium under no body 
forces, a. surface forces are ee are given by 


QU p=(Dp—Fog? +> pay) tan! (y/x)+4(—b,—d,,2—ad,,,x) In (x*+y”) 
Qn, =(by—d,,t+d,,2) tan (y/x)+ ” —d,gt— adn yy) In (%+-y") (28) 
2rru=(b,—dy+d_%) tan™ (y/a)—+2 ouane, In (2?+-y"), 
where «=p/(A+2y). 
‘Here A is Lamé’s elastic constant, related to the rigidity 4. and Poisson’s 


ratio o by the expression 
: A= 2op/(1—2¢). at me se pages. (29) 
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Volterra has discussed in detail the disturbance produced by each of 
these six types of dislocation in a hollow circular cylinder. If we take 
the origin of coordinates at the centre of the cylinder, the axis of the 
cylinder as z axis, and the positive x axis lying in the cut, the six orders 
of dislocation are those sketched in fig. 13. 

The first order, b,, is the edge dislocation. The difference between 
the displacements given by (28) and those given by (9) is caused by the 
neglect of surface forces in the derivation of (28). The second order, by; 
is identical with an edge dislocation along the y axis, as may be seen either 


Fig. 13 


The six simple types of dislocation. (a) Original cut cylinder.- (6) Dislocation 
of type bz. (c) Type —b,. (d) Type b,. (e) Type —d,,. (f) Type —dzz. 
(9) Type day. 


from the form of eqns. (28) or by considering the cylinder cut and 
displaced along the dotted line of fig. 13 (c) instead of the full line of 
fig. 13 (a). This is an example of the theorem that two dislocations are 
equivalent if they represent the same rigid-body displacement of the 
faces of cuts which can be continuously deformed into each other. The 
third order, b,, is the screw dislocation of § 2.2. The components of 
displacement in eqns. (28) and (11) agree. In the last three orders, the 
residue of the displacement integrated round a circular circuit isnot 
independent of the radius of the circle as in the first three orders, but 
proportional to the radius, and the strain does not tend to zero at infinity. 
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We are therefore not concerned with these dislocations in the theory of 
plastic flow. Here again the dislocations of the fourth and fifth orders 
are equivalent. The dislocation of the sixth order is that shown in fig. 1. 


3.3. MORE GENERAL DISLOCATIONS 


Mann (1949) and Bogdanoff (1950) have recently analysed dislocations 
of a more general type than those discussed by Weingarten and Volterra. 
They allow the tangential components of stress to be discontinuous across 
a plane cut, while the tractions across the cut are still continuous. The 
resulting stress system is of the type occurring in Eshelby’s analysis of 
the core of an edge dislocation, or in a slip band in which dislocations 
are still present in the slip plane. It may be regarded as arising from 
the superposition of a continuous distribution of simple dislocations. 
Mann considered stress functions of the form In r or 6 multiplied by 


IIE ee Be | TMA hee Pees (0) 


Bogdanoff showed that the general solution requires also three stress 
functions which were disclosed to him by R. D. Mindlin, namely 


CGTHE SRE Woke atone bee tr ENS, 


and the same factor multiplied by 7 cos @ and r sin @. 

Eshelby (1951) has considered the case in which the medium contains 
a continuous distribution of dislocations. The elastic conditions of 
compatibility may be regarded (Love, p. 50) as expressing the vanishing 
of a tensor 


S,5=Cps, tale treme Cre ots Coh ria a3, bOra ea atte ete Ae) 
If we consider a non-Euclidean space defined by the metric 
9 499 53 +2653; ° : . . . . . . (33) 
its Riemann tensor is given to first order in e;,; by 
Poarst=€ ipr€ ist ix . . . . . . ° ° (34) 


The vanishing of S;; is thus the condition that the space (33) should be 
flat. The diagonal terms of S;; represent the distribution of infinitesimal 
dislocations of the type of fig. 4, the non-diagonal terms the distribution 
of dislocations of the type of fig. 15. 


3.4. BURGERS’S THEORY OF LINE DISLOCATIONS 


Burgers (1939 a) has developed Volterra’s theory of dislocations 
for the case which is of importance for the theory of slip, that is, 
dislocations of the first three orders in solids from which filaments of 
matter have been removed in such a way that the radius of curvature of 
the filament is everywhere large compared with its diameter. He 
considered a closed filament, corresponding to the toroidal hole of 
fig. 12, in an isotropic medium. A dislocation is made by a relative 
displacement b of the two surfaces of a cap whose edge is formed by this 
loop. The loop is known as the line of the dislocation : it is convenient 
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to have a name for the vector b, and we shall call it the Burgers vector 
of the dislocation. The strain at any point depends only on b and on 
the position of the line of the dislocation, and is independent of the form 
of the cap. Burgers illustrates this fact by expressing the displacement u 
at a point x by an integral round the line of the dislocation, together with 
a multiple valued term proportional to the hydrodynamic potential 
at x of a vortex of unit strength along the line of the dislocation. The 
complete solution is given by 
u=u*-+u**- orad ws. Oye aoe 
Here u*— Qb/47, 
where Q is the solid angle (a multiple valued quantity with residue 47) 
subtended by the dislocation line at x, 
b dé 
wh ey leas 
ais Aor | r 
where d& is an element of the line of the dislocation, r= |x—§|, and the 
integral is taken once round the dislocation line in an arbitrary sense, and 
A+p dé 
os oe ale ar r? 
In this solution, the force and acting couple on unit length of the 
filament are of the same order as its radius, and so are negligible. 


> 


Fig. 14 
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ds 


Sign convention for the Burgers vector. The positive direction on the line of 
the dislocation d€ is chosen arbitrarily. The Burgers vector is minus the 
residue of displacement round a right handed circuit f ds. 


The signs of Taylor’s two dimensional dislocations can be established 
by one of the usual conventions of contour integrals. Orowan’s loop 
dislocation appears in section as a pair of edge dislocations of opposite 
sign, and a convention is required to define the sign of its Burgers vector. 
We adopt the convention that the positive direction dé of the circuit 
of the dislocation line (fig. 14) is chosen arbitrarily. 
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Then the Burgers vector is given by b= — J(du/ds) ds, where the integral 
is taken once round a curve which encircles the dislocation line with a 
right-handed screw. Bilby (1951) has discussed this rule. 

Burgers (1939 b) has extended his results to the case of a dislocation 
loop in a cubic crystal. The general formulae are complicated, but in 
the case of screw disiocations with their lines along a cube axis the 
displacements are the same as those in an isotropic medium, while there 
is also a considerable simplification for edge dislocations along a cube axis. 

Eshelby (1949 a) has obtained the solution for an infinite straight edge 
dislocation in a crystal of any symmetry. His general formulae are 
explicit, but necessarily rather clumsy. They simplify in a number of 
important cases for which, owing either to the symmetry of the crystal 
or to the special direction of the dislocation line, the ellipses of stress and 
strain in the plane of the strain have their axes parallel. He has applied 
his result to estimate the energies of edge dislocations lying along various 
lines (cf. § 2.1 of Part III), and their widths and mobilities (cf. § 5.3 of 
Part ITI). 

Nabarro (1951 b) has evaluated Burgers’s formulae for an isotropic 
medium for the special case of an infinitesimal loop of the type of fig. 4. 
Since (§ 4.1) any glide motion may be considered to be produced by the 
creation of such infinitesimal loops, the results may be used to synthesize 
the stress field of any dislocation which could have been caused by glide. 
The method does not seem to be of much value for stationary dislocations, 
but it is the only general method available for moving dislocations. 


3.5. THE ENERGIES OF ELASTIC DISLOCATIONS 


We consider an isolated elastic body containing a number of holes, 
each of the type shown in fig. 12. The body is rendered singly connected 
by a set of surfaces 2, each bounded by one hole, and dislocations 
b;,, dj, (eqn. (27)) are made on each surface. Suppose that when the 
dislocations are made, the force acting on an element d2’, of that surface 


of the cut 2, which suffers the displacement 
Ujs= Di + dys; 


has components X,d2,. Then the energy EH of the system is easily 


seen to be 
H=4 $f (b,,+diyt;)X; dd’, Me ti eee es. 04 (SO) 


where summations over all the repeated suffixes are implied., This 
formula may be employed to calculate the energy of interaction of two 
dislocations. Since b;, and dy, are constants for each surface, they may 
be taken outside the integrals. Each integral then represents a component 
of the resultant force or couple transmitted across a surface ,. 
Alternatively, the total energy may be calculated by integration if the 
components of strain and thus the energy density are known at all points 


of the body. 
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It has already been observed that the energy for unit length of an edge 
or a screw dislocation in an infinite medium diverges logarithmically at 
infinity. In practical applications this difficulty is usually overcome by 
cutting off the integration at a distance which may be determined by the 
average separation of dislocations or some similar consideration. For a 
dislocation loop. the difficulty does not arise, as the integrals in (36) are 
clearly finite. Alternatively, we observe that, at distances large compared 
with the dimensions of the loop, the displacement corresponds to that 
caused by a ‘ double force ’ (Love, p. 187), and so is of order r*. The 
strain is of order 7~*, and the energy is finite. 

Suppose now that the multiply connected body is subjected to external 
forces, which are held constant while the dislocations are formed. If 
small dislocations are formed, the forces acting on d2’, are the sum of 
those X,dX’, produced by the dislocations and the constant X ,’ dd, 
produced by the external forces. These X,’ are derived from the actual 
stresses over Y,, so that even if the stress in the body is uniform far from 
the dislocations, it will be non-uniform over 2’, because of the stress 
concentration produced by the hole surrounding 2. Work is also done 
by the external forces when the dislocations are made, because the outer 
boundary of the body deforms. Let this work be W,. Both W, and the 
work done by the forces acting across 2’, when the dislocation is made 
are linear in the strength of the dislocations, while the self energy HE of 
the system of dislocations is (eqn. (36)) quadratic in their strength, and 
may be neglected. The work done in making the small dislocations 
is thus 

W ASS(O;5 +4j;,0;)X ; dX, = ais Dee Be 
where 6,,’ and d,,,’ are the strengths of the dislocations, and X,; are the 
forces acting over 2’, due to the external forces. 

If ¢ is the strain energy density at any point, ¢ is a function of the 
strains e;; produced by the external forces, and the work done in making 
a small dislocation may be written 


[Jf xe eu’ ar. 2 wl gh ol ee 


Here e;,;’ are the strains produced by the dislocations, and the integral 
is extended over the volume of the body. 

Now consider the same body under no external forces, but with 
dislocations already formed. If arbitrary small strains e;; which do not 
form dislocations are imposed on the body, no work is done by external 
forces or by forces over the surfaces X,, and so by the principle of virtual 
work there is no change in the stored energy. The strain is now Cr 
and we find 


[[fxpewar=o. he Mey ey oe 


In particular, we may take the e,; in (39) proportional to the strains 
produced by the external forces in (38). The function ¢ is homogeneous 


— 
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and quadratic in the strains, and the integral in (39) may be transformed 


by the relation 
0d’ : 
eee) Cry, af 2 es pe Crue 
) ij rs rs 
It follows that eqn. (39) implies the vanishing of expression (38). 
Expression (37) must also vanish, and the work done by the external 
forces when the dislocation is made is 
W = ae ssf (Djs +d 555% 5)X ; dX,. 

This important theorem is due to Colonnetti (1915), and is of value in 
classifying ‘the motions of dislocations (§ 4.1) and in the study of the 
dynamics of dislocations. For the dislocations of the first three orders 
which are important in plasticity we may write 

Vat LO aD poi Piiat ep eden cide Os x con kA) 
where p;; is the stress produced by the external forces, and dZ;, is a 
component of d,, regarded as a vector. The sign of the vector is positive 
if it is derived from the circuit [dé of fig. 14 by a right handed screw rule. 

Following a suggestion of Leibfried (1949), we may interpret 
Colonnetti’s theorem as follows. The vanishing of (37) means that the 
work done in forming dislocations is independent of the external stresses. 
If we express the total strain energy as a quadratic function of the 
strengths of the dislocations and of the applied stresses, there are no 
cross terms between dislocations and applied stresses. This means 
that as long as the dislocations remain constant and Hooke’s Law 
remains valid, we cannot detect the presence of dislocations by observing 
the response of the body to external forces. 


i) 


§4, THe Motion or DISLOCATIONS 


4.1. POSSIBLE MOTIONS OF A DISLOCATION 


The possible motions of a dislocation may be classified (Nabarro 1951 a) 
by means of formula (40). Suppose a single dislocation is formed in a 
solid under a uniform pressure P, so that the stress is given by p;;=—6 ;;P. 
Here §,,=1 if t=j, and 6,,;=0 otherwise. This is the stress at all points 
of X, except those within a distance of order 6 of its boundary. The 
energy of interaction in this boundary region is of order 6°PL, where L 
is the length of the dislocation, and may be neglected if the area of the 
dislocation is large compared with 6?. For a large dislocation loop, the 
work done by the external pressure P when the loop is formed is therefore 


W =P jj 6; d2;, 
See Bid cr ter en gated ec fa, ge (AL) 
and {{b dz is the increase in volume of the body when the dislocation is 
formed. Since (Love, p. 175) in the absence of surface or body forces 


the mean value of the dilatation in the dislocated body is zero, 
f{bd=Z must be the volume of matter inserted into the body when the 
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dislocation was formed. For the loop dislocation of fig. 4, b and d= are 
orthogonal, and f{fbd=Z=0. If an extra layer of atoms of thickness 6 
and having the same boundary had been inserted into the lattice (fig. 15), 
the resulting dislocation would be characterized by [J bd2Z—b x (area of 
loop), while for an equal number of atoms removed [{bd2Z would be 
negative. 

In normal glide motions the dislocation line always moves on a 
cylindrical surface with the vector b as axis, and both b and the position 
and shape of the projection of X on a plane perpendicular to b are 
conserved independently. However (Frank 1951 a), this projection 
may add segments of zero area to its periphery. For example, the 
dislocation of fig. 15 is everywhere of edge type, and each part of it can 


Fig. 15 


Section of a loop of dislocation with Burgers vector b normal to the plane of 
the loop (ay plane). 


only move in the z direction. The parts which have moved are now 
joined to the parts which have not moved by segments of screw type, 
and these can glide freely over any cylindrical surface with b as axis. 
The dislocation is shown in perspective as ABCDA in fig. 16 (a), and its 
projection on the ay plane in fig. 16(b). Its movement to the form 
ABEFDA does not alter its projection. In fig. 16(c) the screw 
segment BE has sent out-an open loop BGHJE, and the xy projection 
is now that shown in fig. 16 (d). The loop GHJ, which projects into 
a single line in the plane perpendicular to b, is similar to a segment of 
the dislocation of fig. 4. 

It is possible for the value of [f b dZ to be conserved without preserving 
the position and shape of the projection of the dislocation line 
perpendicular to b. For example, the added sheet of atoms in figs. 15 
and 16 can move bodily in the ay plane, first gliding over the plane of 
atoms below and then allowing the plane of atoms above to glide back. 
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The process does not use the existing dislocation b with its vector in 
the z direction, but involves the formation of new loops of dislocation 
with Burgers vectors in the wy plane, for which {fbd=Z=0. If only 
a part of the area ABCDA glides in this manner, the shape as well as the 
position of 2 alters. It is also geometrically possible for part of S to 
glide over another part of X, producing a region with Burgers vector 2b. 


Fig. 16 


2 
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Possible glide motions of a dislocation. (a) Loop of dislocation ABCDA, lying 
on a cylinder with axis parallel to oe Burgers vector b. The segment 
BCD can glide to EF, leaving screw segments BE, DF which join it to 
DAB. (6) Projection of (a) normal to the Burgers vector. (c) A loop GHJ 
can spread from the screw segment BE, provided that it lies on a 
cylindrical surface with b as generator. (d) ‘Projection of (c) perpendicular 


to b. 


These movements are not of direct importance in plastic deformation, 
because they do not do work in a uniformly stressed body. In normal 
glide motions the projections of 2 on planes passing through b alter, 
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and so, from eqn. (40), work is done by shear stresses (ij), but in these 
motions all the projections of Y are unaltered. A rigid motion of the 
curve ABCDA up its cylinder is also of this ineffective kind, although 
it does not require the formation of new loops of dislocation. 

The third possible type of motion is a non-conservative motion in 
which the value of {fbdZ alters. This involves the movement of 
segments of edge type in a direction perpendicular to b, and occurs by 
the diffusion of vacant lattice sites or of interstitial ions to the line of 
the dislocations (Nabarro 1948 b). When one vacant site has reached 
the edge of the dislocation in each lattice plane, the dislocation line has 
moved by one lattice spacing perpendicular to b, as may be seen from 
figs. 6and 15. Such a process is slow, and its rate depends on temperature 
in the same way as the coefficient of self diffusion. The opening of a 
crack in a body may also be regarded as a non-conservative motion of 
a dislocation line. Non-conservative processes, which introduce an 
important element of irreversibility into plastic deformation, have been 
discussed in detail by Seitz (1950 b, 1952) and Mott (1952 a). 


4.2. THE FRANK—READ SOURCE 


Frank and Read (1950) have introduced a type of motion of a 
dislocation which has profoundly influenced the development of the 
theory of slip. Previously attention had been concentrated on 
dislocations lying entirely in one plane. The growth of such a loop in 
the glide plane causes slip, but at the same time the dislocation passes 
out of the region under consideration and ultimately passes right out 
of the specimen at its surface. If two portions of the loop happen to 
meet, they annihilate one another. There is thus no explanation for 
the large amount of shear which is experimentally observed to occur 
on a single glide plane. Frank (1948) showed that fast moving 
dislocations might multiply when they struck the free surface or other 
fast dislocations, but it seems clear (Leibfried 1950, Nabarro 1951 ce) 
that such processes can occur only in the upper part of the range of 
stresses which in practice are applied to metals. 

Consider now the loop of dislocation ABCD of fig. 16, which does not 
lie in one glide plane, and suppose that owing to the details of the crystal 
structure, e.g. the possibility of dissociation in certain planes (§ 5.2 of 
Part II) the loop forms a rectangle with corners ABCD. Suppose 
moreover that the segments AB and CD are anchored, while glide can 
occur in the planes containing the slip vector b and one of the lines 
BC, DA. Then each of the segments BC and DA can move in its glide 
plane without interference from the other segments. Figure 17 
illustrates what happens in the glide plane of BC when a stress is applied. 
The original line | bulges to a loop 2, then to 8. At 4, two elements of the 
loop are about to meet. They annihilate, and the next configuration is 
a closed loop 5 together with a short segment of dislocation joining 
Band C. The loop now expands freely, and is ultimately lost, while the 
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short segment again starts the sequence 2, 3, 4... The original 
segment 1, which can give rise to an unlimited succession of closed loops 
of dislocation, is called a Frank—Read source. 


4.3. THE FORCE ON A DISLOCATION : LAW OF CONSTANT RESOLVED SHEAR 
STRESS 


Taylor’s expressions for the stresses round an edge dislocation in an 
elastic continuum require a resultant force to act on any small closed 
curve surrounding the dislocation, whereas those given by Burgers do not, 
and for this reason we have rejected Taylor’s solution. Suppose now that 
a dislocation is in equilibrium in an elastic medium, and a shear stress is 
applied which is uniform in the parts of the medium far from the 
dislocation. There is, in the elastic solution, a concentration of stress 


/ 
Fig. 17 


The Frank—Read source. A segment BC of dislocation lies in the glide plane, 
which is the plane of the paper. At B and C it turns out of the plane, 
and is held fixed. Under an applied stress the mobile segment BC moves 
from 1 to 2, then to 3, 4 and finally to 5. It now consists of a large loop 
together with the original segment BC, which can again repeat the 
sequence 1, 2, 3.:... . 


near the filamentary hole surrounding the dislocation, but the resulting 
stress distribution does not require forces to be applied across the 
surface of this hole, and in this sense the dislocation is still in equilibrium, 
with no force acting on it. Ina crystal lattice, the line of the dislocation 
is no longer fixed, and the energy of the system may alter if the 
dislocation moves. It is now convenient to use the expression ‘ force 
acting on the dislocation’ in another sense. If the displacement of a 
segment of the dislocation through a small distance d with components d; 
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causes the energy of the whole system, including the source of the 
applied stresses, to decrease by an amount U, we say that the force 
acting on the segment has components Ff’; given by 
0U 
Fea ag ee 
There are three contributions to the quantity dU, which correspond 
to the changes in the self energy of the dislocation, in its energy of 
interaction with other dislocations, and in its energy of interaction with 
the applied stresses. The first two are given by —dH, where £ is given 
by eqn. (36), and the third is dW, (eqn. 40). In interpreting these 
formulae, X,, and p,; must be taken as the forces and stresses actually 
present on X,. Because 2X, is bounded in the elastic model by a cavity, 


Fig. 18 


When the segment of dislocation CD moves to C’D’, the total area swept out 
is equal to AA’B’B. 


and in the crystal lattice by a region where Hooke’s law does not hold, 
the values to be taken in calculating the interaction energies differ from 
those which would hold if there were no cavity and no dislocation 
surrounding 2’, and these values alter when the dislocation moves. 
As in the derivation of eqn. (41), we consider only large loops of dislocation, 
and restrict attention to displacements in which segments having small 
radii of curvature move rigidly. As an example consider the dislocation 
of fig. 18, partly bounded by two parallel straight lines. If the part CD 
of the boundary moves by a glide motion parallel to these lines into a 
new position C’D’, the area of the dislocation is increased by 52. Since 
this is a glide motion conserving [{bdZ, the Burgers vector b must lie 
in the plane of the figure. The shift has added two new segments 
AA’, BB’ to the boundary. Their self energy is of order b2ud, where d is 
the distance AA’, assumed much less than the perpendicular distance D 
between the lines. Owing to the increased length of the boundaries 
which cause stress concentrations, the energy of the applied stress P 
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changes by an amount of order 6?P?d/u, and the interaction energy 
of P and the stress field round the dislocation changes by 62Pd. Since 
P is always less than p, the first of these three terms is the largest. The 
dominant term in dW, is of order bPdD, and the change in self energy 
is negligible in comparison with this provided D/b>/P. For a soft single 
crystal u/P may be of the order 5.104, and the self energy of the 
dislocation is negligible only for loops with D/b> 105, that is, for loops 
10° atoms across. In harder materials this figure is reduced to a few 
hundreds of atoms. 

Along CD and C’D’ there is no change of self energy. The interaction 
energy of P and the dislocation stress is of order b?PD, but the 
contributions from CD and C’D’ are nearly equal, their difference being 
of order b?Pd and so negligible. Provided that the stresses would in the 
absence of this dislocation have been constant on the surface 2, we may 
consider the added area to be the equal rectangular area AA’B’B and 
substitute in eqn. (40) the unperturbed values of p;; on this rectangle. 
In the general case of a change in the boundary of a small loop we shall 
again take the unperturbed values of p,,, but we take the actual new area 
62 near the boundary. 

If an element d& of the dislocation (fig. 14) moves a distance dx the 
area of the dislocation increases by 


OPENING ee a aree  P Peres od 44a) 
We write (40) in the form 
SW=bpPot 
=bp . (dx rd&) 
SID) ACS) OX ae es ea) (44) 
The force dF on the element d& is thus given by 
; dF=(b¥P) adé&. Mentioey tas eed) ot»). (40) 


This result was first obtained by Peach and Koehler (1950). 

In the definition of elastic dislocations from which this formula was 
derived, the creation or removal of matter was permitted, and (45) 
therefore gives the total force on the element, including the component 
which would lead to a non-conservative motion. To obtain the force f 
leading to glide, we subtract the component of F normal to the glide 
plane. The normal to the glide plane is the unit vector n defined by 


N= AACS) D AGE. 6. re) fe are pe) (46) 


The component df of the force dF which is effective in producing slip 
is the component lying in the glide plane, given by 


df=dF—n(ndF)=—na(nadF). . . . . . (47) 
Now 
nadF=na[(b xP) dé] 
= (bx)(nd&)—déln(bP)] 
= —d&(bjpn), 
and so df=(nad&)(b Pn). CCEA Os et renee (48) 


X 2 
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Since #n is the force acting across unit area of the glide plane, the 
force per unit length tending to produce glide is 6 times the component 
in the glide direction of the force acting across unit area of the glide 
plane, and acts in the glide plane normal to the line of the dislocation. 
If a loop of dislocation lies entirely in the glide plane, its conservative 
motions are confined to this plane (except for those elements which are 
of pure screw type), and its whole length is subjected to the same normal 
force tending to increase the area of the loop. This is in accordance with 
the experimental result (Schmid 1924) that slip begins when the ‘ resolved 
shear stress ’ b33n/b reaches a value which is characteristic of the crystal, 
but independent of the other components of the applied stress. 

The expression (48) differs from that given by Koehler (1941) by a 
numerical factor. Koehler calculated the work done by a uniform shear 
stress when the dislocation moves by integrating the strain energy within 
a large region, and then slightly displacing the boundary of the region 
over which the integral is taken. He did not allow for the difference 
between the displacements produced by the dislocation at a free boundary 
and the displacements at the corresponding points of an infinite medium 
(cf. § 4.4). The correct result was stated for a special case by Mott and 
Nabarro (1948) and proved by Leibfried (1949) and Read and Shockley 
(1950). When one considers the force between two dislocations or 
between a dislocation and an inhomogeneity of a different kind, the 
chances of error and confusion are enhanced. The position has been 
summarized by Eshelby (1951) and Peach (1951). 


4.4. EXPLICIT CALCULATION OF THE FORCE ON A DISLOCATION 


Eshelby (unpublished) has given detailed calculations of the forces on 
an edge dislocation lying along the axis of a cylinder of infinite length 
and finite radius, and on a screw dislocation lying along any line parallel 
to the axis of such a cylinder. Nabarro (1951 b) has considered the 
work done in forming an infinitesimal loop of dislocation at the centre 
of a homogeneously stressed sphere. The results are all in agreement 
with (48). 

4.41. Force on an Edge Dislocation 

The displacements produced by an edge dislocation lying in an infinite 
block along the z axis are given by (8), and the corresponding stress 
function x° is given by (1) and (4). In polar coordinates the stresses are 
given by 


x°=— Dr Inr sin 0 ee ee ty 

P+, =— Dr sin 0 
P,e°=Dr- cos 8 NS 6 ia oa nee ee 

and Pop’ = — Dr sin 6 
where D=bp/2n(1—o). Seat. Su A ee 


The tractions across a cylindrical surface of radius R with the 
dislocation as its axis do not vanish, and if this cylinder is cut out of the 
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block and its surfaces are freed from stress the displacements alter. The 
new displacements are derived from the stress-function y°+ x1, where 
Vacs iy oy BIO ee eee? s,s (BR) 
Over the surface r= the tractions produced by x1 annul those produced 
by x°. 
If now the dislocation is moved to the point x=d£, y=0, where 
d€<R, the stress function y° changes to x°+ 7, where 
x?= — SE 0° /0x 
=1D8é sin 20. ie te sates ica Pees (5 } 
The corresponding stresses are 
Pr = —2Dd5Er- sin 26 
Pre =Dsdsér- cos 26. 
The stress function which annuls the tractions over r=R is 
2=— DSER~r? sin 20+ DdER-4r4 sin 20, . . . (55) 
and the stress function for the cylinder with free surfaces and an edge 
dislocation slightly displaced from its axis is y°+y!+y’?+ y%. The 
displacement of the surface when the dislocation moves from (0, 0) to 
(5€, 0) is the displacement due to y*+ y3, which can be shown to be 
i= VOeein 20 mi withg—= Oe) oY 158) 
A shear stress p,,—P produces a normal force P sin 20 on unit area 
of the cylinder, and the work done in the displacement (56) is 


(54) 


li TE 
W.==> | | bP8E sin? 2. RdO=bPBE. 


The force on unit length of the dislocation is OW ,/0€=bP. 
4.42. Force on a Screw Dislocation 
The displacement and stress produced by a screw dislocation along 
the z axis are given by (15). For a parallel dislocation passing through 
the point x=€, y=0, the displacement is 
b 
uo= stan — z ere ee gsircal AO Td) 
while, on a circle of radius R>é€ about the origin, 
ES sin 6 
Pre~—~ D7 R2—2RE cos 0+ 2° 
The stress on the circle is annulled by a dislocation of the same 
strength and opposite sign at the inverse point (R*/é, 0), for which the 
displacement is 


(58) 


ui=— “a tan- (59) 


ip 
e Qa x—h?/E° 
Since on this circle 
R sin 6 


—2RE cos 0+ &’ 
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the total displacement of a free boundary when the dislocation moves 
from (€, 0) to (€+6€, 0) is 
bsé R sin @ 
bus — Re ORE cos 0 Le 
If a uniform shear stress P, for which p,,=P sin (@Q—«), acts across 
the surface of the cylinder, the work done in this displacement is 


(60) 


ee | " P sin (0—a)8u,Rd0=bPSE cos «. 


The force is bP cos «, again corresponding to (48), whatever the initial 
value of €. If the dislocation moves across the cylinder at a constant 
rate, the external forces work at a constant rate. There is no sudden ~ 
displacement as the dislocation enters or leaves the specimen. 


4.43. Force Expanding an Infinitesimal Loop of Dislocation 

The calculation is similar to that of §§ 4.41 and 4.42. The stress field 
of a loop in an infinite medium is obtained. To the corresponding 
displacements are added the displacements produced by forces which 
annul the tractions on the surface of a sphere of radius R large compared 
with the radius of the loop. If the loop is formed in the presence of a 
homogeneous external stress, the only component of this stress which 
does work is the component which gives a traction across the glide plane 
in the glide direction. 


4.5. FORCE ON A DISLOCATION IN A VARYING STRESS FIELD 


Suppose a dislocation is introduced into a body in which the stress 
field is not uniform. Then we wish to know if the force on each element 
of the dislocation is still given by (42), where $ is the applied stress in 
the neighbourhood of this element, or whether it also depends on the 
‘value of P at other points. The discussion of § 4.3 depends on the 
assumption that the stress is uniform, and if we calculate the energy 
of interaction of the dislocation stresses and the applied stresses we find 
that the contribution from regions far from the dislocation line are not 
small. Kshelby (unpublished) has verified that in certain special cases (48) 
gives the correct result, and has also given a general proof (1951). Here 
we consider two examples. 

The first verification is for the edge dislocation of § 4.41. If the 
stresses p,, and p,, on the boundary of the cylinder are expanded in 
Fourier series, it is clear from eqns. (56) that only that component of P 
corresponding to a uniform p,, shear stress does work when the 
dislocation moves. 

For the eccentric screw dislocation of § 4.42, the work done is given 
by (60) as 

~n ft” —2KE cos 0+ € 


where V(R, O)=aypoet YD yo. «| hee ree i0 ek 
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Now (§ 2.2), w, is a harmonic function in the zy plane. It follows 
that p,, and p,, are harmonic. V is not harmonic, but the quantity 
V=V+3(R?—a?—y’)p,, is harmonic, and is equal to V on the boundary 
a+y?= R?. 

We then have 

it bog i V'(R, 6) dé 
em J_, R?—2RE cos 6+ &’ 
which, by the use of Poisson’s integral theorem for a plane harmonic 
function, may be written 


W .=2bdEV'(é, 0)/(R?— &). 
Since V(&, 0)=0, this becomes 


“W=b8ép.,,(€, 0). eget ea (63) 


The force on the dislocation again depends only on the shear stress 
in its immediate neighbourhood. 


4.6. FORCES BETWEEN DISLOCATIONS 


The idea that dislocations exert forces on each other was introduced 
by Taylor, and forms the basis of his theory of work hardening (§ 5.31). 
He assumed that the force between two edge dislocations could be 
calculated from (48), where 38 was the stress produced by one at the 
centre of the other. Koehler (1941) calculated the force between two edge 
dislocations of opposite sign, not necessarily lying in the same glide 
plane, by integrating the strain energy density over the whole body. 
Since the dislocations are of opposite sign, the strains are proportional 
to r~ at large distances, and the energy is given by a convergent integral. 
His results differ from those of Taylor, and in particular he finds that the 
force is not proportional to the stress at the centre of each dislocation, 
although the numerical discrepancy is not very great. Cottrell (1949) 
considered a specimen in which one dislocation was already formed, and 
estimated the energy of interaction by considering the force required 
to make the second dislocation over a cut extending from its axis to a 
distant stress-free cylindrical surface with the first dislocation as its axis. 
The energy of interaction is the work done by the forces produced by the 
first dislocation acting through the displacements produced by the 
second dislocation. This method, which is essentially that of Volterra, 
has been extended by Bilby (1950a). Cottrell found the force of 
interaction to be the same as that which would be produced by a uniform 
stress equal to the stress at the centre of each dislocation. Eshelby 
(unpublished) has verified this result. ) 
- The error in Koehler’s calculation arises from the fact that the finite 
total energy is expressed as the difference between the infinite self energy 
of two dislocations and their infinite interaction energy. In evaluating 
these two energies the passage to infinity is made in the one case by 
considering a circle of large radius and in the other case by considering 
a strip of infinite length and large width. Eshelby has repeated the 
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calculation using bipolar coordinates, which allows him to avoid infinite 
energies. He has also treated the problem in a manner which involves 
no infinite integrals and no infinite or discontinuous integrands by a 
generalization of the treatment given by Dean and Wilson (1947). The 
problem is treated as one of plane strain. An infinite sheet contains 
two equal circular holes joined by a cut along their line of centres. By 
forces applied only to the surfaces of the cut, these surfaces are given 
equal and opposite parallel displacements in an arbitrary direction in 
the plane, and the work done by these forces in producing the dislocation 
is calculated. If we take axes along and perpendicular to the arbitrary 
direction, the calculated energy is seen to be that of a pair of edge 
dislocations in arbitrary positions with equal and opposite Burgers vectors. 
The result again agrees with that obtained by Cottrell. 

Since the total energy of a pair of dislocations of the same sign is 
infinite, it is difficult to calculate their interaction. We assume it is 
equal and opposite to the interaction of a pair of opposite sign. 


4.61. Force between Edge Dislocations 

The stress produced by an edge dislocation at the origin having a 
Burgers vector with Cartesian components (b, 0) is given by eqns. (50) 
and (51). The force on a dislocation of opposite sign at the point (r, 6) 
is given by eqn. (45), and has components 


F,.=—bDr 
Fy=—bDr— sin 20. | Me 
The interaction energy may be obtained by integration, and is 
W=bD In r—4bD cos 20, «th Sele eens 


where D=bp|27(1—o). 


If the second dislocation cannot move out of its glide plane (fig. 19), 
the component of force producing conservative motion is 


F,=—bDr™ cos 0 cos 26=—bp,,, . . . . . (66) 
or F=—}boDy™ sin 48. 43 30 ee 


The symbol 1 represents an end view of an edge dislocation. The long 
stroke is the glide plane, and the short stroke represents the tip of the 
half plane of atoms which ends on the glide plane. The same dislocation 
is occasionally represented by the symbol v. 

The glide plane is a plane y=const., and F,, is plotted as a faneuoe 
of x in fig. 20 (a). There is stable equilibrium for 6=7/4 and unstable 
for 6=0 and 6=7/2, with symmetrical positions on the other side of the 
y axis, 

For two dislocations of the same sign, these forces are reversed. The 
position of stable equilibrium is given by 6=4n and the position of 
unstable equilibrium by @=}7. In this case the position of stable 
equilibrium does not have the lowest energy of all positions which can 
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be reached by glide, for there is a long range repulsion between the 
dislocations. When 0=37, r=y, and when 6=47, r= /2y. Changing 
the signs in eqn. (65), we find 
W (0=4377) =—bD In y—4b6D 
and W (0=47) =—bD In v/2y. 
Fig. 19 
# 


glide plane 


Force on a dislocation at (r,@) produced by a dislocation of opposite sign at 
the origin. Only the component F’, is effective in causing glide. The 
dislocation at the origin would more usually be represented by the 


symbol 1. 


(a) The component of force F, acting in the glide direction of the dislocation 
of fig. 19. The unit of force is bD/y, where y is the (constant) distance of 
the glide plane from the dislocation at the origin. (6) The component of 
force Fy in the new glide direction if the moving dislocation of fig. 19 
has its Burgers vector turned through a right angle while the dislocation 
at the origin is unaltered. The unit of force is bD/x. 
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The activation energy per unit length if the dislocation is to escape from 
its bound position is given by 
W ,=bD(4—4 In 2)=0-153 BD. (68) 


4.62. Force between Screw Dislocations 

Consider two cylindrical holes of radius 79, with axes a distance r apart. 
If a slit is made along a plane joining the axes, and the sides of the slit 
are displaced +3) parallel to the axes, the result is a pair of screw 
dislocations. Since stresses round screw dislocations, like those round 
edge dislocations, vary as 7-1, we may use Eshelby’s result that, in the 
limit r)<r, it is not necessary to annul the normal forces over the two 
cylinders, and the stress may be taken as a simple superposition of the 
stresses produced by each dislocation in a medium with no other 
singularities. With this approximation, the interaction energy is 


W=(b?p/2z7) Inr, ae er ete 
and the force is 
F620. 9 a Oe oe 


4.63. Force between Parallel Dislocations with Burgers Vectors not Parallel 

It is obvious from the form of the function giving the energy density 
at any point that there is no interaction between edge and screw 
dislocations lying along parallel lines. It remains to consider two dis- 
locations lying along lines parallel to Oz, and having Burgers vectors 
(b, 0, 0) and (0, 6, 0). We may use the method of Cottrell and Bilby, 
considering the first dislocation already existing at the origin, and the 
second formed by displacements over the surface R>r>ryp, 0=@). The 
energy of interaction is easily shown to be 


W=2)D sin 20, « sk igpls acount a eae 
and the force in the glide direction of the second dislocation is 
F,=—bD2x"! cos? 6 cos 20=bp,,_. «= eae ana 


When the integration is made over a radial cut, the angular dependence 
of W arises entirely from the stress function x1 of eqn. (52). The force F, 
is plotted as a function of y in fig. 20(b). There is stable equilibrium 
for =} and for 6=—}n, and unstable for 2=—4n and 6=}7. The 
. position of lowest energy is given by 6=—}7; in this the zone of 
compression of the dislocation (0, b, 0) overlaps the zone of dilatation of 
the dislocation (b, 0, 0). The activation energy per unit length to move 
the dislocation from y= 0 to y=—=2 is 


Wa=tbD.. OO er ees 


In general, if a dislocation lies along a line with direction cosines | we 
may resolve its Burgers vector b into an edge component e and a screw 
component s given by 

e=—I[a(lab) 
and s= Ib. 
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It two dislocations lie along parallel lines, the shortest path between which 
1s a vector r, their interaction energy per unit length is given by 


W= 5 [5182 nr+(1—o)-eye, nrteyr er)... (74) 


4.64. Forces between Dislocations along Perpendicular Lines 

In §§ 4.61-4.63 we have considered the general interaction between 
dislocations lying on parallel lines. If the dislocation lines are 
perpendicular, the force between them is not constant along the length 
of each, but is concentrated in the region where they approach closely. 
If the lines are initially straight, their interaction may cause a bulge to 
appear in each where it approaches the other. In § 5.22 this difficulty is 
avoided by considering the interaction between an infinite grid of screw 
dislocations and a single screw dislocation lying parallel to the plane of 
the grid and perpendicular to each dislocation of the grid, and far from 
the grid. In the present section we shall neglect this complication, and 
shall also neglect the difficulties which arise on an atomic scale if one 
dislocation cuts another, so that the regions of misfit in their cores are 
superimposed through a small volume. The total interaction energy 
from this last cause is cb°~, where c is a numerical coefficient probably 
of the order of 0-1 or 0-2, and the force is small unless the dislocations 
approach within say 36 of each other. 

We may take the first dislocation to lie along the z axis, while the - 
second lies parallel to the y axis and passes through the point (X, 0, 0). 
Then the essentially different configurations which are possible are those 
listed in the following table : 


Burgers vector of Burgers vector of 
Serial dislocation along dislocation along 
[001] [010] 


[100] (edge [100] (edge) 


1 ) 

2 [100] (edge) [010] (screw) 
3 [100] (edge) [001] (edge) 
4 [010] (edge) [010] (screw) 
5 [010] (edge) [001] (edge) 
6 [001] (screw) [010] (screw) 


4.641. Two Edges, both Burgers Vectors along Shortest Distance between 
Lines.—Although the energy density contains cross terms in p,, and 
in the dilatation, there is no interaction between the dislocations, as 
may be seen by forming the second dislocation by a displacement over 
the cut z=0, w>X. This configuration occurs in pencil glide (Taylor 
and Elam 1926) and cross slip (Maddin et al. 1948). 


4.642. One Edge with Burgers Vector along Shortest Distance, One 
Screw.—No interaction, by the argument of § 4.641. If the edge passes 
through the screw it develops a jog (see footnote to Mott 1951), so that 
it is still a pure edge, but no longer lies in one glide plane. 
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4.643. One Edge with Burgers Vector along Shortest Distance, One 
Edge with Burgers Vector Parallel to Line of First—No force on second 
in its own slip direction, because the field of the first does not vary in 
this direction. Considering the first formed in the presence of the second 
by a displacement over the cut y=0, «<0, we find the interaction vanishes. 


4.644. One Edge with Burgers Vector Parallel to Line of Second, One 
Screw.—No interaction, by the argument of § 4.641. The screw can pass 
through the edge, which then has a kink, but still lies entirely in one 
glide plane. 


4.645. Two Edges, each with Burgers Vector Parallel to Line of other.— 
No force in glide direction of either, by the first argument of § 4.643. 
There is a finite interaction leading to non-conservative motion. We 
estimate this by considering the first dislocation to be formed in the 
cylinder 2?+y°<R?, 2<k?R?, R>X, k>1. The stresses are given in 
cylindrical coordinates by the stress functions (49) and (52), and the 
cylindrical surface is free from stress. There is an axial stress given by 
Pez=N Por +Poo)/2(A+p), and a state of plane strain can only be 
maintained if normal stresses act over the ends of the cylinder z= +kR. 
Now let. the second dislocation be made over the cut z=0, a>X. The 
energy of interaction is proportional to the integral of p., over this cut, 
which does not vanish, but contains a term linear in X leading to a 
constant force between the dislocations. During the formation of the 
second dislocation, the ends z=-+-kR of the cylinder move. Since k>1, 
they cannot distort but must move rigidly. The work done by the 
stresses acting on them therefore depends only on the total force and 
total couple on each end, each of which can easily be shown to vanish, 
so that the interaction depends only on the integral of p,, over the cut. 


4.646. Crossed Screws.—Finite interaction, which is discussed in 
§ 5.22. Two screw dislocations cannot glide through one another, for 
after they have crossed each contains a jog which is an element of edge 
dislocation with its Burgers vector perpendicular to its direction of 
motion. Any further motion must be non-conservative. 


4.647. Summary of the Laws of Force between Dislocations—The 
results of §§ 4.61, 4.62, 4.63, 4.64 and 5.22 may be summarized in the 
following rules, which apply to dislocations on parallel or on perpendicular 
lines: two screw dislocations always cause each other to glide along the 
shortest distance between them. Two edge dislocations cause each other 
to glide if, and only if, their lines are parallel. If two perpendicular 
dislocations each have their Burgers vector parallel to the line of the 
other, there is an interaction which does not lead to glide. In all other 
cases there is no long range interaction. 


4.65. Dislocations in Thin Plates 
Kshelby and Stroh (1951) have considered dislocations in thin plates. 


The force between two screw dislocations lying normal to the plate falls 
off exponentially with the distance between them. 


Mathematical Theory of Stationary Dislocations 311 


4.66. Force between a Dislocation and a Solute Atom 

The interaction between dislocations and foreign atoms dissolved in 
the lattice forms the basis of the modern theory of the yield point which 
is observed in mild steel and other materials (Sutton and Willstrop 1937, 
Cottrell 1948, 1952 a, Nabarro 1948 a, Cottrell and Bilby 1949, Cottrell and 
Jaswon 1949, Bilby 1950b). If the dissolved atom is regarded simply 
as a centre of dilatation occupying a volume 8V greater than that of a 
normal atom (in the case of a substitutional solute) or causing an 
expansion 6V (in the case of an interstitial solute), it does not interact 
with the shear stresses in the matrix, but only with the hydrostatic 
pressure P. A screw dislocation produces no hydrostatic pressure in 
regions where linear elastic theory is valid; the pressure produced by 
an edge dislocation is given by (50) as 


P>=—Pyy— Poo Pez 

3A+2u bur sin 6 

~ A+2u° T ; 

The energy of interaction is psV. If one looks into the physical inter- 

pretation of this mathematical centre of dilatation (see references to 

§ 4.3), it becomes clear that the problem needs to be defined in more 

detail. The dissolved atom must be regarded as a sphere having definite 

elastic properties. Changes in the assumptions alter the interaction 
energy by a numerical factor. 

Frank, in discussion of Cottrell (1948), pointed out that in some cases, 
such as that of an interstitial solute atom in a body-centred cubic lattice, 
the solute also interacts with shear stresses in the matrix. In these 
cases there is also an interaction between solute atoms and screw 
dislocations, and these are the systems which show the most clearly 
marked yield phenomena. Crussard (1950) has given a quantitative 
theory of this interaction. 


(75) 


§5. ARRAYS OF DISLOCATIONS 


We shall now discuss the properties of various arrays of dislocations. 
Most of the arrays contain dislocations uniformly spaced, but in § 5.1 
we consider a more general case in which dislocations moving on the same 
glide plane reach equilibrium in a pattern which is not uniformly spaced. 


5.1. A ROW OF DISLOCATIONS ON THE SAME GLIDE PLANE 


Consider a set of m similar edge dislocations in the same glide plane, 
each subject to a force f tending to move it along the plane. Suppose 
the dislocation lying furthest in the f direction is prevented by some 
obstruction from moving across the glide plane. Cottrell (1949) has 
pointed out that the force per unit length acting on the obstruction is nf. 
For if the obstruction moves a distance d in the direction of f, each of 
the n dislocations moves the same distance, and the work done by the 
external forces is nfd. Thus a glide plane containing mobile dislocations 
produces stress concentrations analogous to those produced by a crack. 
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Taylor (cf. Starr 1928) and Zener (1948 a) have suggested that such a 
glide plane may have no resistance to shear, so that the stress 
concentration at its end is actually equal to the stress concentration 
produced by a crack. For a crack of length L, with a radius of 
curvature p at the end, this stress concentration has a maximum value 
of (2L/27p)1, at a point a short distance beyond the tip of the crack 
(Starr). 

Consider a glide plane of width L, and suppose that slip begins at one 
end and cannot continue right across because the glide planes of the 
neighbouring crystal are inclined to those of the first crystal. The 
repulsion between similar dislocations means that under a given applied 
stress P only a definite number n(P) of dislocations will enter the glide 
plane. The force acting on unit length of the leading dislocation is 
then bPn(P). Since n(P) is clearly an increasing function of P, the stress 
on the leading dislocation increases more rapidly than P, whereas the 
concentrated stress produced by a crack is proportional to the applied 
stress. This problem has been treated by Eshelby, Frank and Nabarro 
(1951), who find that, to a close approximation, 


n(PJ=PL2D, 3 | eee 


where D is given by (4). This is the stress concentration factor K ,. 
We may compare it with the factor for a crack K,=(2L/27p)!” by putting 


p=b. We find 
Ke 27 Le P 
x- (=) m(1—a) mere tire | fe 


which is equal to unity when 


ye 1 2b \1/2 
z= maya) . eK a a me (78) 


The concentrated stress is then given by 
P 2 
mm oF 27r(1—oa) 
If the concentrated stress is less than this, K,<K,, while for higher 
stresses K y>K,. 

For some purposes the important quantity is not the stress on the 
leading dislocation, but the stress distribution in a fairly wide region 
ahead of the tip of the glide plane. It can be shown that, in a region 
ahead of the leading dislocation which is a distance more than D?/4LP? 
and less than L from the leading dislocation, the stress is nearly equal 
to that which would be produced by a crack of length 4L. Within this 
distance D*/4LP* of the leading dislocation, the stress of the leading 
dislocation itself predominates, while at large distances the stress 
exceeds P by the stress field of a group of n dislocations at the origin. 
The relative displacement of the two sides of the glide plane near its 
tip oscillates around the value of the relative displacement of the sides 
a crack of length 4/ under the same stress. Koehler (1952) has estimated 
the tensile stresses near the leading dislocation. 


“0-035; |.) > ance 
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The method used to obtain (76) depends on a result of Stieltjes (1885) 
its application to this problem was suggested by H. Heilbronn. 
Suppose the leading dislocation is held up at the point x,=0, and the 
remaining dislocations are in equilibrium at points %o....%, such that 
@y<%2<....<%,. Then the ith dislocation is acted on by a force 
f=—Pb tending to reduce its coordinate, and by a force of interaction 
bD|(x;—x,;) with each other dislocation j. The equilibrium values of 2, 
are determined by the set of equations 
a De StU Le Meee iEO fue 1s Es Oyun, & on 80) 


a 
jJ#U 


and we require the largest value of n for which x, <L. 
Put 


(x)= (e—a@,)(e@—ag)....(@—2,), . « » 281) 


so that the shear stress produced by all the dislocation at any point « 
on the glide plane is DF’(x)/F(x). Put also 


F(x) Pe Dest ae - 
9 (x)= Sia) A oe Ae ae 7; (82) 
pHi 
Then eqns. (80) may be written 
g(x ;)=a, 
where a=P/D=27(1—c)P/bp. ae Me sle en (O33) 


As x tends to 2;, so g,(z) tends to $F" (x,)/F’(x,;). The problem is thus 
reduced to the determination of a polynomial of degree n which vanishes 
at «=0 and at n—1 distinct positive values v,....2,, and for which 
4F"(x;)/F'(x;)=a at each zero. The polynomial will certainly satisfy this 
condition if it satisfies the differential equation 


Ole 1 Nae Oi ot a (84) 


and it is easily shown that this differential equation has in fact a 
polynomial solution of degree n satisfying the required conditions. For 
small x, the second term in (84) may be neglected, and the equation 
solved to show that the first v dislocations, where v is small compared 
with n, are contained in a distance proportional to v?/na. The total 
number n of dislocations occupies a width w, of the glide plane, and it 
will be shown that for large n the value of x, is given very closely by 


L,=2n/a. eer soi) | iy. (80) 


Equating this to L, we find that the maximum number of dislocations 
which can be held in the glide plane is given by (76). 

To justify (85), we proceed as follows. First suppose that all the 
dislocations except the last are constrained to lie at the origin. The 
force they produce on the last dislocation is (n—1)b*y/X,,; where Ae eis 
the equilibrium position of the last dislocation in this constrained system. 
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Equating this to Pb, we have X,=(n—1)/a. But it is clear that if all 
dislocations but the first are now released, all will move away from the 
first, and so v, >X,—=(n—l)/a. Heilbronn has obtained the closer result 
that a, >(2n—yn")/a, where y is an absolute constant. The determina- 
tion of a lower limit to 2, is based on the differential equation (84). 
We write 
F(2)=ue™; 

and (84) becomes 

u" + (2nax-1—a?)u=0. oh. Nha ted ee 


Suppose 2, >2n/a, and choose the sign of the leading term in the 
polynomial solution of (84) so that F and wu are positive when x>z,,. 
Then w’(x,)>0, and u’(x)>0 for 2>2,. Integrating, u’(x)>u’(x,), and, 
integrating again, u(7)>(a—z,)u'(a,) for <>a,. But in the polynomial 
solution of (84), F(a)=O(x") as x+oo, and so u(x)-+0 as x>oo. This is 
a contradiction, and the assumption 2, >2n/a is false. Hence 
2n—yn'¥8 <ax,<2n. 

A further case of interest is that in which no external stress is applied, 
and the first and last dislocations are held by obstructions. This problem 
and a number of cthers were discussed by Eshelby et al. Nye (1949), 
in experiments on crystals of silver chloride, met a related problem, in 
which the dislocations were held up not by internal obstacles in an 
effectively infinite medium, but by the surface of the crystal. If the 
glide plane y=0 is bounded by the edges of the crystal z=-+JL, the 
problem must be solved with the additional condition that the stresses 
Pre and p,, are zero on the surfaces x= +L. The solution was obtained 
by Mann on the assumption that dislocations are distributed at equal 
spacings across the glide plane. In all of these problems the glide plane 
becomes bent in the way shown in fig. 2. It does not seem possible to 
extend the method of Stieltjes to the case of a slip plane containing both 
positive and negative dislocations. Following Leibfried (1951), we may 
consider the limiting form of the distribution when the total number of 
dislocations is very large. We suppose that the number of dislocations 
lying between x’ and 2’+dz’ on the glide plane is D(x’) dx’. At the 
point x these produce a stress 

| DAC T OS i 
so «2 
where the Cauchy principal value of the integral is understood. The 
condition of equilibrium of a dislocation at the point 2 is thus 
eh) ie) ae ie 
| te eR aye 


where P is the external stress. This equation must hold at all points for 
which Y(x) is not zero. The derivation of (88) is purely formal, for the 
stress (87) has a singularity at each dislocation in the distribution, and 
these are the only points for which (88) should be true. However, 


(88) 
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comparison of (88) with § 2.1 of Part III shows that (88) is the condition 
that the glide plane containing dislocations should have no resistance 
to shear, and the analysis of Eshelby et al. shows that, in a particular 
case, this is correct. 

Using (88), Leibfried has solved a number of problems, including all 
those previously treated by Eshelby et al. He has also treated the case 
of a slip plane of length 2a containing equal numbers of positive and 
negative dislocations, forming a model of a slip plane containing a 
Frank—-Read source. The solution of (88) for this problem is 


Ot) — Rela eee a sw. (89) 


The total number of dislocations of either sign is 
| eee ris pee ee G0) 
and the shear strain is proportional to 
bf LD) dxu=bPa?/2D=r(1—o)a* Pip, . . « (91) 


which is proportional to the applied stress. Segedin’s theory (1951) of 
a penny-shaped crack in a sheared elastic body may be interpreted 
similarly. 

Leibfried also considered a glide plane containing less than Pa/7D 
dislocations of each sign, representing a Frank—Read source which has 
become blocked before producing the maximum number of loops which 
the glide plane can accommodate under the applied stress. 


5.2. THE GRAIN BOUNDARY AS A WALL OF DISLOCATIONS 


5.21. A Wall of Edge Dislocations 

The first calculation of the stresses produced by joining two crystals 
whose axes are inclined to one another at a small angle is due to Taylor 
(1934). He considered two simple cubic lattices with lattice spacing 6, 
inclined to one another at a small angle %, so that the z axes in the two 
crystals are parallel, while the surface of separation is parallel to this 
axis and bisects the angle between the crystals. We call this a tit 
boundary. Take the y axis in the boundary surface. Then before the 
crystals are joined together they have the form shown in fig. 21, where 
a step occurs after every 2/5 atoms. The join is made with least energy 
if the step M on the left is brought near to the gap N on the right, when 
the step P on the right falls in the gap Q on the left. If the opposing 
surfaces are now pressed together, each suffers a saw-tooth deformation 
of amplitude+4), and of wavelength 6/¢. 

By a Fourier analysis of this deformation it is easily shown that the 
stress P at a point distant « from the surface of misfit is of order 


P=2uiJm e-”, where n= 2m} | x | /b. taageet ee 192) 


P.M. SUPPL.—JULY 1952 
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The stress is only appreciable at distances from the surface of misfit 
less than about 2b//, which is the distance between successive steps 
in the face of a crystal. Taylor suggested that these regions of high 
stress were the boundaries of the mosaic blocks which were postulated 
by Darwin to account for the intensity of x-ray reflections, and later 
identified with various microscopically observed structures in crystals. 


Fig. 21 


C ef aI > > OT, © _ € . oo 5 a 7 od “4 - 
[wo grains inclined at a small angle #. They join with least energy if the 
corners P and M fit in the gaps Q and N. 


Dislocations near the middle of a mosaic block ean move under an 
external stress equal to that given by (92), where the value of 7 is that 
for the middle of a block, and this should be the elastic limit Taking 
the width of a block as L=2x=—10-4 em, b=3x 10-8 em b=3x 10-4 
radians, and p= 2-6 x 101! dyne em~*, the calculated elastic limit becomes 
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260 g. wt. mm™, which is a reasonable order of magnitude. Taylor 
assumed that dislocations could not cross the mosaic boundaries at low 
temperatures, but could penetrate regions of higher and higher stress in 
the boundaries as the temperature was raised. 

Burgers (1939 a) pointed out that the composite crystal could be 
regarded as a single crystal containing edge dislocations of vector (0b, 0), 
spaced at regular intervals b/s along the y axis. The value of 7 for the 
middle of a mosaic block is 7 times the number of dislocations in a wall 
of the block. This point of view makes obvious a result which is implicit 
in Taylor’s analysis, that the stress is concentrated along a row of parallel 
lines in the xz plane, and falls off rapidly between these lines. Bragg 
(1940) has also considered boundaries of this type. 

It is easily verified that if a crystal contains the infinite row of 
dislocations specified in the previous paragraph, the regions far from 
the x plane and lying on either side of it are unstrained, and rotated 
through angles +4 about the z axis. At these distant points the 
uniform row of dislocations may be regarded as a continuous distribution 
of dislocation spread over the x plane with density %/b. The stresses 
and displacements ® may be obtained from those ¢ given in § 2.1 by 
integration, using the formula 


L 
d= lim $[ bay, ree ek (93) 


I+ 
The stresses p,, and p,,, are odd in y, and their integrals vanish. ‘The 
integral of p,,, converges to zero as Loo, and the crystal is thus free 
from stress. The rotation produced by a single dislocation is given by 
| Omae® 
on Py? 
Substituting in (93), the rotation produced by a sheet of dislocations is 
found to be —4% sgn w, where sgn x is 1 if w>0 and —1 if x<0. 

It is easy to make a rough estimate of the energy of such a wall of 
dislocations (Shockley and Read 1949), by using the formulae of § 4.41 
in conjunction with the result that the stress is confined to a strip of 
width roughly equal to the distance between dislocations. We regard 
the total energy as the sum of equal contributions from each dislocation, 
this contribution arising from a circular cylinder of radius 6/2 with the 
dislocation lying along its axis. Within this cylinder the stress is 
approximately that produced by a single dislocation in an infinite medium. 
Since p,,=,, the energy density is $u~"p,0° +3(A+H) Pre’ If this is 
integrated from b/2% right down to r=0, the result is infinite, and we must 
cut the integration off at some radius 79, of the order of 6. We expect 
that r, will be independent of %, provided that is reasonably small, 
that is to say, that the energy in the core of a dislocation is little affected. 
by the presence of other dislocations more than a few times b away 
from it. A discussion of this assumption and an estimate of 77 depend 
on assumptions about the atomic forces in the core (§ 2.4 of Part III). 


NeZ 


CHO =} 
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The energy associated with unit length of each dislocation is then 
bu b 
47(1—o) ne 2rob ; 


which we may write in the form 


bu bo 
r0 as 4 
tml aie meee. Bee (94) 
with 
p=] rs. 


The distance between dislocations in the wall is b/s, and the energy 
of unit area of the wall is thus 


by if 
imino)? Oy: | OLY Se) en 


A more detailed discussion of the energy of a grain boundary is given 
in § 5.23. 

The formation of a wall such as that shown in fig. 22 in a region 
containing a number of dislocations of the same kind does not require an 
activation energy. Any two dislocations arranged as in fig. 19 with the 
angle 6 lying between 7/4 and 37/4 will move until the line joining them 
is perpendicular to their glide plane. Suppose that the process of building 
up a wall has continued until the wall of fig. 22 has been formed, which 
contains dislocations at a uniform spacing of w stretching from (0, L) 
to (0, —L). This is the process of polygonization observed by Lacombe 
and Beaujard (1947), Lacombe (1948), Cahn (1948, 1949, 1950), Guinier 
and Tennevin (1948, 1950) and Dunn and Daniels (1951). The theory 
was first discussed by Cottrell (1949). The force in the glide direction 
on a dislocation at a point (x, y) at a distance from the wall large compared 
with w is given by eqn. (67) with sign reversed as 


sin 40 dy 
Fea) es 
where tan 6=(y—n)/x. 


L 
F,=16D | 


This reduces to 

2LbD x(e?—y*?+ L?) 
The denominator of this expression is always positive (except at the ends 
of the wall), and the force tends to pull the dislocation into the plane 
of the wall or to drive it out according as (#2—y?+L*) is positive or 
negative. The region of attraction is bounded by the hyperbola sketched 
in the figure. 

When the dislocation comes within a distance w of the wall, it is no 
longer justifiable to represent the sum over all the dislocations in the 
wall by an integral. If its glide plane coincides with the glide plane of a 
dislocation already in the wall, the position of equilibrium will clearly 
be one in which neither dislocation lies exactly in the plane x=0, while 
a dislocation gliding midway between two dislocations in the wall may 


f= (96) 
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reach this plane. We shall consider in detail the force on a dislocation 
of an infinite wall if it glides out of the wall. The force when it is at the 
point (x, 0) is 
© sin 40 
F,=4bD ZX 
n=1 Ww 
where tan 6=nw/x. 


Fig. 22 


> 


<&> 


1 ization. A wall of dislocations stretching from y=—L to y=L exerts 
ie bie range attraction sweeping into the line of the wall dislocations 
above and below the two branches of the hyperbola. The dislocation 
shown, which lies between the branches, is repelled. 
This may be written 
y 26D 2 1—n*w? |x? 
ion HL guy (1? /227)? 


2 pd ; 
= 2 [1-S [sine SZ]. ee ee (97) 
x 
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For x small compared to w the expression in square brackets tends to 
—7*z?/3w?, and the restoring force is —7*bDa/3w*, while for large w the 
expression tends exponentially to —1, and the force tends to —bD/x. 
The force is plotted in fig. 23 (a). 

The force on a dislocation moving in the presence of an infinite regular 
wall along the glide plane of one of the dislocations in the wall is obtained 
by adding the repulsion F,,=bD/x of a single dislocation at the origin. 
It tends exponentially to zero for large x in accordance with Taylor’s 
formula (92) for the stress. The last formula in Burgers’s (1939 a) paper 
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Force on a dislocation approaching a wall of dislocations spaced at equal 
distances w. The unit of force is bD/w. (a) There is one dislocation 
missing from the wall, and the mobile dislocation is moving in this glide 
plane. (The sign of /, should be reversed.) (6) Dislocation approaching 
a complete wall along the glide plane of a dislocation in the wall. 
(c) Dislocation striking the wall midway between two dislocations. 


gives the stress p,, produced by a complete wall at any point, and 
the force on a dislocation is bp,,. For a dislocation moving on the glide 
plane of a dislocation in the wall the force is 


x 
ae mee ie 


in agreement with eqn. (97), while in the case of a dislocation which 
would strike the wall midway between two dislocations the force is 


es, bD w?x? mu 
f.=- — [ cosh ™. arch! eee ae eel 


x ww 
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In the second case the dislocation is attracted into the wall. Its binding 
energy when it reaches the wall is 


WDD Ato ueeee me. >. (100) 


The forces are shown in curves (b) and (c) of fig. 23. Seeger (1951) has 
evaluated the forces on dislocations lying on intermediate glide planes. 


5.22. A Wall of Screw Dislocations 


There is a radial repulsion between two similar edge dislocations, but, 
because the interaction energy of a pair of edge dislocations depends on 
the direction of the line joining them, and because they can move only in 
parallel glide planes, it is possible for their interaction to lead to a 
diminution of the distance between them. Neither of these factors applies 
to the interaction between two screw dislocations, and they will always 
repel one another. Nevertheless, it is possible to form walls of screw 
dislocations which are stable provided that the dislocations at -the 
boundaries of the wall are held fixed. These walls consist of grids of 
screw dislocations lying along two perpendicular directions (Frank 1948, 
van der Merwe 1950), and may be called twist boundaries. 

We begin by calculating the effect at distant points of an infinite row 
of screw dislocations having Burgers vectors (0, 0, 6), spaced at equal 
distances w along the y axis. From eqns. (14) and (93) we see that the 
only stress which does not vanish is 


bu (? «dy by. 


Pyz= Fay | ey see choy beget key #) 


The rotation is 


w= sgn 2. ite te (LOZ) 
Now consider a screw dislocation (0, 6, 0), passing through the point 
(X, 0, 0). The stress (84) exerts on this a force per unit length 
F,=—6*u/2w sgn X, which pulls it towards the wall. If, instead of a 
single screw dislocation (0, 6, 0), we have a sheet of similar dislocations 
spaced. w apart in the plane «=const., the force between the two sheets 
is an attraction of 62./2w? per unit area. If the two sheets coincide, a 
single dislocation of the second sheet which breaks away and reaches 
the point (X, 0, 0) is subject to a restoring force 


L ow 4 x2 Poy 
bu TA aX \—1 
=o Ee cotk (=) i pedo. Aelia: 


The function in square brackets decreases from 1 when X=0 to 0 when 
A = 00. 
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If the whole second sheet is in the plane =X, it produces a stress 
Py — gH sen (e—X) 
and a rotation 
w= ~ sgn (~a—X). 


When X=0, the stress cancels that given by (101), so that a wall of 
crossed screws produces no stress at large distances, while the total 
rotation is w,—b/2w sgn zx. Rotations of this kind have been observed 
in plastically deformed crystals by Heidenreich and Shockley (1948) 
and Wilman (1950, 1951). van der Merwe has illustrated the formation 
of such a wall of crossed screws by considering a crystal split along a 
plane in which the atoms lie in a simple square array. One half of the 
crystal is rotated with respect to the other, and the cut surfaces are again 
placed in contact. The resulting configuration is shown in fig. 24 (a). 
The forces acting across the boundary tend to pull neighbouring rows of 
atoms on the two sides into alignment, which leads to the configuration 
of fig. 24 (b), in which the grid of screw dislocations is indicated. 


5.23. The General Grain Boundary as a Wall of Dislocations 

As an example of a boundary containing edge dislocations of two 
kinds, we may consider fig. 25. It follows from the discussion of § 4.63 
that this array of dislocations is in stable equilibrium, and it produces a 
stress field tending exponentially to zero at distances greater than w from 
the wall. The lattices on either side of the wall are rotated through 
angles of +b/2,\/2w, and this mixed wall thus represents a boundary 
lying unsymmetrically between two grains. 

A more general case is shown in fig. 26. We use the notation of Read 
and Shockley (1950), but write ys where they write 6. The line OP, inclined 
at an angle ¢ to the x axis, forms the boundary between portions of the 
lattice which are rotated through small angles of +4. Unit length of - 
OP intercepts 6-1 sin (6—3}) w-planes and b~! cos (6—2#) y-planes of 
the lattice on the left, and b~! sin (6+ 43) x-planes and b~! cos (6+ 4%) 
y-planes of the lattice on the right. The total number of planes intercepted 
is, for small ys, b~' cos ¢ w-planes and bly sin ¢ y-planes, corresponding 
to b-' cos ¢ dislocations of type (b, 0, 0) and b-ks sin ¢ dislocations of 
type (0, 6, 0) in unit length of the wall. 

Shockley and Read (1949) assumed that each etch pit in the sub- 
boundaries observed by Lacombe (1948) marks the emergence of a 
single dislocation. They compared the spacing of etch pits with that 
calculated from the observed difference of orientation between neighbour- 
ing sub-grains, and found satisfactory agreement. (See also a footnote 
to Read and Shockley 1950.) 

Frank (1950 b) has shown how the boundary between two grains may 
be analysed into dislocations when the angle % between the grains is not 
small, and the axis of their relative rotation meets the boundary 
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obliquely. He introduced the idea of the median lattice @ from which 
the lattices 3B and € of the two grains may be derived by equal and 
opposite rotations, and considered ® and € to be derived from @ by 
inserting dislocations characteristic of the lattice A along the grain 


Fig. 25 
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A stable wall lying at 45° to each glide direction, with two types of dislocation 
equally spaced. 


Fig. 26 


A general tilt boundary. The two er 
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boundary. Since the rotation need not be small, the characteristic 
dislocations of @, ® and € are different. He defined corresponding 
points B and C in the two lattices to be points which are derived by 
rotation from the same point A of the original lattice. Take the origin 
of coordinates O in the boundary, and consider a point A which in 
the original lattice had coordinates r. Let the lattices B and € be derived 
from @ by rotations about a line | through the origin through angles 
+3%. The coordinates of the point B referred to the lattice @ are 
rg=r+lar sin $4-+1 a (lar)(1—cos 44), 
while those of C referred to the same lattice are | 
ro=r—lar sin $4+1,a (lar)(1—cos $y). 

The separation of the corresponding points B and C is thus 
BC=2ralsin3%. This is a dislocation composed of the last three of 
Volterra’s orders of dislocation (§ 3.2), but it is convenient to regard it 
as a continuous uniform distribution of dislocations of the first three 
orders. The points which, before the dislocation was made, coincided 
at r, are now separated by 2ral sin 4%, while neighbouring points at 
rtdr are now separated by 2(r+dr)alsind%. The dislocation 
introduced across the grain boundary between the points r and 
r-dr is 2dral sin $4. We may introduce the vector s=2I sin $4, which 
defines the axis and amount of the relative rotation of the two lattices, 
and is independent of the position of their boundary. Then the density 
of dislocations in any-part of their boundary is defined by the statement 
that a line r in the boundary crosses dislocations of total strength ras. 

Read and Shockley have given a detailed analysis of the energy of a 
grain boundary of the type shown in fig. 26, and Read (1952) has 
considered a general grain boundary. The calculated dependence of the 
energy on % and ¢ agrees well with the measurements of Dunn and 
Lionetti (1949) and Dunn, Daniels and Bolton (1950). A review of the 
experimental results has been given by Fisher and Dunn (1952). The 
energy is most conveniently calculated by integrating the tangential 
stress across the slip plane of each dislocation from infinity to the actual 
position of the dislocation (§ 4.6). It is given by 


eee) ee eo,» (104) 

where 
By= ps (ein $-+008 4) Rie ate 7s) .(105) 

and 


(sin d+_cos 4)A=(sin ¢-+cos ¢)Ay—sin ¢ cos ¢ (sin d-+-cos ¢) 
—(sin d Insin¢d-+cos¢Incos¢). . . . (106) 
The constant A, corresponds to the constant In 4 of (94). It is a measure 
of the energy in the core of a dislocation, and occurs with a factor 
7% (sin d+ cos 4), which measures the total number of dislocations in 
unit length of the wall. The value of A—A, given by (106) varies 
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between —0-15 and +0-13. Sinee the estimated values of A, lie between 
0-5 and 0-8, the variation of A with ¢ may usually be neglected. (The 
validity of this argument is now in doubt. See § 2.42 of Part III.) 
If we neglect the variation of A with ¢, (104) and (105) show that the two 
boundaries B,LB, and B,L’B, of fig. 27 have the same energy. Read 
and Shockley also show that a movement of all the dislocations along 
their glide planes from one boundary to the other does not produce a 
relative displacement of the two sides of the crystal, provided the 
angle js is small. A relative displacement is produced only if the 
dislocations move non-conservatively by diffusion. 


A small angle tilt boundary (after Read and Shockley). The y dislocations are 
not shown. The wall can move by glide alone from B,LB, to B,L’B,. 


If the angle % is not small, other considerations arise. Sometimes 
more than one pair of angles (w, 4) will provide a reasonable description 
of the same boundary. The analysis of particular cases shows that the 
predicted properties of such a boundary are independent of the 
description employed. Moreover, we have tacitly assumed that 
dislocations are spaced evenly along the boundary. If, for a given ¢ and 
a given small %, the necessary mean spacing is 15} lattice planes, we shall 
actually have groups of dislocations with spacings 15, 15, 15, 16, forming 
a ‘superlattice’. This small deviation will not alter the energy much, 
but if % is larger, so that the spacing should be, say, 34 planes, the actual 
spacings will be alternately 3 and 4, The energies given by (104)-(106) 
are calculated for those values of y% and ¢ for which both x and y dis- 
locations have integral spacings. Read and Shockley show that small 
deviations from this orientation cause the energy to rise steeply. The 
plot of energy against orientation shows a cusp similar to the # In x cusp 
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at the origin, but of smaller size. Since the integral spacings of 
dislocations of each type are destroyed by changes in ¢ as well as by 
changes in 4, boundaries of this kind, in particular twin boundaries, have 
strongly preferred orientations in the two lattices. 


5.3. LATTICES OF DISLOCATIONS 


The idea of a lattice of dislocations was introduced by Taylor (1934), 
who considered that the distribution of positive and negative edge 
dislocations in a strained single crystal could be approximately 
represented by a regular lattice. Since the stresses used by Taylor were 
incorrect, the lattices he found to be stable differed somewhat from those 
we shall consider. The present treatment is in fact mathematically 
simpler, and the results are similar. However, our present ideas of the 
mechanism of work hardening suggest that arrays of dislocations which 
could reasonably be represented by such a lattice occur only in kink 
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A lattice of edge dislocations. (The top dislocation at the left would more 
usually be represented by the symbol 1.) 


bands, and in these bands the dislocations are probably no longer free to 
move in their original glide planes. A model of more interest is that of a 
Frank—Read source in a finite crystal embedded in a matrix into which 
the dislocations cannot penetrate (§ 5.1). 
Figure 28 represents a lattice of positive and negative dislocations in 
which the positive dislocations lie on glide planes separated by d in 
columns at distances a. The negative dislocations lie in glide planes 
midway between these, in similar columns displaced a horizontal 
distance x. If the positive dislocations are held fixed, there will clearly 
be equilibrium when x=0 and when x=3a. it follows from eqns. (98) 
and (99) that there can be no equilibrium if all the dislocations are 
released, for the forces between any two columns are repulsive. 
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Possibly other arrays, in which positive and negative dislocations form 
close pairs, may be stable. When the positive dislocations in fig. 28 are 
fixed, the force on a negative dislocation is given by eqn. (99) as 


bD 
Fy 7 ba 
(ana) (107) 
be x+na 7 
where f= 2 , a / cosh? ei 


This force is plotted as a function of x/a for various values of d/a in 
fig. 29. The position z=0 is always one of unstable equilibrium, and 


Fig. 29 
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The force f;, between the positive and negative dislocations of fig. 28, for various 
values of d/a. 


a—=}a one of stable equilibrium. The maximum force fi, and the slope 
of/ex at x= a are given in the following table. 


dja 0-1 0-5 1 2 
ae 0-45 0-45 0-38 0-097 
—a Of/dx 0-00 0-49 1:89 0-47 
(a/d)f? max 2-0 0-80 0-14 0-003 


If a shear stress P is applied to the crystal, each negative dislocation 
is acted on by a force —Pb, and it comes to rest at a point « determined 
by the equation 


FS Philo 16 ee ee ay 
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If Pd is greater than F,,,,, this equation has no root, and the lattice of 
negative dislocations moves through the crystal. 


5.31. Taylor's Theory of Work Hardening 
Taylor assumed that the elastic limit of a work-hardened single 
crystal was given by 
Pb=f 


max 


or eee buf aeiod (=o) eeu oil, ey ho ed 09) 


He also assumed that at low temperatures the lattice of dislocations would 
be formed by dislocations moving into the crystal from mosaic walls 
separated by a distance LZ. The plastic strain during this movement is 


e=Lb/ad. Se Metre shen ee LLO) 


At higher temperatures the mosaic walls would be partly permeable to 
dislocations, so that some of the dislocations would have travelled 
distances greater than L before forming the array of fig. 28, and the plastic 
strain would be correspondingly greater. 

Eliminating ad from eqns. (109) and (110), we obtain 

Tere, OPS fe 
e 4(1—o)*Ld’ ™* 

The parabolic relation between stress and strain is in good general 
agreement with experimental results on cubic materials. The number 
of atoms in the side of a mosaic block is at low temperatures equal to L/b, 
or p</P? multiplied by the numerical factor af*,,,/4(1—o)*d, which 
for moderate values of d/a lies between 0-1 and 1. For ie at 
—185°c, the results of Boas and Schmid (1931) correspond to 
p«/P?=1-56X 10°, giving L/b-+6x 104, which is of the same order as 
the observed spacing of various imperfections in metals. The same 
agreement in order of magnitude is found for Al, Ca, Au, Fe and NaCl, 
at room temperature. ; 

This agreement in form and. order ai magnitude is very satisfactory. 
However, the calculation considers the metal as having only one glide 
plane, whereas face-centred cubic metals have four glide planes. Single 
crystals of hexagonal metals have only one glide plane, and they show 
an increase of elastic limit with strain which is very much smaller than 
that calculated, and is approximately linear in the strain. 

5.32. Influence of Dislocations on the Shear Modulus 

Masing (1944, 1948) has pointed out that the presence of a lattice of 
dislocations reduces the shear modulus of the crystal. The displacement 
of the lattice of negative dislocations corresponding to eqn. (108) produces 
a reversible shear strain! If P is small, the. re of the lattice 
is given by 6x Bilge Pb, and the strain is <’=bédaz/ad. This may be 
written 


(111) 


«’=2(1—o)P/(ua Of/ex), 
while the elastic strain in a perfect crystal would be 
e=Piu. 
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The fractional change in modulus is 

e’ /e=2(1—a)/(adf/dx). < Pe) 4 ents hoe 
This quantity is independent of the size of the mosaic blocks L, and of 
the density of dislocations in the lattice, which is proportional to l/ad. 
It depends only on the ratio d/a. Its value is plotted in fig. 30. 
Masing’s results for large d/a are similar, but for small d/a he finds a 
value of ¢’/e which decreases steadily as d/a decreases. This is because 
he considers a lattice which would be stable only if the stress around a 
dislocation were that given by Taylor. This lattice, corresponding to 
a=0 in fig. 28, remains stiff if a increases with d held constant, whereas 
our lattice, corresponding to =a, becomes soft as a increases. 
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Fractional change in modulus ¢’/e produced by the presence of the lattice of 
dislocations of fig. 28, for various values of d/a. 


Equation (103) gives the restoring force on a screw dislocation escaping 
from a crossed grid of screws. Adding the repulsion +6?u/27X of a 
dislocation occupying the missing place in the grid leads to a force 
between a crossed grid in the plane a=0 and a screw dislocation (0, b, 0) 
passing through the point (X, 0, 0) equal to 


2 
Fo=— ce E coth =| : 
2w w 


This is a repulsion diminishing exponentially with distance. On this 
basis a theory similar to that just given for a crystal containing sheets 
of edge dislocations may be developed for a crystal containing parallel 
crossed grids of screw dislocations, all right handed or all left handed. 
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Such a system forms an approximate model of-a crystal which has been 
twisted about an axis perpendicular to the glide plane. 

Since the maximum value of adf/éx in (112) is only about 2, the 
fractional change in modulus produced by a lattice of mobile dislocations 
is a large part of the modulus of a perfect crystal. In practice the moduli 
measured statically are changed only a few per cent by cold work. 
(They also agree closely with the moduli estimated from the speeds of 
ultrasonic waves, but this may not be significant because the frequencies 
of these waves are probably not high enough for the movement of 
dislocations to be suppressed.) We deduce that the system of dislocations 
in an ordinary metal is held in equilibrium by its interaction with 
immobile lattice defects and not merely by the forces between mobile 
dislocations. However, deviations of the order given by (112) have 
been observed by Crussard (1945) in zinc, by Crussard (1947) in 
magnesium, and by Andrade and Henderson (1951) in single crystals of 
gold and silver, using in all cases shears of the order of 0-001. 


$6. THe INTERACTION OF DISLOCATIONS WITH FREE SURFACES 


6.1, ELASTIC INTERACTION 


If an elastic dislocation lies near to a free surface, it is less constrained 
than a dislocation in the interior of a specimen, and its energy is lower. 
There is thus a force attracting a dislocation to a free surface. This 
force was calculated for an edge dislocation in a circular cylinder by 
Koehler (1941); we shall treat in detail the simpler case of a screw 
dislocation, and then give an outline of Koehler’s method. In both 
cases the force is approximately that due to an image dislocation of the 
same strength and opposite sign. 


6.11. The Screw Dislocation in a Circular Cylinder 

We use Eshelby’s method of § 4.42. The sum of the displacements (57) 
and (59) is constant on any circle passing through the points (€, 0) and 
(R?/é, 0). Therefore no stress acts across the orthogonal trajectories of 
these curves, which are the coaxal circles 


[(a—€P-+y?]—Al(e—B?/E)+-y"]=0, 


with these two points as limit points given by A=0 and A=oo. The 
circle x?+-y?—R? given by A=€/R? is free from stress, and may again 
be taken as the outer boundary of the dislocated body. For the inner 
boundary required to form a dislocation we take A=<*é*/(R4*— &*) where 
e is small. This is a circle of radius « around the point (€, 0), with its 
centre displaced from that point by a distance of order <*. By considering 
the work done on the cut surfaces if the dislocation is made over a cut 
from a=é+e to =f, the energy is easily shown to be 

bu, Re—&—E£e 
vs ln - Sete: ° 
P.M, SUPPL.—JULY 1952 Z 


We (113) 
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This is equal to (b%/47) In (R/e) when =0 and the dislocation is at the 
centre of the cylinder, and is zero when £=R—e and the inner and outer 
boundaries meet. The force on the dislocation is given by 
pe OW by 2&+<« 
o-— 5b ty PAE 
So long as the dislocation does not lie within a distance of order « from 
either the centre or the boundary of the cylinder, this may be written 


gine fe 


2 Or Ree 


where w is the distance between the dislocation and its image. 


F b21/2rrw, 


6.12. The Edge Dislocation in a Circular Cylinder 

Koehler evaluated the energy in this case by working in bipolar 
coordinates. A small cylinder surrounding the line of the dislocation is 
again excluded from the range of integration, although in this case the 
surface of the small cylinder is not free from stress. The force is 
approximately 

F=b*p/27(1—c)w=bD/w, 

where w is again the distance between the dislocation and its image. 


6.2. INFLUENCE OF THE SURFACE ENERGY 


6.21. Energy of an Exposed Ledge 

When an edge dislocation enters a crystal from a free surface, it 
exposes a ledge of width b. If the surface energy is y, the energy by per 
unit length must be provided, in addition to the energy of the dislocation 
itself. An equal surface energy appears when the dislocation leaves 
the crystal. If the energy of the dislocation exceeds by per unit length 
it is free to reach the surface ; if less, energy is required. For a complete 
dislocation in a close packed lattice b is the distance of a closest approach, 
equal to a/\/2 in the cubic lattice and a in the hexagonal. If the 
dislocation is moving slowly and is a few atoms away from a free surface, 
its energy is about 36D per unit length, where D=bu/27(1—o). 

This is compared with by in the following table. The value of y is 
that found experimentally for the liquid metal. Fisher and Dunn (1952) 
give the following ‘fairly reliable’ values for the surface energies of 
solid metals (units of dyneem): Ag 1180, Cu 1430, Au 1550. 
These lie 30-50% above the values for the liquids. . 


rs o y 3bD by 
Metal 101 dyne 10? dyne 10-erg 10-5 erg 
l0=§cm em em-t em-! om 

Ag 2-88 2-8 0:37 8-0 2-94 2°31 

Au 2-88 2-8 0-42 5-8-10-0 3:18 1-70-2-88 

Cd 2:97 2:1 0:30 6:3 2°11 1-88 

Cu 2-55 4-4 0:34 11-0 3°45 2-81 

Pb 3:49 0-6 0-44 4-4 1-04 1-29 

Pt 2°77 6-3 0-39 18-2 6-30 5-04 


Zn: 2-66 3°6 0:23 7:4 2-63 1-98 
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Except in the case of Pb, the energy of the dislocation is always greater 
than that of the exposed surface. An isolated partial (§ 5.21 of Part IIL) 
has 3 the energy of a complete dislocation, and exposes a ledge 1 /\/3 times 
as wide. The value of by now exceeds that of 40D in all cases. Such a 
crude calculation does not justify quantitative conclusions, but it is 
clear that a dislocation may require energy in order to break through 
the surface of even a perfect crystal. 

A screw dislocation does not expose a complete ledge as it enters or 
leaves the specimen, but gradually exposes a ledge as it moves. There is 
therefore no energy required to form or destroy it, but each end is acted 
on by a force by tending to shorten the exposed ledges. 


6.22. Hollow Dislocations and Capillary Equilibrium 
Frank (1951 b) has examined the condition that a.dislocation should 
have less energy when its core is hollow than when it is filled with strained 
material. If a screw dislocation of Burgers vector b lies along the axis 
of a cylinder of radius R, the elastic energy of unit length is of the form 


E,=(b'u/47) In (R/O). . 2 2... (TA) 


The value of 0, which is of order unity, depends on the law of force in the 
core of the dislocation ; Frank considers various laws. If the dislocation 
is hollow, with radius r so large that Hooke’s law applies accurately 
beyond this radius, there is an elastic energy (b7u/47) In (R/r), and a 
surface energy 27yr, where y is the surface energy, giving a total energy 


H,=(b'u/477) In (R/r)+ 27yr. ees bis ROY Ey) 
As r varies, #, passes through a minimum for 
r=fo=b"u/8n*y. air -werse Pletosete (116) 


The quantity y/u, which is a measure of the range of intermolecular 
forces, varies rather little from one substance to another, and is about 4 A. 
Thus r9>6 only if 6>20A, and unless this condition is satisfied the 
dislocation will certainly not be hollow. If it is satisfied, the dislocation 
will be hollow only if the minimum value of #, given by (115) and (116) 
is less than H,, which requires that b>8m7°6y/u. Since usually 0<1, 
this condition is satisfied if r»>b. We conclude that dislocations will 
be hollow if b exceeds about 20A (Frank obtains 10A as a more reliable 
estimate), but not otherwise. Verma (1951) has found hollow dislocations 
in carborundum (b>=15A). In a long chain paraffin (b5=47 A) Dawson 
and Vand (1951) found no hollows. Frank has pointed out that in the 
latter case the energy of the core of the dislocation is very low, because 
two long chain molecules have many positions with energy only slightly 
higher than that of their proper positions in the crystal. This means 
that @ in (114) is large, and the condition b>87°@y/u is not satisfied. 
Frank also considered the way in which a hollow dislocation widens 
out when it meets a free surface, and the shape of the depression in the 
surface where an ordinary solid-cored dislocation emerges. If the 
surface energy is independent of orientation, the depression is shallow 
in form, but of infinite volume and depth. If the surface is a natural 


Z2 
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habit face of the crystal, so that a slight change in the orientation of the 
surface exposed makes a first order change in the energy of the surface, 
the surface depression may be completely eliminated if 6 lies in the 
surface, and is small in almost all other cases. 


§7. THe Link ENERGY OF DISLOCATIONS 


The energy of an elastic dislocation of arbitrary shape may be 
calculated by the methods of Volterra and Burgers. We shall illustrate 
this in § 7.1 by considering an infinite body containing a dislocation of 
which the singular line is a circle of radius R. The body is made doubly 
connected by a toroidal hole of radius «<# enclosing the singular line, 
and the Burgers vector b of the dislocation lies in the plane of the circle. 
This corresponds to Orowan’s dislocation of fig. 4. In this treatment, 
the energy increases indefinitely as ¢ is diminished, and it would be more 
satisfactory to have a treatment based on an approximation such as that 
of Peierls. As will be seen from §§ 2 and 4 of Part III, these closer 
approximations have been applied only to a straight dislocation and to a 
dislocation which deviates only slightly from a straight line. 

Most theories of plastic deformation are based on models in which a 
dislocation is formed and travels through the lattice under the influence 
of an applied stress until it meets some type of obstacle. Either in the 
process of formation or in overcoming this obstacle, the dislocation can 
act only with the help of thermal vibrations. Different models give 
different dependences of the necessary activation energy on the applied 
stress, and an analysis of creep tests may indicate which model is 
appropriate to a particular type of deformation. In $7.2 we calculate 
the activation energy for a number of models. A related problem which 
is still unsolved is the calculation of the appropriate mean stress which 
opposes the motion of a dislocation through a lattice (such as that of a 
metallic solid solution) containing stresses which fluctuate over very 
short distances. Two approaches, leading to different results, are 
discussed in § 7.3. 


7.1, THE ELASTIC ENERGY OF A CIRCULAR DISLOCATION 


We take the origin of Cartesian coordinates at the centre of the 
dislocation of § 7, with the z axis normal to the plane of the circle, and 
the w axis parallel to the slip direction. Then the Burgers vector is 
b=(b, 0,0). Let the polar angle @ of a point (£,7, 0) on the circle be 
given by tané=7/f. Then the displacement u at any point is given 
by (35) of § 3.4. 


We have 
u*= (Q/47)(b, 0, 0) 
2 
ute i (0, 0, b cos 0)R do/4nr 
0 
Abita fe 
and —— motes! 
pb Im(A+2n) J bz cos 0 R dé/r, 


where r’=(a—R cos 0)?+(y—R sin 6)?+-22. 
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The elastic energy of the dislocation is 
W=—4) || p,, da dy, ere oe oa (LAT) 


where the integral is taken over the circle x?+y?<(R—c)?, z=0. 
Since p,,=(du,/dv+du,/dz), the energy may be written 


W=— 5-6 [[S dex dy— = = Rb [fo (S*) ao ae ay 


Measles z COS 
apne Re [|| cas (SS ) a0 ae dy. mies danke 


We consider first the integral involving 2, which is equal to the 
scalar magnetic potential of a unit current flowing in the circle of 
radius R. It is easily seen that the integral is the mutual inductance of 
the circles of radii R and R—e, taken with negative sign, and this mutual 
inductance may be evaluated as the line integral of the vector potential 
of the larger circle, taken round the smaller circle. This is 


an cos 6 Rdé 
9 [R?+(R—c)?—2R(R—e) cos 6]1/?° 


—27(R—e) 


To evaluate this when «<f, put sin 46=t, when the integral becomes . 
US ef dt 
In the region 0<t<(«/2R)!/* the integrand is approximately (<?+ 4R7t?)-1?, 
while in the region (¢/2R)'*<t<1 the integrand is approximately 
(1—2¢?)/2Rt(1—#)1/?. 
Evaluating these two separately, we obtain a contribution to the 
energy of 


— 87k? fac 


W ,=Rb*u(4 In R/e+ In 2—1). 
The integrands of the second and third terms in W reduce in the 
plane z=0 to —(x—R cos @) cos 6/[(~—R cos 6)?+(y—F sin 6)? )?”. 
Integration with respect to x gives cos 6/[(7—R cos 0)?-+(y—R sin 6)?}!”, 
taken between the limits x= +[(R—<«)?—y?]!?. Write x=(R—e) cos «, 
y=(R—e) sin «, and the integrals become 
1 


2 
a | Whe | ay | ea cos (—a) }1? 


Pennants see ar ouee ] cos 6 cos a dé da 
~ {RP (R—e)?+ 2R(R—e) cos (6-+-«) }3/ 
a cos ¢ dh 
ea tires) i [R?+ (R—e.)*—2R(R—e) cos ¢]??° 
The contribution of these terms to the energy is thus [(A+)/(A+2)—3] 
times the contribution of the first term, giving a total energy of 


W=2rR wel a+ a (In R/e+2Mn 2/e), . . . (119) 
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The factor in square brackets is }[1+1/(1—o)], the mean of the 
corresponding factors for a straight screw dislocation and a straight 
edge dislocation. To estimate the energy of a circular dislocation in a 
cubic crystal lattice we set «=2¢/e, where (=b/2(1—a), by analogy with 
the value found in § 2.2 of Part III for a pair of straight dislocations. 


7.2. ACTIVATION ENERGY PROBLEMS 


In most models of the creep process, a dislocation will travel without 
any activation energy if the applied stress reaches some value Pj, 
the critical stress. Under a lower applied stress P, it will move only with 
an activation energy U(P,, P1). Since U=0 when P)>=P,, it is convenient 
to express U as a function of P, and the activation stress P;—Py. If we 
make the reasonable assumption that in similar tests at different 
temperatures 7’ the creep rate always becomes noticeable when 


U(P,—P,, P;)=mkT, i eee 


where m is a constant independent of P,, P, and 7’, and k is Boltzmann’s 
constant, we see that creep experiments at different temperatures give 
some information on the form of the function U. 

Wyatt (1951) has shown that if creep occurs by a pure exhaustion 
process (Mott and Nabarro 1948, Smith 1948), only dislocations with 
a narrow range of values of P, contribute to the creep rate at any given 
time ¢ after the start of the test, and that the activation energy U of 
the active group obeys the law 

dU kT ae 

Fiekr re 
Suppose the specimen is allowed to creep under a load Py, up to a 
time ¢,, and the load is then slightly increased to Py, The activation 
stress and activation energy for the group of dislocations active at this 
time are decreased, and the creep rate increases to the larger value 
which it had at some earlier time ¢;. We now assume that U depends 
very little on P, in the range considered, and that the total number of 
dislocations available in a range from P, to P,+dP, and also the total 
amount of slip caused by each dislocation when it moves are almost 
independent of P,. It follows that the activation energies, and hence 
the activation stresses, are the same when a load P, is applied after the 
specimen has crept for a time ¢, under a load P, and when a constant 
load P, has been applied for a time ty. 


Thus P(t) —Pog=P (ts) —Prow 
or Pi (le) — Py (tg) = Pop — Pes on ead ee ee ee 
a known quantity. It also follows from (121) that 

U(t,)— U (ts)=kT'(In t,—In ty). » «, . ©... (128) 
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Comparison of (122) and (123) gives the form of the function U(P,—P,). 
The difficulty in applying the method in practice is that the experi- 
mentally available range of P, is so small that the function is almost 
linear in this range. 

We now calculate the function U(P,—P,, P,) for a number of models. 
In § 7.21 we take account of the fact that the line energy of a dislocation 
is not strictly constant, but in other cases we assume a constant energy H 
for unit length. 

7.21. The Models of Becker and Orowan 
_ The original model of Becker (1925) and Orowan (1934 a) does not 
provide a detailed mechanism for the nucleation of slip, but assumes 
that slip will begin if the shear stress in a volume V reaches a critical 
value P,. In the absence of an external stress this requires an energy 
VP,?/2u, where w is the shear modulus. In the presence of an external 
stress Po, this is reduced to 

U=V(P,—P )*/2p. al oh se sah cee are ae 

If observable slip begins when U=mkT, the strength at temperature 7 
is given by 

eS ATUL ede A ee se i ie (125) 

We may also consider the specific mechanism proposed by Orowan 
(1934 b), in which the nucleation of slip occurs by the formation of a 
loop of dislocation in a region of the crystal previously perfect. We 
assume that the region in which the loop is formed is determined by the 
presence of a residual stress P* in the same direction as the applied 
stress Py, and also of a stress concentration factor g. The total stress 
on the loop is P=P*+qP,. When the radius of the loop is &, the total 
energy of the specimen is W—7-h?b(P*-+qP,), where W is given by (119). 

In the next approximation, the circular dislocation cannot. be in 
equilibrium under the stress P, for this stress exerts an equal outward 
force bP on each element of its length. Since the energies per unit 
length of the edge and screw elements are in the ratio 1: 1—o, their 
radii of curvature must be in the same ratio. The equilibrium shape of 
the dislocation is roughly elliptical, with the minor axis in the glide 
direction, and shorter than the major axis by a factor (l—o)/*. If we 
assume the circular form, the total energy of the specimen increases by 
W—7xR*bP, which is positive for small R, reaches a maximum, and 
becomes negative for large R. The radius R,, of the circle of maximum 
energy, and its energy U, are given by 


a 1 1 Oe war 
_= — a eet 126 
Be 8a a ers )me c a 
U1 Ly he (228, 
and ma ( i+ —) 2 (m= - ). cy testo (127) 


Numerical values have been obtained by Frank (1950a). He replaced 
our factor of 1+1/(1—c) by 2, and considered the values (=2b/e and 
~=4b/e where we take (=b/2(1—o). His results are shown in fig. 31. 
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Since 63. is of order 5 ev, and thermally activated processes occur 
very rarely at room temperature if U is greater than about 1} ev, loops 
will be formed spontaneously only if P/w is close to the critical value, 
which is about 1/20 or 1/30. For a soft metal such as sodium, where b®u 
is less than 1 ev, loops may be formed when Pj is about 1/80. 

7.22. The Model of Mott and Nabarro 

Mott and Nabarro (1948) considered a straight element of dislocation 
of length L moving in a periodic stress field, so that the total force on the 
dislocation when it has moved a distance x is ; 

P,bL cos (272A). Oe Ay Bae ene 


Fig. 31 


Energy of Critical 
Dislocation loop 


C= 2be 
U 
Bu = 4b/e 
Reciprocal stress 

Radius of Critical 

Dislocation loop 
R 
id. (6) 


Reciprocal stress 


(a) Energy and (6) radius of a metastable loop of dislocation as a function of — 
the applied stress, for two values of the energy in the core of a dislocation 
(after Frank). 


An applied stress Py exerts a force P,bL on the element, which thus has 
stable and unstable positions of equilibrium given by alternate roots of 
the equation 
P, cos (27a/A)=P,. ON er bate EL) 
If Py is nearly equal to P,, a pair of roots is given by 
2rra/A=+[2(P,—Pp)/Pi}*. . . . . . (180) 
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The activation energy is 


2/2, (P,—P,)8? 
v= oa Men h(1St} 


and this is equal to mkT’ when 
Po=P,—(37mkP/2,/2bLA)T2®, . | . . (132) 
7.23. The Models of Orowan, and of Frank and Read 
Orowan (1948) proposed a model of an over-aged precipitation- 
hardened alloy in which the nucleation of slip occurred by the bursting 
of a loop of dislocation held fixed by particles of precipitate at two 
points a distance 2R apart. Frank and Read (1950) introduced a model 
which would give an unlimited amount of slip from a single source. 


Fig. 32 


) 


A stable circular are of dislocation AMB and an unstable arc AM’B, each of 
radius 7 and closing a gap AB of width 2R. 


The mechanism by which a nucleus of slip starts to grow is the same 
as that suggested by Orowan. This type of nucleus has the remarkable 
property that its behaviour is determined solely by the one parameter hi. 
Assuming that the line energy £ is constant, the critical stress P, is 
given by 
Ae hy Die RA ole 2. (133) 
For this stress, the loop of dislocation in equilibrium forms a semicircle 
of radius R. At a smaller applied stress P, there are two positions of 
equilibrium (fig. 32), both circular arcs of radius r given by 
reatiOleas we arta. ee ( 134) 
The shorter, stable arc subtends at its centre an angle 2 sin“'(ft/r). 
Its line energy is 2Hrsin-“(R/r). As it has moved from its original 
position as a straight line joining the two points to its present poston 
the applied stress P, has done work 6P,|7? sin“! (R/r)—R(r?—#*)'"]. 
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The total energy of the stable arc is thus 
2Er sin-! (R/r)—bPolr? sin (R/r)—R(r?—R*)"”]. 
By means of (133) and (134) this may be reduced to the form 
2 2)1/2 
a [sin 5 + re | hit ee Pee as 
The line energy of the larger, unstable, arc is 2H7r[7—sin™ (R/r)], 
while the work done in its formation is 
bP [r2{a—sin- (R/r)} + R(r?—R?)?]. 
The total energy reduces to 
EH? Po el eee] 


a7—sin-! 


OP, P, PP 
The activation energy U is the difference of (136) and (135), given by 


(136) 


2H? PRS CE E  is p aie 
Eee Pe eee 
U= BP, cos P, PE |: (137) 
If P, is very nearly equal to P,, this becomes 
2(pP._ P.)s/2 
_ 8 2B? (Py— P,P? re 


3b Pe 
This model, like that of Mott and Nabarro, leads to U «(P,—P,)*. 
The essential reason for this is that Mott and Nabarro take the length of 
dislocation which is moving to be a constant L, while in the present model 
the length again tends to a finite value 7R as Py, tends to P,. 
If slip is observed when U=mkT, we have 
Po=P— (8bmk/84/2E?)? 8P P8728, in.) eee 
Thermal agitation allows slip to be observed when P, is less than P, 
by a fraction . 
Pi—P, _ (3bPymkT\2!8 
Por \e8a/ ar 
Taking as a lower limit for H the value 0-16 63. and assuming m=40 


leads to 

Pi Po sg P, kT\28 

Pl \ ee 

If we take P,=lkg mm, p=4x10"dynecm, b=3x10-® om, 
T'=300° k, this becomes 

(Py—-P,)/P,=55 X10" iS Lae ee 
Thus, on this model, the strength should be very insensitive to 
temperature, and the assumption that P,—P, is small is always justified. 
7.24. The Model of Cottrell and Bilby 

In the model of Cottrell and Bilby, a long straight dislocation 

approaches a region of adverse stress, similar to that considered by 
Mott and Nabarro, In the application which they considered, a 
dislocation is bound to a row of solute atoms by a force which for small 


(140) 
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displacements increases linearly with the displacement, and with 
increasing displacement reaches a maximum and decreases to zero at 
large distances. We consider only the behaviour of the dislocation 
when it lies near to the region of maximum restoring force. The 
nucleation of slip occurs by a bulge of dislocation spilling over the line 
of maximum stress by thermal activation. There is a form of bulge 
which is in equilibrium at all points, the equilibrium at some points 
being unstable. Once this is formed, further slip can occur without 
activation. ‘The width of the critical bulge depends on the activation 
stress in a manner to be determined. 

The analysis is tractable only if we assume that the critical bulge 
always has a small slope. That is to say, if x is the coordinate in the 
glide plane along which the internal stress is constant, and y is the direction 
of motion of the dislocation, 


dy|dx <1. Raat ar are neers pr antis tte ge | 


The curvature of the bulge at any point is then d?y/dz?. We take the 
maximum value of the internal stress to be P,, while the stress at points 
near the maximum is P,—x«’y?. If an external stress P, is applied, the 
form of a dislocation which is in equilibrium at all points is given by 


Ed?y/da?=b(P,—Po—kKy*).  . 2. « ((148) 
Two solutions of this are obvious. The first, 
Y=—Yo=—(Pi—Po) lk, . . » « (144) 


is a straight dislocation in stable equilibrium at all points, which has 
not crossed the region of maximum stress. The second solution, y=+Y, 
represents a straight dislocation entirely: in unstable equilibrium. The 
critical form coincides with (144) at points far from the centre of the 
bulge, and has y>0 near the middle of the bulge. To integrate (143) 
we multiply both sides by 6 dy/dx, and fix the constant of integration by 
the condition that dy/dz=0 if y=—y). The result is 


3H (dy/dx)? = 6b(P,—Po)(ytYyo)—2b(y? +40"). - + (145) 


In general, the solution of an equation of this form is a Weierstrassian 
elliptic function. In this particular case the solution required is 
elementary, 

Yy=y,(3 sech* x/x,—1), ee ee ee: LO) 


where y= (4H? /b?x?(P,—Py)}"4, ees 32 a 14) 


as may be verified by substitution. The maximum value of dy/dx is 
4y,/\/3%, which is proportional to (P;—P))°*. As Po tend to Py, this 
tends to zero, and (142) is thus satisfied. 

The width of the bulge, a), tends to infinity as the»activation stress 
(P,—P,) tends to zero. We therefore expect that the activation energy 
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U will tend to zero with (P,—P,) more slowly than in the models of §§ 7.22 
and 7.23. The line energy of a bulged dislocation exceeds that of a 


straight dislocation by 


For the critical bulge (146) this is 


(18H y 92/297) [ sech4(a/x) tanh? (a/a) du=24Hy,"/5ay. . (148) 


The work done by all the stresses, internal and external, on an element of 
the dislocation as it moves forward is 


bf! (Py—Pr teen?) dqp=b[(Po— Ply + yo) + HOV + U0): 


Integrating over the length of the dislocation, we find the total work 
for the critical bulge to be 


— 9bK?y,° { sech4(x/x) tanh? (a/x,) dx 


== — 1260 y,°/5—=— 24 by /5ay. ss 
The activation energy is the difference of (148) and (149), 
48Hy,2 24 Gay 
0 = SS 
5X 5 \k 
If slip becomes observable when U=mkT, the strength is given as a 
function of temperature by 
Po=P,—(25m7k?«3/1152Hb)? 5745, a Rl ye ge ee 
In the case considered by Cottrell and Bilby, the stress is given in 
terms of the displacement y’ from the position of tightest binding by 
P=2Apy'/(y’+p?), og ely 45 4 Oc Un 
where p is a length estimated to be 2x 10-§ cm. The maximum stress P, 
occurs when y’=p/4/3, and is 
P,=3/3A/8p?. Mr a eiame es toe SBN 
At this point dP/dy’ vanishes, and d?P/dy”=—9P,/2p?. If we 
approximate to this law by a parabolic law giving the same maximum 
stress and the same value of d*P/dy? at the point of maximum stress, we 


find 
Ke QP /4p*. ce tn Se See Ee 
The best estimate of P, is (Cottrell 1952 a) P,=6 x 101° dyne cm. 
Taking »=8-28>10! dyne cm, b=2-48x 10-8 em, we estimate H= 
0:166%1=8-13 x 10-° dyne. If we assume m=50, eqn. (151) becomes 
Po=6 X10°— 4K 108 TS a 
where Po is the resolved shear stress in dyne cm (half the yield 


strength for a. crystal with a glide system favourably oriented), 
and 7’ is the absolute temperature in degrees Kelvin. At room 


PoP, 84. ee 


Catia at cnt 


a Se ee 


+ a Ou” ee 
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temperature, 7=300°K, and P)=2-2x10!dynecm-*. At 80° x, 
Po=47 x 10! dyne cm, which is already 0-78 of its limiting value at 
low temperatures. The present approximation, which assumes that 
dy/dx<1, is always valid when P,>P,. The greatest value of dy/dz, 
namely 4y9/4/3a , occurs when Py=0, and is 4 (bpP,/E)'?=0-8. Even 
in this limiting case the error is not likely to be very serious, since the 
greatest value of dy/dx is multiplied by the vanishing factor d2y/da. 


7.25 A Model of an Over-aged Alloy 

Orowan’s model of an over-aged alloy considers a dislocation held 
rigidly at two points, and a bulge forming and bursting between the 
points. This model should apply to an over-aged alloy in which the 
particles of precipitate have become incoherent with the lattice of the 
matrix, and there are no longer any internal strains. Where the misfit 
between the lattices of the matrix and the precipitate is small, over-ageing 
may occur by particles of precipitate growing to a large size, while still 
remaining coherent. Dislocations moving in the resulting stress field 
will move most freely in the regions of the lattice furthest from particles 
of precipitate, but they are not held at fixed points close to the parts 
which move. We are thus led to consider a model in which a dislocation 
moving in the y direction is held up near the line y=0 by an adverse 
internal stress: this adverse stress is somewhat lower in one region of 
this line (which we take to lie near x=0) than elsewhere. The exact 
form of the field of internal stress cannot greatly affect the general 
physical behaviour, and our choice is governed by mathematical 
convenience. The way in which a suitable function was found is indicated 
at the end of this section. We actually take the internal stress to be 


Py —12y? + (4B? /2b?x,)[1—}(9—y2) sech! (x]x,)], . . (156) 
where y and z, are arbitrary constants. 
Assuming as before that dy/dx<1, the shape of the dislocation is 
given by 
Ed®y|da?=b{P,—P )—ry?+ (4B?/x7b?a,*)[1—4(9—y”) sech4 (x/a1)]}. (157) 
We now look for approximate solutions of this equation which will be 
exact in the limit P,>P,. The solution corresponding to (144) is 
Y= (2E]«*ba")[4(3—y) sech? (a/a,)—1], . . . . (158) 
while the solution corresponding to (146) is 
Yo (2H /x*bx,?)[$(8+y) sech? (x/z,)—l]. -« . . . (159) 


If y=3, the dip in the potential ridge vanishes, and (158) and (159) 
reduce to (144) and (146). In general, the critical stress is not equal to Es 
and there is a positive activation energy when P»>=P,. There is no simple 
way of determining the actual critical stress, or of finding the activation 
energy when the applied stress is just less than the critical stress. In the 
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particular case y=0, the solutions (158) and (159) coincide, which means 
that the critical stress in this case actually is P,. For Py=P,, the 
equation 
Ed?y/dx?=b{P,—P,.—y? + (4E?/«*b*x,*)[1—2 sech* (x/x,)]} (160) 
has two coincident solutions . 
y*= (2H /x*bx,”)[8sech? (x/x,)—1]. She ae a ee 
If the applied stress Py is now taken slightly less than P,, the solutions y* 


separate into two solutions again representing the stable and the unstable 
equilibrium forms of the dislocation. These solutions may be written 


yy hay ten Ce ee 
where 7 is small. We substitute (162) into (160) and find that 7 satisfies 
the equation 
En" +-212by*n=b(P,—P y—x*y*),” = >--- - « (163) 
or, if we write 
gest, PO ee ee 
d’nid&’-+ (6 sech? £—4)n=bx,?(P,—Po—k’n’)/H. . . . (165) 


Since the solutions y* coincide, they represent states of neutral 
equilibrium. We therefore expect that the dominant term in 7 will be of 
order (P,—P,)'”. This means that the terms on the left of (165) are of 
larger order than those on the right, and must vanish to order (P,—P,)!/?. 
Thus the dominant term in 7 is a solution of the homogeneous equation 
7” +(6 sech? €—4)n=0. This equation has the bounded, even, solution 
n=sech? €, and an unbounded, odd, solution sech? € [° cosht w du, which we 
reject. We thus try for 7 a solution of the form 


n=C(P,—Po)'? sech* €+(P,—Po)e(é), . . . (166) 


where the constant C and the function ¢ remain finite as Pyp>P,. To 
order (P,;—P,), eqn. (165) reduces to ' 


("-+(6 sech® £—4)¢=(bx,2/B)(1—12C? sech4 €). . . (167) 


Again the homogeneous equation has a solution ¢=sech? €, and we look 
for a general solution of the form 


C=u sech? €. 0 eu WAR oe UN el ae 
Then wu must satisfy : 
w” sech? €—4u’ sech? € tanh €=(ba,?/H)(1—x®C? sech €). 


This equation has the integrating factor sech* €, leading to the integral 
é 
u'=(be,2/E) cosh! é | (10? sech* y) sech? x dy. . (169) 


We are looking for solutions in which ¢, and therefore uw, is an even 
function of €. It follows that we must take €;=0. Moreover, for large 
| €|, ¢ must be bounded, w must be O [exp 2| |], and w’ of the same order. 
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The integral in (169) is elementary, and it is easily seen that w’~exp 4| é | 
unless «°C?—15/8. This choice leads to 
u’=— (ba,?/8H)(4 sinh € cosh +3 tanh €). 
Integrating again, we find 
U=— (bx,?/8E)(2 cosh? €+3 In cosh €)-+u, | 


C= — (bx,?/8H)(2+3 sech? é In cosh €)+-w, sech? &. 
Since we have postulated that ¢ remains finite as P,>P,, the last term 
in (170) is of the same form as the dominant term in 7, but of smaller 
order. We therefore neglect it, and obtain the pair of solutions 

n= +[15(P,—P)/8x?]!/2 sech? € 

—[bx?(P,—P)/8H](2+3 sech? € In cosh €). Se Ae ee LT) 

The two solutions differ by a bulge of width about 2x,, which moves 
forward a distance [15(P,—P,)/2x?]!?. This behaviour is of the same 
form as in the models of §§ 7.22 and 7.23, and the activation energy is 
again proportional to (P,—P,)?”. 

We now indicate the argument by which a suitable form of stress field 
was found. In the limit Py»=P,, the shape of the dislocation is governed 
by an equation of the form 

ee oe eee eae he ee ee ey 
where X and Y are constant multiples of x and y respectively. We wish 
to know two solutions of this equation. Letthem be Y,and Y,+ 4. Then 


(170) 


Veet Meena ee ooo (178) 

and eo yeteay Ao At A Xy.1 bln (174) 
Subtracting (173) from (174), we find 

Vo (ee aA ia A ete Pal ana (175) 


Thus if 4 is known, Y, may be determined explicitly and without 
quadratures. Substituting this value of Y, in (173) similarly determines 
f(X). The function 4 measures the shape of the bulge which moves. 
forward as the dislocation moves from its stable to its activated position. 
It is convenient for this bulge to be even in X and limited to the region 
where X is small. Moreover, we know from (144) and (146) that in the 
casef(X)=2, eqn. (173) hasthesolutions Y;=—1and Y,+ 4=8sech’ X —1, 
corresponding to J4=3 sech? X. We are thus led to consider the problem 
arising from the choice 4=y sech? X, which leads to a stress field of the 
form (156). 
7.26. A Model of a Localized Barrier 

Professor Cottrell has suggested to me a class of models* which differ 
essentially from those already considered. The critical stress P, is large, 
but acts in a very small region, so that the activation energy may be 
of order mkT even if the applied stress Py is much smaller than P,. 
As an example we may consider that the active glide plane is threaded 
by stationary dislocations corresponding to another glide system. 


* See Cottrell (1952b). 
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Suppose there are no long range forces between the two sets of dislocations, 
but that an atomic derangement of high energy is produced in a small 
region with linear dimensions of order 6 if the lines of both dislocations 
pass through this region. As a specific model, we may suppose that the 
two dislocations are both of pure edge type, and would pass freely through 
each other, but for the fact that each is dissociated into two partial 
dislocation (§ 5.21 of Part III). The four partial dislocations all have 
screw components, and cannot pass through one another pair by pair, 
The complete dislocations can only pass through one another if their 
components recombine over a short distance near the crossing point. 
This idea is due to Heidenreich and Shockley (1948). An activation 
energy U, is needed for this, and in the absence of thermal vibrations a 
moderate stress P, will not appreciably alter the separation of the 
components of the two dislocations. 

Suppose a dislocation is held up in this way at a succession of points 
separated by distances 2R. The position of stable equilibrium depends 
very little on the applied stress Py. If the activation energy U is supplied 
at one of the locking points, the dislocation at this point moves forward 
to a position of unstable equilibrium, a distance 6b ahead, where @ is of 
order unity. The loops of dislocation joining this point to its neighbours 
each move forward by an average distance }6b, and the applied stress 
does work 26b°>RP,. The activation energy U for the stressed system 
is thus 


U=U,—200°RP 4. & ay oe ee 
If creep is first observed when U=mkT’, we have a strength given by 
P =U, |26b?R—(mk/26b?R)T. ~ sae ae 


If we compare two similar systems with the same value of U,, but 
differing spacings R and R’ of the locking dislocations, they will begin 
to creep at the same temperature under different stresses Py) and P,’, 
where 

RPU=R Py ew sh ne 


When the linear relation (176) holds between U and P, incremental creep 
tests (§ 7.2) give directly the quantity dU/dP,=—20b?R. It follows 
from (178) that in two such tests at the same temperature on different 
specimens the quantities P)dU/dP,) and P,' dU'/dP,’ must have the 
same value. Wyatt's observations (1950) satisfy this relation well. 


7.3. THE AVERAGE STRESS ON A DISLOCATION WHICH IS NEARLY STRAIGHT 


Consider a dislocation lying in a medium in which the internal stress 
varies from place to place about a mean value zero. Suppose the average 
amplitude of the varying stress is P,, while the distance between 
neighbouring regions of maximum stress is A. The stress P, will bend 
the dislocation into ares of radius p given by 


p=E/bP. 
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If p <A, the internal stresses can bend the dislocation so that it lies almost 
entirely in regions of low potential energy where the internal stress is 
much less than P;. If an external stress is applied, the dislocation moves 
forward by a process in which individual loops of length A move forward 
by a distance J from one potential trough to the next. The analysis 
of § 7.22 was developed for this process. No quantitative analysis has 
been attempted for the case p==A, which seems to represent a 
precipitation-hardened alloy treated to develop maximum hardness. 
We now consider the case p>A, which may be written H>bP.A. Since 
F is of order 46*u, this inequality may also be written 


ATS Dy Dipole Dil wr elnic phe MT 


When (179) is satisfied, the internal stresses are not large enough, in 
comparison with the tension of a dislocation, to force the dislocation into 
the regions of low stress. Lengths of the dislocation large in comparison 
with A run almost straight across the crystal, being bent by individual 
regions of stress into arcs with radii of curvature p=E/bP,. The 
amplitude of the wavy curve into which the dislocation is bent is of 
order A?/p=2A?P,/bu. 

If the dislocation showed no large scale deviations from straightness 
it would move under the influence of any external stress, however small. 
For a dislocation of a length ZL which is much greater than A is composed 
of L/A elements, on each of which the internal stress exerts a force of 
order of magnitude bP.A and of random sign. The resultant force is of 
order bP.A(L/A)1’?, while the force due to an external stress P is bP.L, which 
can always be made larger than bP,A(L/A)! by taking L large enough. 

In fact, we must not assume that the mean path of the dislocation 
will not depart from a straight line, however large L becomes. We can 
find a length L so large that a dislocation of this length is effectively 
flexible under the internal stresses. The motion of a dislocation under an 
external stress takes place by the movement of individual segments of 
this length L. No rigorous way of estimating L has been discovered. 
Two plausible methods have been suggested, and these lead to essentially 
different results. Masing (1948) has given an exact treatment for an 
extremely simplified model. 

Mott and Nabarro (1948) pointed out that the mean force on a 
length L, equal to bP,(A/L)'?, would bend the mean position of the 
dislocation into a curve of radius 


BLN bp (L\1? 
Piss bP, (7) SIF (3) . . . . . (180), 


If L is so large that this radius p, is comparable with L, we may say that 
the dislocation is effectively flexible under the internal stress. The 
condition for this is 
E2 bu? 
E= PIA 4PA? 


P.M. SUPPL.—JULY 1952 2A 


(181) 
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and the average stress on a length L is 


(=) PO i (182) 
L : E bu. 

The same argument may be put in another form. Each segment of 
length A is bent into an arc of radius E/bP;,. Its ends are therefore 
inclined at an angle bP.A/E. The average angle between the mean 
directions of the dislocation at points L apart is (L/A)!*bP,A/E. If this 
is of order one radian, we again have the condition (181). 

Mott (1952 b) has suggested another way of estimating the length L 
of the loop of dislocation which moves independently. He assumes 
that neighbouring positions of equilibrium for such a loop will be a 
distance A, apart near the middle of the loop, while coinciding at the 
ends of the loop. We may take these positions to be circular ares bulging 
a distance 4A to either side of the chord of length L joining their ends. 
Then L?=4p,A, giving 


4B \203 Qby\28 
a pees {3 —~ [| — 1/3 
ra (Yon (Ye, 5 89 


and an average stress 


AN 1/2 AbP, 1/3 a AP, 1/3 


Part III. Errect oF THE CRYSTAL STRUCTURE 


§1. Exastic DisLocATIONS AND CRYSTAL DISLOCATIONS 


The theory of elastic dislocations can be applied to determine the 
stresses far from a dislocation line, and to establish general theorems, but | 
there are important differences between elastic dislocations and 
dislocations in crystal lattices. The elastic dislocation is made in a 
multiply connected body containing certain fixed holes. If we try to 
continue the solution into the interior of these holes the strain becomes 
very large. For a given radius of hole the strain in the medium can be 
kept less than some predetermined small value ey, by reducing the 
vector b. In the crystal dislocation, an analysis of atomic displacements 
is required at distances of the order of one interatomic spacing from the 
dislocation line, while the vector b must be a lattice vector. Estimates 
of the strain are therefore required where the strain is no longer small. 
It is an essential part of the picture of crystal dislocations that for these 
finite strains Hooke’s law is no longer valid, and in particular that a 
relative displacement of two layers of atoms by a distance b produces 
zero restoring force. Because of this, and because the toroidal holes of the 
elastic theory are reduced to lines threading the lattice, a crystal 
dislocation is mobile. 
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In defining the Burgers vector b of a crystal dislocation, it is convenient 
to use the method of counting steps round a closed circuit which was 
described at the end of § 2.3 of Part IT, instead of the circuit integral (17) 
at the beginning of that section. In making this count, no attention is 
paid to the fact that the steps are not all precisely parallel. The circuit 
must lie entirely in regions of moderate or small strain where local 
crystal axes can be identified. The local axes at one point of the circuit 
may be rotated with respect to those at another point, and the labelling 
of the individual axes is determined by continuous comparison along 
the circuit. The result is expressed in lattice vectors of the lattice at 
the arbitrary initial point. It will be seen that the residue of steps is —b, 
for an increase in strain increases the length of each step, and decreases 
the number of steps. This method, which is based on the fact that the 
displacement in a crystal lattice is only well defined at each lattice point, 
is not only more convenient, but also more consistent. In an elastic 
dislocation the idea of a continuous distribution of dislocation lines is 
meaningless, as each line must be enclosed in a hole, but a continuous 
distribution of dislocations arises in some theories of crystal dislocations 
(e.g. § 2.1). For such dislocations the value of the integral (17) is not the 
same for all circuits enclosing the ‘ middle’ of the dislocation, whereas 
the count has always a residue of a whole lattice vector, equal to the value 
of the integral taken round a very large circuit. Frank (1951 a) and 
Read and Shockley (1952) have given detailed discussions of the 
definition of the Burgers vector of a crystal dislocation. 

The general connection between elastic dislocations of the first three 
orders and crystal dislocations is apparent, and we have seen that 
dislocations of the last three orders do not lead to simple glide when they 
move across the body. What are the properties of a crystalline body 
which contains a dislocation of the latter group? The crystal lattices 
on either side of the cut surface are inclined to one another. The surface 
is thus a definite discontinuity in the structure. For isotropic media, 
and also for dislocations of the first three orders in crystalline media, the 
surface is arbitrary so long as its boundary remains fixed, but this result 
is no longer true for dislocations of the last three orders in a crystal. 
For example, if the ring of fig. 1 is crystalline, its final state will depend 
on the orientation of the line AA’ with respect to the crystal axes. Just 
as the rotation of Y which generates the dislocations of orders 4, 5 and 6 
may be regarded as composed of displacements of its elements proportional 
to their distances from the origin, so the dislocations may be regarded 
as continuous distributions of dislocations of order 1, 2 and 3 over 2. 
Thus it is apparent from eqn. (28) of Part II that the dislocation d,,, is 
closely related to a uniform distribution of infinitesimal dislocations 
b, and b,. This analysis takes a slightly different form in the crystal 
lattice, the continuous distribution of infinitesimal dislocations being 
replaced by a grid of localized dislocation lines, each of strength one 
lattice vector (cf. § 2.31). 

2A2 
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Bilby (1950 c) has published some instructive photographs of models 
of dislocations in a simple cubic lattice. We shall discuss later (§ 5.21) 
the relation between the dislocations in a two dimensional raft of bubbles 
and the dislocations of a close packed crystal. 

According to § 2.1 of Part II, the strain at a distance r from a line 
dislocation in an elastic medium obeying Hooke’s law is of order r~?, 
and therefore infinite at r=0. If the energy of the dislocation is expressed 
as an integral over r, this integral diverges at both limits. The divergence 
at infinity is inevitable; that at the origin arises solely from the 
mathematical approximation of assuming Hooke’s law to hold for large 
strains. In fact, if one plane slides over another, the shear stress between 
them must be periodic with period b. Approximate solutions have been 
derived for the displacements in a body containing a dislocation on the 
assumption that the shear stress across the glide plane in any region is 
a sinusoidal function of the amount of glide which has occurred in that 
region. The approximations are therefore only valid if the relative 
displacement of points above and below the glide plane changes by 
amounts much less than b in passing from one pair of atoms to a neighbour. 
Two types of approximation have been based on this idea. In the first, 
due to Peierls (1940), the material on either side of the glide plane is 
assumed to behave as an isotropic elastic continuum. The second type 
of approximation, due to Frenkel and Kontorova (1938), assumes the 
material below the glide plane to be rigid, while that above the glide plane 
consists of a single row or sheet of atoms connected by springs. The 
second model is only of illustrative value in the study of slip, but its 
properties are more easily studied, and it has other important applications. 
We discuss these approximations in §§ 2 and 3, while §4 deals with the 
more delicate problems which arise as soon as we take account of the 
atomic structure of the surfaces which face one another across the glide 
plane. Finally, § 5 treats problems which depend on the detailed structure 
of the lattice under consideration. 


§ 2. THE APPROXIMATION OF PEIERLS 
2.1, THE SOLUTION FOR A SINGLE DISLOCATION 


Figure 33 represents a crystal containing an edge dislocation lying 
along the axis Oy. The glide plane z=0 is taken midway between the 
sheets of atoms A and B which slip over one another. The natural 
spacings of atoms in the # and z directions are 6 and a respectively, and 
the material is considered to be elastically isotropic. In some applications 
of the model we shall assume a=b. This does not materially affect the 
analysis. The atoms of A in the neighbourhood of O are spaced at 
distances less than b. They are subject to two forces. The elastic 
reaction of the material above A tends to increase the spacing to 6, while 
the attraction of the atoms in B tends to reduce the spacing still further. 
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We calculate these forces separately as functions of the positions of the 
atoms, and determine the positions of the atoms by requiring the two 
forces to balance at all points of the glide plane. 

Since the dislocation is long and straight, the problem is one of plane 
strain. Let the displacements of an atom of A be (u, w), while those 
of a neighbouring atom in B are (@, @). The origins of these displacements 
are chosen so that when all the displacements vanish the atoms in 
A and B are their correct distance a apart in the z direction, but are 4b out 
of step in the x direction, as is shown near the point O. Then it is easily 
seen that i= —u, ©=w, so that A and B, though curved by the presence 
of the dislocation, still fit together without any stress p,, between them. 


Fig. 33 
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Peierls’s model of a dislocation. The material above AA’ and below BB’ is 
regarded as forming an elastic continuum. The force between the rows 
AA’ and BB’ is a periodic function of the displacement. 


To calculate the elastic force, we consider A and the rows above to 
form an elastic continuum, to the free surface of which is applied a 
tangential stress p,,. Then formal elastic theory leads to the relation 

fe du(x') dx" (1) 
Paz ea) vax" x’ —x’ 
where the Cauchy principal value of the integral is taken. 

Eshelby (1949) has given a direct proof of (1) using the Veal 

results for a dislocation in a continuum. The variation of u(x’) with 
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position x’ is considered to be caused by a continuous distribution of 
elementary edge dislocations along the # axis, with density 6’(«’) dx’. 
Since the relative displacement of the two surfaces is 2u, we have 
b'(v’)=—2 dujdz’. The shear stress which this infinitesimal dislocation 
produces at the point 2 on the glide plane is given by (9) of Part IT as 
b'p/27(1—o)(w—a’), and (1) follows by integrating over all a’. 

The stress given by eqn. (1) is to be equated to that caused by the 
atoms in B, which we assume to be a sinusoidal function of the relative 
displacement u—a—=2u, with period b. The amplitude of the function is 
determined by the condition that a small relative shift of A and B from 
their equilibrium positions produces a restoring stress equal to that given 
by Hooke’s law for the same shear. This leads to 


bp 
Paz — 5— Sin: fo ed jae 
Prez a 2m, ) (2) 
where d= 2u/b. 
Comparing eqns. (1) and (2), we obtain the integral equation 


© dd(x’) da’ l—o , ; 
[ He) = ane Choa ae 


The solution corresponding to fig. 36 is 


1 2(1—o)a 
jee bye reas 
7 a 
(4) 
= Stantts 
or Uu= on i 


‘where (=a/2(1—o). 

The width of the dislocation is small, for the displacement falls to half 
its extreme value at a distance ¢ from the centre. The original assumption 
that the dislocation is spread over a large number of atoms has thus led 
to a contradiction, and, if the sinusoidal law of force (2) and the neglect 


of stresses p,, are valid approximations, a dislocation must in fact be’ 


only a few atoms wide. We shall discuss this question again in § 2.4. 
The displacements at all points of the crystal above A may be 

calculated (Leibfried and Liicke 1949). The singularity which in the 

pure elastic theory arises at the origin is now found at a distance € below A 


outside the upper half crystal, while the displacements in the lower half’ 


crystal have a corresponding singularity at a height ¢ above B. We may 
regard the displacements in A as caused by an elastic dislocation passing 
through the point (0, ja—{), together with a uniform dilatation along 
this line, or alternatively we may regard the displacements at all points 
of the crystal as derived from elastic dislocations spread continuously 
over the x axis with density b¢/m(a*+-@). 

A screw dislocation is not confined to a single glide plane, and the 
assumption of a special law of force between a particular pair of planes 
of atoms does not at first sight seem to provide a satisfactory model. 


eo 
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However, HEshelby (1949 b) has pointed out that if this assumption is 
made, the analysis proceeds exactly as for the edge dislocation, except 
that o in the formulae is replaced by zero. The corresponding value of 
¢ is 3a, so that the stresses on both sides of the assumed glide plane 
coincide with those produced by an elastic dislocation at the origin, and 
are independent of the initial choice of glide plane. 

If the crystal is not infinite in the z direction, but is bounded by the 
planes z=+4)D, the elastic solution for a screw dislocation may be found 
by the method of images. Leibfried and Dietze (1949) have shown that 
if a dislocation lies in the middle of the plate, and we write }D’=1(D—b) 
for the width of the portion of the plate lying above z—4b, eqn. (1) is 
replaced by 
Bf? dv(z’) da’ ‘ 
Py ple Welder auialy [iam \ 0) | goed en. on (5) 
Just as in the case of the infinite crystal, the corresponding Peierls 
equation is satisfied by the displacements produced by an elastic 
dislocation a suitable distance ¢(D/b) below the plane z=4b. The total 
energy is finite. 

For the edge dislocation, the energy of misfit is 

2 (oe) b2 
ae fat cos 27d) dar eon ew Le ON 


The elastic strain energy is | p,,wdz. If this integral is taken between 
infinite limits it diverges. Between limits +R, it yields 
bu Rly tan y dy 
ah, (HETIL ee 
2n*(1—a) Jo Uae 
which for large R is asymptotic to bu In (R/2¢)/47(1—o). 

Lomer (1952) has pointed out that the energy is reduced if the stresses 
are relaxed across the surface of a cylinder of radius R. The stresses 
corresponding to this relaxation are those derived from the stress 
function y! given by (52) of Part II. The reduction in energy is easily 
obtained by integrating along the line 6=—0 as 

Bee bu 
30 (i Pro dr= 877(1—o) ° 


2.2. THE SOLUTION FOR A PAIR OF DISLOCATIONS 


A single straight dislocation cannot be formed within a crystal by 
localized glide motions. In Orowan’s original picture (fig. 4), glide is 
initiated by the formation of a closed loop of dislocation. We take as 
a model of this process the formation of a pair of edge dislocations of 
opposite sign on the same glide plane. This system cannot be in stable 
equilibrium, but if a shear stress is applied to the crystal it is possible 
(Nabarro 1947) to find a position of unstable equilibrium in which the 
attraction between the dislocations balances the force with which the 
_ applied stress tends to separate them. The energy of this unstable state, 
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referred to the homogeneously sheared crystal as zero, is the activation 
energy for the initiation of slip. The configuration is shown in fig. 34. 
It is now convenient to measure displacements in the rows A and B 
from zeros such that ¢=0 represents perfect fit. This changes the sign 
of the sine function in eqns. (2) and (3). Consider the function 
x* sin 6 


b= = cot-t( re —cot#) +7, a fate eed 


where @ is a constant. 
Substituting in eqn. (1), we find that the tangential stress required 
to maintain this displacement is 
by. 


— (sin 27d—sin 4 
Ing (SiB 27p—sin 36), 


Fig. 34 
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A pair of dislocations in a sheared crystal. The angle of shear shown is ten 
times as great as that required to maintain equilibrium in the planes A and B. 


whereas the stress produced by the atoms in the other half of the crystal is 
bu. 
Ing SB 2rd. 


The actual stress exceeds that required to maintain equilibrium by a 
constant quantity (bu/27a) sin 40. The upper and lower halves of the 
crystal can only remain in equilibrium if external stresses of this amount 
are applied to them. This stress would produce a constant relative 
displacement of the planes A and B given by 66/47, which agrees with 
the strain given by eqn. (8) for large values of 2. In practical cases the 
applied stress P is usually much less than », and @ is therefore small. 


_ Mathematical Theory of Stationary Dislocations 355 


For small 6, the two dislocations are widely separated. We may define 
their centres as the points v=+2¢/0, and the distance between them 
w=2x is given by 


w= 2a/(1—o)0, 
while P=bp6/4zra, 
leading to P=by|2r(1—o)w. 


Figure 34 corresponds to 6=n/10. 

These results for small @ agree with those given by elastic theory. 
The stress P exerts on each dislocation a force Pb=b%u6/47a, while the 
attraction between the dislocations is b2u/27(1—c)w, which is also equal 
to b'u6/47a. 

The activation energy to form this pair of dislocations is the increase 
in the elastic energy of the two half crystals, plus the increase in the 
energy of misfit between A and B, minus the work done by the external 
forces. 

If @ is small, the activation energy for unit length of the dislocation 
pair is 6*u In (2/@)/27(1—c). Adding the work bPw done by the external 
forces, we find the energy associated with each dislocation to be 


pie: w 

tn(1—o) (1+ as zi) 
This is the work done by the stresses in the elastic solution when the 
range of integration is cut off at a distance 2¢/e from the centre of each 
dislocation. This cut-off radius agrees with that obtained for a single 
dislocation. For a screw dislocation, the elastic stresses lead to a similar 
expression in which the factor (l—c) is replaced by unity. Dean and 
Mann (1949) have treated the problem of two edge dislocations by the 
methods of elasticity theory. They neglect the boundary condition at 
infinity, and their result therefore depends on the shape of the outer 
boundary. 


2.3. ARRAYS OF DISLOCATIONS 


van der Merwe (1950) has extended the Peierls approximation to three 
problems involving uniformly spaced rows of dislocations. We outline 
his treatment in §§ 2.31, 2.32 and 2.33. Koehler (1941) has calculated the 
energies of lattices of positive and negative dislocations, assuming that 
the elastic solution holds at all points distant more than r, from a 
dislocation, while the energy stored within r, is independent of the 
presence of other dislocations. The lattices he considers contain pairs of 
positive and negative edge dislocations with their joins lying perpendicular 
to the glide plane. The energy is obtained as a lattice sum using 
potentials similar to (65) of Part II. A pair of dislocations is taken as 
reference unit, and the potential of interaction between this and a similar 
unit at a distance r is of order r~?. This is summed over all units lying 
on the axis of the reference unit. The interaction between the reference 
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unit and all units lying on a parallel line at a distance a diminishes 
exponentially with a according to a formula of the type (92), (98) or (99) 
of Part II, and the summation over all lines is rapidly convergent, as in 
the sums of § 5.3 of Part II. The energy W for each pair of dislocations 
is of the form fe 3 
a Ine 

where d is the separation of dislocations in a pair, and y is a constant 
depending on the lattice structure and on the law of force assumed to 
hold in the core of each dislocation. Koehler’s formulae are incorrect 
in detail, as his expression for the interaction energy of a pair of edge 
- dislocations differs from (65) of Part IL by numerical factors. 


- 2.381. Edge Dislocations on the same Glide Plane 

van der Merwe first considered a row of edge dislocations on the same 
glide plane in a simple cubic lattice. This configuration occurs between 
the bent lamellae of fig. 2, or in the interface between a thin sheet of crystal 
and an oriented layer of equal thickness deposited on it. The initial 
stages of the analysis are similar to those of § 2.1. The natural spacings 
of atoms in the layers A and B of fig. 33 are no longer equal, but are 
b, and bg respectively, having a vernier P such that (P+1)b,s=Pbg. 
It is convenient to refer the positions of atoms in each lattice to a 
symmetrical reference lattice C having a spacing bg such that 
(P+1)bs=(P+4)bp=Pbz=p. Then each atom of B corresponds through 
an atom of C with a unique atom of A, and the force between atoms of 
A and B is a periodic function of their relative displacement, with 


period by. The displacement itself has period p, and the analysis is — 


therefore conducted with Fourier series instead of Fourier integrals. The 
governing equation is 


Qn ee 128 ee | [sin 276(X-+t)—sin 276(X—t)] cot Ht dt, 
(10) 
where p=2u'/be 
X = 2ra/p 
and tan d6=7p/(P+4)u(1—<). 


Here / is an elastic constant approximately equal to p, and w’ is the 
displacement measured from the reference lattice. The solution of (10) 
is shown to be 


1 
$(X)=}-+ =tan™ (tan $3 tan £X). oS, Oe ee 


2.32. Screw Dislocations Forming a Grain Boundary 

The resolution into two grids of screw dislocations of the boundary 
between two grains which have been rotated about the normal to a 
common low index plane is illustrated in fig. 24. When the angle of 
rotation +4 about this normal is small, the only important effect of 
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the interaction between the two grids is to reduce to zero the shear at 
points far from the interface, and we shall consider the screws along the 
z axis of fig. 24 (b) separately. Their separation is p=4b cosec 3%, and 
the period of the force of interaction between the adjoining surface is 
b sec $y. Assuming that the force is derived from a nearest neighbour 
interaction, its maximum potential is independent of , and this leads 


to a force law 
_ Bo COS SY . (2:rv cos Fy 
y= 5 sin —— : (12) 
The governing equation reduces to the form (10), with the solution (11). 
From this solution the surface energy 7(%s) of the boundary is calculated, 
the figure first obtained being doubled to allow for the second set of 
dislocations. The result is 


b 
wb) = gos {1+y—(1+y?)¥2—y In [2y {+27}, - (13) 


where y=27pp.. 1 tan 44 sec dy. 

For small % this is of the expected form (cf. (95) of Part II). If we 
attempt to extrapolate (13) to larger values of 4, we find that it does 
not have the necessary property of symmetry about .s=47, which is the 
greatest possible relative rotation of two simple square plane lattices. 
This happens because (12) represents only the interaction between those 
rows of atoms above and below the boundary in fig. 24 which are parallel 
when J=0. When % lies between $7 and #7, it is rows which were 
originally perpendicular which tend to align. van der Merwe proposed 
the interpolation formula 7,()=7(%)+7(47—)—7($7), and the two 
formulae are plotted in fig. 35. 

When 4=3r, all the atoms of one lattice coincide with atoms of the 
other lattice. When tan t¢ is rational, say m/n, some atoms again come 
into coincidence, the atom originally at (nb, mb) falling on that originally 
at (nb, —mb), and this may lead to a reduction in 7(#) for this value of y. 
The following argument (based on a suggestion by Dr. G. A. P. Wyllie) 
shows that the reduction is likely to be small. For any fixed value of %, 
the energy of the surface of misfit may be varied by a rigid translation 
of the two half crystals. During this translation the grids of screw 
dislocations move rigidly across the interface. The energy of each 
dislocation varies as it moves with respect to the atomic lattice, but the 
calculation of § 4.1 shows that the variation is small in a metallic lattice. 
In this case an irrational value of tan 44, which leads to a distribution of 
dislocation energies about their mean, will not lead to a significantly 
greater surface energy than a neighbouring rational value for which a 
suitable translation allows many dislocations to find their positions of 
lowest energy. The surface energy is abnormally low only for values 
of ;s close to those which make ail the atoms on one side of the interface 
correspond with atoms on the other side because we consider only nearest 
neighbour interactions across the interface. If an atom A, on one side 


358 F. R. N. Nabarro on the 


is placed above an atom B, on the other side, the energy is low only if ¢% is 
so chosen that a next neighbouring atom A, lies very nearly above some 
nearest neighbour B, of By. Thus the 7(s) curve for a twist boundary, 
unlike that for a tilt boundary (§ 5.23 of Part IT), should not be markedly 
cusped. Against this, the experiments of Wilman (1951) indicate a marked 
preference for certain rational values of tan $y. 
2.33. Edge Dislocations Forming a Grain Boundary 

The analysis of a wall of edge dislocations forming the grain boundary 
of fig. 21 is similar to that of the preceding sections. The elastic constant 
corresponding to ji) of (12) is shown to be y«/(1—2c), and the width of the 


Fig. 35 
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Energy of a twist boundary (after van der Merwe). A: Elastic energy. 
B: Energy of misfit. C: Total energy. D:Symmetrized formula for 
total energy. E: Total energy according to a formula 7t=7 xh In (yo/x) 
with y% adjusted to give a good fit for small x. 
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dislocation, measured along the line joining two dislocations, is less 
than 6. This means effectively that the solution reduces to Burgers’s 
solution in which Hooke’s law holds everywhere except in the 
neighbourhood of a dislocation, and in this neighbourhood the approxima- 
tion breaks down. The results, however, illustrate an important effect, 
which is shown still more clearly by the calculations for a twist boundary 
(fig. 35). Suppose that in the approximate formula 7(i)=r oi In (es o/b) 
we choose %, in order to obtain the best fit to (13) for small values of ps. 
Then for larger values of ys the energy given by (13) is greater than that 


given by the approximate formula. We shall discuss the interpretation 
of this effect in § 2.42.: 
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2.4, DEVIATIONS FROM A SINUSOIDAL LAW OF FORCE 


If one plane of atoms is laid upon another, its positions of equilibrium 
are separated by lattice vectors 6. It often happens that between two 
neighbouring positions of equilibrium there is another position of stable 
mechanical equilibrium corresponding to a twinned or faulted structure. 
The.curve representing the potential energy of the two layers as a function 
of their relative displacement then has subsidiary minima between the 
principal minima, which are separated by distances b. Even if there are 
no minima, the maxima may lie much lower than those given by the 
sinusoidal law of force (2). Zener (1948 b) has pointed out that the 
common occurrence of phase changes from close packed to body-centred 
cubic structures in metals indicates that the maximum restoring force and 
potential energy are usually much less than (2) would require. 

If the potential maximum is very broad and flat, the energy of misfit V 
between two surfaces is much less than that corresponding to Hooke’s 
law and a sinusoidal law of force, and is almost independent of the 
relative displacement w (§ 2.1) except when w/b is small. If the width 
of the dislocation is J, the surface energy is of order /V, while the elastic 
energy is of order b°u/1. The total energy is thus 


Ged VC aD eee ok to faye ae. «!) s, YGLA) 
which is smallest when Pe OU eee ees Wale, Set LD) 
corresponding to u=2b(byV)1/?. 


The width of the dislocation is proportional to V~1/?. Since the 
curvature near the potential minima is fixed by Hooke’s law, it is easy 
to see that the maximum shear stress is of order p,=(uV/b)'”, and is 
attained at a distance of order u,=(bV/u)' from the position of stable 
equilibrium. We note that lp,,=by and lu,,=b?. 

It is not in general possible to solve in closed form the equation which 
replaces (3) when the sinusoidal law of force is replaced by an arbitrary 
law. 


2.41. The Method of Foreman, Jaswon and Wood 


Foreman, Jaswon and Wood (1951) have given a parametric solution 
for a family of laws of force. We use a notation which corresponds to 
that of § 2.1, except that the distance a of that section is assumed equal 
to b. The parameter a, determines the law of force, which reduces to (2) 
for a)=1, while % is a parameter in terms of which a, wu, and p,, are 
expressed for each law of force. The solution is 


L=A6 tan $y ] 
u=— (b/47a9)[ao+(ao—1) sin #] - 
Dy2=(u/ 4774 9")[2(2a9—1) sin $+ (ap—1) sin 2] 
—7p<n, a>l. 
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For all values of do, the relation between u and p,, agrees with Hooke’s 
law for small strains, but the maximum shear stress p,, is less than that 
given by a sinusoidal law of force. We now define the half width / 
more precisely as the shortest distance in which w changes by +b. Values 
of J and p,, are given for various choices of a) in the following table : 


Ay 1 2 3 4 5 10 ora) 
U/2¢ 1-000 1-255 1-653 2-076 2-507 4-697 0-442a, 
QP m/e 1-000 0-787 0-594 0-471 0-389 0-207 2-202a,-1 
Bubble diameter in mm — 1-9 — 13 — 0-7 — 
rolb 0:56 1:04 149 1:93 2:36 4-53 0-43a, 

iy (degrees) 45 26 17 13 11 aaa) 58a, 1 


The product lp, is almost exactly constant. The fourth row of figures 
in the table gives the diameter of the soap bubbles (Bragg and Nye 1947, 
Bragg and Lomer 1949, Lomer 1949) which give a two dimensional model 
with a law of force approximating to that given by (16) with the given 
value of a). According to Lomer, the law of force for copper is well 
represented by bubbles of diameter 1-2 mm, corresponding to a, rather 
greater than 4. The alkalis are best represented by bubbles about 2 mm 
in diameter, corresponding to a)=2. 

Foreman (unpublished) has proved that the surface energy is 
independent of the law of force in the glide plane, not only for a law of 
the type (16), but for any law. We give here a form of proof due to 
Foreman and Nabarro. Let the potential energy for a displacement wu be 
V(u). Then the total energy in the slip plane is 


~ dV oo dV du 
f= ee V(u) dx=[Va]— [equau=— [#3 ie: <a \tdoy 
The value of dV/du=2p,, is given by (1). Substituting this value in (17), 
we obtain 
a) Qu ee Ree ee Yet! he 
b-—- | =| oo Ta at AD... Se SRE 
Interchanging x and x’, we have 
oye" (Pere ae ae , 
fa | {| = Te qi dvdx’. . . . (19) 


Adding (18) and (19), 


9 (oa) 9) , 
p= |" Gt OM eens 
He {Loot ee eee a meer t 
giving H=b*u/47(1—c) in agreement with (6). 
If we apply a similar treatment in the case of a screw dislocation, the 
solution represents a dislocation extended along the glide plane originally 


selected. It is no longer true, as it was for the elastic theory and for 
the sinusoidal law of force, that the dislocation has circular symmetry 
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and all possible directions of motion are equivalent. If two glide planes 
both containing 6 are crystallographically equivalent, a screw dislocation 
lying in one can only turn into the other by first shrinking until it has 
circular symmetry. This requires an energy of activation. 


2.42. The Structure and Energy of a Tilt Boundary 

According to eqn. (95) of Part II, the energy of a tilt boundary of 
angle % is roughly proportional to % In (%,/ss). This energy rises to a 
maximum when =7,/e, and vanishes when s=,. If %, is not a small 
angle, it is obvious that the approximations have broken down when 
exceeds y,/e. Read and Shockley (1950) have calculated ys, in terms of 
the energy of a single dislocation. If the energy of unit length of a 
dislocation, integrated from the axis out to a radius R (without relaxing 
the tractions across the cylindrical surface r—R) is expressed in the 
form const. In (R/7r,), then 


peetiony ek kA Ue se) Se C20) 


Their numerical application to the case of Peierls’s model seems to 
contain an error: it follows from §§ 2.1 and 2.2 that for this model 
fo=2C/e—b/e(1—oa), giving J)»=e?(1—o)/27=45° for c=}. Lomer (1952) 
has evaluated 7) for the more general law of force used by Foreman, 
Jaswon and Wood. He finds for the energy of unit length 
bu R(1—c) Ayp—1\2 

eS 477(1—c) E Ad ee ( 209 ) ik Bel) 
The corresponding values of r,/b and % are given in the table of § 2.41. 
For copper (a)=4), %y is only 13°, while the maximum energy occurs at 5°. 

The value of r, corresponding to the approximation of van der Merwe 
(§ 2.33) is somewhat different. If we assume his parameter py to be 
equal to », we find r9,=b(1—2c)/2(1—o). 

It is to be expected that for sufficiently large values of ys the energy 
will lie above that given by the simple formula y¥ In (9//), as indicated 
by the calculations of van der Merwe. The approximate formula was 
derived by integrating the energy density round each dislocation from a 
distance of order r, to a distance of order 6/25, where the stress field of 
the next dislocation is of equal importance. If the dislocation is of a . 
width 7 considerably greater than 6, this approximation breaks down 
when 6/2 becomes comparable with }3/, that is to say, % becomes as 
great as b/l. For Peierls’s approximation this is $}=1—o = 40°: for the 
law of force (16) with a,=4 it is J=20°. When + is greater than this 
value, the energy associated with each dislocation is approximately that 
lying in an ellipse with semi-axes b/2 and 31, which may be split into a 
core energy C and a strain energy D/is. The energy of the grain boundary 
in this region would thus be of the form Ct+D. 

The form Os+D is based on the assumption that the length and 
energy of the core of each dislocation remain unaltered even when the 
dislocations are separated from each other by distances less than their 
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width. The experiments of Lomer and Nye (1952) on grain boundaries in 
bubble rafts show that this is not correct: as % increases to a value at 
which the spacing of dislocations is comparable with their length, the 
dislocations become narrower, so that the region of misfit is concentrated 
in the grain boundary. That this is bound to happen can be seen 
theoretically in two ways. Firstly, the width of an isolated dislocation 
is determined by a balance between the strain energy and the energy 
of misfit. In a tilt boundary of large angle, the energy of misfit 
predominates if the dislocations have the same width as when they are 
isolated. The configuration will therefore adjust itself to reduce the 
energy of misfit, the corresponding increase in strain energy being less 
than for an isolated dislocation. Alternatively, we may regard each wide 
dislocation as being a distribution of infinitesimal dislocations along the 
glide plane. If the dislocations form a tilt boundary, the same forces 
which draw the dislocations as a whole into a line tend to draw the 
infinitesimal components of each dislocation into the same line. 

The experiments on a raft of bubbles show that as % increases through 
a fairly small range the form of the grain boundary alters from one in 
which the misfit is concentrated in strips around the glide planes of 
individual dislocations to one in which the regions of misfit lie 
perpendicular to the glide plane. The regions of misfit of individual 
dislocations rapidly overlap, so that it is impossible to distinguish the 
individual dislocations, and the structure reduces to two almost perfect 
grains fitting badly along their boundary. It is probably this region 
which gives a grain boundary energy almost independent of ¥&. 
A theoretical discussion has been given by Shockley (1952). 


2.5, BARRIER PROBLEMS 


The calculation of § 5.1 of Part II gives the stress field surrounding 
a procession of dislocations on the same glide plane when the leading 
dislocation is held up by a barrier. However, it does not give a detailed 
analysis on an atomic scale of the stresses near the leading dislocation, 
nor does it provide any physical model of the barrier. This problem has 
been treated on the Peierls approximation by Jaswon and Foreman (1952). 
They treat first the problem of a single dislocation meeting a barrier. 
Their model is similar to that of Peierls, except that for «>s the law of 
force across the glide plane alters, while the elastic properties of the 
medium on either side of the glide plane remain unaltered. This change 
in the law of force constitutes a barrier at the point x=s to a dislocation 
approaching from the left. If the sinusoidal law of force is replaced on 
the right of the barrier by Hooke’s law, we have a fair model of a grain 
boundary, the slip directions permissible in the grain on the right 
differing from the slip direction of the dislocation in the grain on the left. 
The barrier repels the dislocation, which is held at the origin by an applied 
stress P. The position of the dislocation is defined to be the point of 
maximum misfit, at which w=0. If there is a crack along the glide plane 
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to the right of z=s, the shear stress vanishes, P is negative and the 
equilibrium of the dislocation at «=0 is unstable. 
Equation (1) is now replaced by 


| er ee dul’) dx’ en 
1 Dae) am(1—a) ibe da’ uv’ —o9 +P, . . . . (22) 
while (2) is replaced by 
be 4rru 
| Sy es L<S 
Poa —2uu/b xs (grain boundary) a3) 
{ 0 x>s (crack). 


If the solution u(x) which satisfies (22) and (23) is known, the shear 
stress P may be found from it by a simple virtual work argument. 
We consider a solution representing a single dislocation, so that 
u(co)—u(— co)=—$b. Then in a virtual displacement in which the whole 
pattern of atoms in the planes z=-+4b moves forward in the x direction 
rigidly by a distance 6, the applied stress P does work bP$. The strain 
energy of the elastic continua above and below these planes is unaltered, 
while the relative displacement of a pair of atoms opposite to one another 
in the glide plane increases by —26 du/dz. The energy of the interatomic 
forces in the glide plane thus increases by 

a du 
— 26 ie ‘ D(C) ae dx. 


By the principle of virtual work, this is equal to the work done by the 
external stresses, so that 


ayia du 
ae aes Oe 2) 
P= 5 |, Pes(®) yu 5 RE eee ee © 9) 
This may be written in the form 
2 cur) 
P=—; PolU) du, (25) 
bJ f— 0) 


where the integral is taken over the branch of the function p,.,(w) 
determined by the corresponding value of x. Thus we cannot evaluate 
the integral (25) exactly until w(a2) is known. 

However, if P is small, and the law of force on the right of the barrier 
does not deviate greatly from the sinusoidal law for small displacements 
from the position of stable equilibrium, we may obtain a useful approxima- 
tion. Consider, for example, the case of Hooke’s law, illustrated in 
fig. 36. The value of P is 2b" times the positive area between the wu axis 
and the curve ABCDOFG, shown shaded. Since, for the sinusoidal law, 
the area under the curve LKDOFG vanishes, the area required is that of 
the curve ABCDKLA. We note that AB is of order P, while AL is of 
order P?, The areas ABKLA and HBKH are thus of order P*, Since. 


P.M. SUPPL.—JULY 1952 2B 
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the required area is of order P, we may approximate to it by the area 
HBCDKH. Calling the difference between the two laws of force (23) 
Ap(u), we may write this result 


us) 
pp=2| Apu) duc’... 
—zb 
For Hooke’s law, 
Ap(u)==167?u(u-+ $b)3/30%. i 2= ee 
Fig. 36 
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The sinusoidal law of foree HKDOFG approximates to Hooke’s law NHC or 
QFR near the positions of stable equilibrium H and F. The force on the 
dislocation is proportional to the area ABCDKLA, which for small applied 
stresses OP is nearly equal to HBCDKH. 


Tras small, s is large, and we may assume that u(s) is given closely by 
the solution (4) which holds in the absence of a barrier. This gives 


U(s)+- 3b = b¢/27s = b*/47(1—o)s. vigts ) weed 
Substituting (27) and (28) in (26), we obtain 
P=pb4/907°(I—o)*es. 5 5 eee atte oa) 


As a typical numerical value, s=5b when Paxil: 
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We may now extend the result to find an approximate solution for a 
procession of n dislocations held up by a barrier. Our treatment will be 
valid provided that the distance from the leading dislocation to the 
barrier 4s appreciably less than its distance from the second dislocation. 
According to (29), the repulsion of a barrier falls off as s~4, and our 
assumption about the spacing of dislocations means that the forces 
exerted by the barrier on all dislocations other than the first are 
negligible. In this case, s is given by (29) with P replaced by nP. 
If the distance between the first and second dislocations is d, the 
total stress on the second dislocation, which must vanish in equilibrium, 
Is (n—1)P—by/2r(1—c)d. Thus d=by/27(1—c)(n—1)P, and the 
condition d>s which must be satisfied if this simple treatment shall be 
valid may be written 


nP <0-85p. a es ae ee Ei), 


We have here neglected the difference between n and n—1. 

If, for the displacements actually occurring in this region, the law of 
force beyond the barrier does not differ greatly from the sinusoidal law, 
an approximate solution of (22) and (23) may be found by a perturbation 
method. We take as zero order solution uw, the solution (4), and write 
* U=Uy+U,, Where uw, and its derivative are small. This leads to a linear 
integral equation for u,. Jaswon and Foreman have obtained the 
general solution of this integral equation in series form, using the sets of 
functions 


C,,(*)—(—1)"=cos (2n tan~1 x)—(— 1)” 


and S,(v)=sin (2n tan“! x). 


§ 3. THz APPROXIMATION OF FRENKEL AND KONTOROVA 


Frenkel and Kontorova considered a mechanism of slip differing from 
that of dislocations. Their basic slip process is the propagation of a 
compressional wave along a line of atoms. The atoms in the front portion 
of the line are still in their original positions, while the atoms of the rear 
portion have moved forward by one lattice distance 6. Formally, this 
may be considered as the motion along the line of atoms of a small closed 
loop of dislocation with its Burgers vector b lying along the line, and 
the motion of a straight dislocation line may be considered as the 
resultant of many of Frenkel’s processes. The motion of a Frenkel 
‘ caterpillar ’ does not relieve a shear stress applied to the crystal, and 
caterpillars cannot be set in motion by uniform external stresses; the 
motion is of the kind described in § 4.1 of Part II as ‘ ineffective’. The 
physical process considered is probably not of importance in plasticity, 
but the mathematical analysis provides a useful simplified model of a 
dislocation (Dehlinger and Kochendorfer 1940). Paneth (1950) has 
suggested that this mechanism may sometimes be important in diffusion. 


2B2 
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All atoms ‘other than those in the moving chain are considered fixed. 
Each atom in the moving chain has an energy of interaction with the 
substrate which is a sinusoidal function of its position x, and there are 
elastic forces between atoms of the chain proportional to the deviation 
of the distance between nearest neighbours from its equilibrium value a. 
The potential energy of the chain is thus of the form 


AX(1—cos 27a;,/a)+402(%,41—2%,—4)*. » ike ele 
k k 
In the state of lowest energy, each atom lies at the bottom of a potential 
trough, and we shall assume that the chain is very long and that 
deviations from this lowest state occur only near the middle of the chain. 
The summation in (31) may then be taken from k=—o to k= oo. 
Write , 
y= (k+ €;)a, 

and (31) becomes 

U=A2(1—cos 27 €,)+4aa72(E,,,—E,)?»  - « « (32) 

k k 


The condition of equilibrium for kth atom is dU/d&;,=0, or 


(S44 — 26+ Sp-1)= 277A sin 27; - + + (33) 
The general solution of this difference equation is not known, although 
Frank and van der Merwe (1949 a) have illustrated some simple special 
solutions. In the special cases for which the displacements of neighbouring 
atoms are nearly equal (€,—&,_,<1), the difference equation may be 
replaced by the differential equation 


att® dEldk*==2arA. gin Dore 4. ee 
which integrates to give ~ 

dé\2 2A 

Tk ram Ea 008 2). 6. ee ea 


The constant of integration is determined by the fact that for a single 
dislocation € and dé/dk both vanish in the unstrained parts of the chain. 
The solutions obtained by other choices of the constant of integration 
represent either stable infinite sequences of equally spaced dislocations all 
of the same sign, or unstable sequences of dislocations of alternating sign. 
They are expressible in elliptic functions. 

Integrating again, we find 


tan $7g=exp [ta(k—x)/lb], . . . . . . . (36) 
where Ly?=aa/4A. RRR A i ns 


The signs + correspond to dislocations of opposite sign, while k=« 
defines the position of the middle of the dislocation. The solution may 
be written in the more symmetrical form 

tan [Jr(¢—4)]=tanh [$a(k—x)|lo]. 


In this one dimensional model the displacement produced by a single 
dislocation vanishes exponentially at large distances. 
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Frank and van der Merwe, and also Kochendorfer and Seeger (1950), 
have extended this treatment to the case in which the natural spacing 
in the chain 6 differs from the period of the substrate a. Equation (34) 
and its solutions are unaltered, but the expression for the tension in the 
chain is altered. The behaviour of a chain with free ends is shown to 
depend on the parameter of misfit m=(b—a)l,/a. If the average lattice 
spacing in the dislocated chain is 6, we define m by m—=(b—a)l,/a. 
A function A is calculated such that in the state of lowest energy 
m=—=A(m). A second function B(m) is found such that if m<B(m) more 
dislocations form spontaneously without activation energy, and a third 
function C(m) such that if m>C(m) dislocations spontaneously leave the 
chain. If A(m)>m>B(m), the energy is reduced if the number of 
dislocations in the chain increases, but this increase requires an activation 
energy. The functions A, B and C are shown in fig. 37. 


Fig. 37 


Density of dislocations in a monolayer m as a function of the misfit m between 
the monolayer and the substrate (after Frank and van der Merwe). 
A: State of lowest energy. B: Below this curve, dislocations form 
spontaneously in the monolayer. C: Above this curve, dislocations 
spontaneously leave the monolayer. 


The theory has important applications to oriented overgrowths on 
crystal faces (van der Merwe and Frank 1949, van der Merwe 1949, 
Frank and van der Merwe 1949 b). However, Smollett and Blackman 
(1951) have objected to some of the applications of the theory, showing 
that a monatomic layer strained to the extent which Frank and van der 
Merwe postulated would in some cases be mechanically unstable. They 
also disagree with Frank and van der Merwe on the interpretation of 
some of the experimental data. Further points have been raised by 
Frank (1951 c), Blackman (1951) and Lucas (1951). 
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3.1. INTRODUCTION OF A SECOND HARMONIC 
Frank and van der Merwe (1949c) have given an explicit solution of 
the differential equation for the model of Frenkel and Kontorova when 
the potential A cos 27x/a of (31) is replaced by the more general form 
A (cos 2ma/a+r? cos 47a/a). If r?<}, the effect of the added term is to 
sharpen the minima and broaden the maxima of the potential, thus 
providing a better approximation to the potential in a real crystal. The 
limiting case of r?= } is shown in fig. 38 (a). Provided that r<i, the 
total range of the potential is still 24; if A is determined from the 
curvature at the potential minimum, its value is 1/(1+4r°) of that 
corresponding to a simple sinusoidal approximation. The introduction 
of a second harmonic term with 7?<} has little effect on the properties 


of the model. 
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(4) (0) 
Effect of a second harmonic term in the potential between two planes of atoms 
(after Frank and van der Merwe). (a) Amplitude of second harmonic 
a quarter that of fundamental. (b) Amplitude of second harmonic four 
times that of fundamental. 


If 7?>4, the potential curve has a subsidiary minimum, as shown 
in fig. 38 (b), which corresponds to 7?=4. The properties of the model 
are now more complicated, for in addition to the dislocation in which 
n particles occupy a space of approximately (n--1)a, there may be partial 
dislocations in which the particles to the left of a region of misfit lie 
near to deep minima, while those on the right lie near to subsidiary 
minima. ‘The model approximately represents the partial dislocations 
of § 5.21, although these are not plane deformations and involve shears 
which the model neglects. The dislocation lies between two close packed 
planes ; the deep minima are sites in the crystal lattice, while the shallow 
minima are sites in the twinned or faulted lattice. 


3.2, PERTURBATION THEORY 


Seeger and Kochendorfer (1951) have developed a perturbation theory 
for the model of Frenkel and Kontorova. They apply it to a dynamical 
problem, and also to the static problem of an array of two or more 
dislocations not in equilibrium. As they point out, no true equilibrium 
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solution exists for the problems they examine; their solutions found by 
perturbation theory have the property that they fail to satisfy the 
governing equation less grossly than do the zero order solutions. These 
zero order solutions are obtained by joining together portions of different 
exact solutions. In the cases studied in detail (pairs of dislocations of 
the same or opposite signs, three dislocations of the same sign), the 
effect of the perturbation is mainly to increase the distance between the 
dislocations. It is not clear whether any significance can be attached to 
the perturbed solutions other than the obvious one that when the 
dislocations are farther apart they interact less, and therefore the 
condition of equilibrium can be more nearly satisfied. 


§ 4. [INTERACTION OF AN EDGE DisLocaTION WITH THE Larrick PoTENTIAL 


In treating a straight edge dislocation in a simple cubic lattice, we 
have hitherto assumed that the dislocation lies symmetrically in the 
lattice, so that the added half plane of atoms lies midway between two 
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Atomic displacements as a dislocation moves through the lattice by a 
distance $b. (a) Symmetrical position (~=0) in which the added half 
plane ends at the centre of the dislocation. (b) Unsymmetrical position 
(w=4). (c) Symmetrical position (a=}) in which two rows of atoms above — 
the glide plane match one row below, this row passing through the centre 
of the dislocation. 


half planes of atoms in the other half crystal. As the dislocation moves 
through the lattice, it passes through unsymmetrical configurations to a 
different symmetrical configuration in which one half plane of atoms 
on the expanded side of the glide plane lies midway between two half 
planes on the compressed side (fig. 39). Further motion passes through 
unsymmetrical configurations back to a state equivalent to the original, 
During this displacement the energy of the dislocation varies, and. the 
dislocation can only move through the lattice if a finite force acts on it. 
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In the next section we calculate the variation by a method due to 
Peierls (1940). The positions of minimum energy for a dislocation in a 
lattice are straight lines. If the shape of the dislocation departs slightly 
from a straight line, its energy increases. We estimate this increase 
in § 4.2. 
4.1, CHANGE IN ENERGY OF A STRAIGHT EDGE DISLOCATION AS IT MOVES 
THROUGH THE LATTICE 
The treatment is an extension of that used in § 2.1. We assume that 
the displacement of any atom in the interface is given by (4) of that 
section. Instead of calculating the energy of the interface by the 
integral (6), we take the corresponding sum over the positions of the 
centres of the atoms above and below the interface. The result now 
depends on the displacement «b of the centre of the dislocation from the 
position shown in fig. 39 (a). We shall not make the corresponding 
substitution of a sum for the integral which gives the elastic strain 
energy. The inconsistency of using displacements calculated for a 
continuum and then using a summation prevents a clear decision whether 
this term should be included or not; its inclusion alters only the 
numerical coefficient in (40). The energy of misfit is 
3 ioe) 
v Ata S [1+ cos 2{tan-! 2(1—o)(a+4n)b/a}]. . . (38) 
oO 


== 
| 81rd 


This may be transformed by using the relation 
Pi dea oe 5, f(z) cos 2ras dz . ... . Saige) 


into a rapidly convergent series, as was done by Peierls and by Nabarro 
(1947). However, it may be summed exactly, for we have 


47?(1—o)?v e ; 
eT 2 . {[a/b(1—o) ?+[n+ 2a}? } 1. 
This sum is readily evaluated by considering the poles of the function 
{[a/b(1—o) }?+[z+ 2a]?! cot zz, and we find 

bu tanh (47¢/b) 


ne 4ir(1—o) 1—cos 47a sech (4rZ/b)" © © * * (40) 


Since exp (47¢/b)>1, this may be written, to a good approximation, 
bu 
~~ 4(1—o) 
On account of the approximate method of calculation employed, the two 
symmetrical configurations of the dislocation (z=0 and «=4) have the 
same energy. The argument of the exponential is 27 times the number of 
atoms in the half width of a dislocation. Ina metal, the deviation of the 
law of force between two planes from the sinusoidal form will in general 
increase the half width of the dislocation (§ 2.4), and the amplitude of v 
is reduced. Foreman, Jaswon and Wood have obtained the formula 


v [1-+2 cos 4zra exp (—47¢/b)]. . . . (41) 
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corresponding to (40) for the family of laws of force (16). Their result 
confirms the exponential decrease in the amplitude of » with increasing 
width of the dislocation; the amplitude is reduced by the factor 
(1+ K+ 6K") exp (—K), where K=27(a,—1)/(1—o). 

In a lattice such as that of diamond which is held by directed bonds 
the concept of width of a dislocation loses its meaning, and the 
approximations on which (40) is based break down completely. It is 
clear that in such a case the amplitude of v no longer contains this very 
small exponential factor. 

The force acting on unit length of the dislocation is 

ldv 2bu | dort 

Pam 5 = [asin da exp| —=E I. Sau (42) 
If the dislocation is to move continuously through the lattice, there 
must be an external shear stress P producing a force bP exceeding the 
maximum value of #. This leads to a shear strength 

2 Anl 

For a partial dislocation in a close packed lattice, a=b./2, and the 
value of y/P is given in the following table for various values of o: 


a 0-2 0-25 0-3 0-35 0-4 
fP 27,000 52,000 110,000 280,000 810,000 


For the simple cubic lattice with a=b and c=}, p/P is given as a function 
of the parameter a, in the following table : 


Ay 1 2 4 10 
u/P 5000 2x108 5x10% 4x 10% 


Although the simplest assumptions lead to a strength greater than 
that of a pure single crystal, a flattening of the law of force corresponding 
to a)=3 is sufficient to make this force negligibly small. Moreover, as 
Shockley (1947) has pointed out, arguments such as that of the following 
section show that a dislocation will not usually lie exactly along a line 
of minimum energy, but will ride across the crests and troughs of the 
potential v. Then, by arguments similar to those of § 7.3 of Part I, 
we see that the stress required to move a whole dislocation is much less 
than the maximum stress acting on any one element of it. 

Dietze (1952 a) has suggested that to regard the interaction across the 
glide plane as acting between lattice points is as crude an approximation 
as to regard the interaction as acting between continuous surfaces. 
It is more realistic to consider the centres of force as spreading for some 
distance around each lattice point; the value of v then lies between 
the value zero given by the approximation of a continuum and the 
value (40) given by the approximation of lattice points. Dietze considers 
in detail several cases in which the spread of a centre of force is given 
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by a Gaussian function s~! exp (—2*/s?), where s is a parameter defining 
the extent of the spread. For s=0 the model reduces to one of lattice 
points, and for s= o toa continuum. A reasonable value of s is probably 
about 4b. The amplitude of v is reduced by a factor exp [—(zs/b9)"], 
where 6, is a constant which depends on the nature and orientation of the 
dislocation. Thermal vibrations are regarded (Dietze 1952 b) as increas- 
ing the value of s. 

Dietze solves the problem by considering the case of a continuum as a 
zero order approximation, the difference between the integral equation 
for the real problem and that for the continuum being regarded as a 
small perturbation, as is the applied stress P. The condition that the 
perturbed equation should have a solution is that the perturbation should 
be orthogonal to the zero order solution, and this condition is expressed 
as a relation between P and «a. 


4,2. THE ZIG-ZAG DISLOCATION 


In considering a dislocation which approximates to a given curve, we 
have first to decide whether its equilibrium form is a smooth curve or a 
series of steps lying parallel to lattice directions and approximating to 
the smooth curve. We have then to decide by how much the energy of 
the equilibrium form exceeds that of a straight dislocation. 


Fig. 40 


<---> 


A single tooth of a zig-zag dislocation passing through the points X, Z and Y. 
The full line of slope 5 represents the equilibrium form: the broken lines 
XTUZ and ZY are possible extreme forms. 


We consider only the case in which the tangent to the given curve is 
inclined at a small angle to the axis of a free straight dislocation (Mott 
and Nabarro 1948). 

Figure 40 represents one complete tooth of a zig-zag dislocation which 
is constrained to lie along one equilibrium line in the lattice at a series 
of points such as X and Y, separated by a distance A much greater 
than 6, and to lie along the next equilibrium line at a distance 4b at a 
series of equidistant points such as Z. We allow it to take any form 
intermediate between that of the dotted line XTUZ and that of the 
dashed line ZY. We vary the angle % at which the dislocation line crosses 
from one equilibrium line to the other, and minimize the energy. In 
practice we carry out the calculation only for small values of %, and 
observe whether the energy increases or decreases as ys increases from its 
least possible value, b/A. 
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The energy of the zig-zag dislocation exceeds that of a straight 
dislocation lying along XY because the energy of misfit in the glide plane 
is increased, and because the elastic energy of the blocks above and 
below the glide plane is increased. As % increases, the region in which 
the misfit is unnecessarily bad becomes smaller, and the energy of misfit 
diminishes, but the same elastic displacements are concentrated into a 
smaller space, and the elastic energy increases. 

We first estimate the change in the energy of misfit. We assume that 
an element of length dy at a distance x from the line XY has the same 
energy of misfit as an element of the same length in a dislocation lying 
parallel to the line XY (the y axis in fig. 33), and passing through the 
centre of the element. This energy exceeds that of an equal element 
lying in the line XY by an amount given by (40), and equal to 


[2a(1—o)]-1b’u exp (—42¢/b)(1— cos 472/b) dy. 


This expression vanishes over the parts of the zig-zag parallel to XY. 
The sloping parts of the zig-zag occupy a fraction b/Ay of its length, and 
the average misfit energy per unit length of the zig-zag dislocation 
exceeds that of a straight dislocation lying along XY by 


W j= [27 Ap(1—o)]-1b63u exp (—47Z/b). . . 1. (44) 


Next we estimate the elastic energy of the zig-zag dislocation. For a 
straight dislocation along the y axis the displacements uw and w are 
given by (8) of Part Il. There is no v component of displacement along 
the y axis, and w and w are independent of y. We assume that the 
zig-zag dislocation also produces no v component of displacement, but, 
near the sloping parts of the zig-zag, wu and w now depend on y. The 
shear strains between neighbouring planes perpendicular to Oy are given by 


Cay = sy OU/OX, Coy=ay Owjdx, . . . . . (48) 


The strain energy per unit length near the sloping parts exceeds that 
of a straight dislocation by 


Qu { i (yy? Lee?) dx de= dpul? | | | (Gs) ae () | Ao Oy 


(This corrects a numerical error in the original paper.) If this 
integral were extended to infinite values of a and z it would diverge 
logarithmically. It is, however, clear that our approximation of 
considering that the displacement of atoms in any plane perpendicular 
to Oy depends only on the point at which the dislocation meets that 
plane, and not at all on the displacements of atoms in neighbouring 
planes, must break down in regions of the plane distant 4 or more from 
the dislocation. The variations from plane to plane become smaller at 
these distances, and the shear strain energy (46) is reduced. An estimate 
of the rate at which the shear energy falls off with distance is given in the 
Appendix to the original paper. The result is to multiply the integrand 
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in (46) by a factor exp [—47(2?+-2*) )7/A]. This new integral is taken 
over a fraction b/A; of the total length of the dislocation, and the average 
elastic energy per unit length of the zig-zag dislocation exceeds that of 
the straight dislocation by 

W ,=(ubp/A)I(A/b), 


where 


1(4/b)= | | | (ss) a (5) | exp | - wo] dedz. . .* (47) 


The values of du/dx and dw/dx must be inserted from § 2.1 of Part IT. 
The approximate value of the integral I(A/b) is derived in the paper 
quoted, and we finally obtain 
kbeup A : 

W=TeaMGe °° et (48) 
where (3—20)?>k>(1—2c)?, and €=a/2(1—oc). Here the numerical 
factor 1.3 used in the original derivation has been replaced by the more 
accurate value $7. We shall take k=1 as an estimate of the true value. 

The total energy per unit length is 


We Ween 


exp (—47¢/b) 1 In (A/4x€) 
27 A| hana es +]. Se 
Varying %, we obtain a minimum of W ae 
8 exp (—47¢/b) 
as a oes I—o) In (Jan) (50) 


To obtain a numerical value, we take o=4, a= /2b, A=100b, giving 
p=2-44 exp (—37/4/2)=3-1x 10-3. 

This value of % is small, so that our original assumption is self-consistent. 

However, the minimum value of % allowed by the geometry of fig. 40 

is #=b/A=10~, and this value, which is also small, is the angle giving 

the lowest possible energy. If A is rather larger, % takes the value 

given by (50). The energy is then clits en Ae 


In (A/4z 
wae [4nt) exh aan, 2 ee 


8(1—c) 
which, with the numerical values taken above, is 
W=0-61(b%u/7A) exp (—27¢/b). . . . . . (52) 


A straight dislocation of length A has become a zig-zag of length 
A+b (cosec 4—cot )==A+4byp, and its energy has increased by AW. 
The tension # in it is, therefore, 24W/by, which is‘ 


sal aad’ 4 | 
H= ry Vee Pee tte Ko mace onc (53) 


With the same numerical values, this gives 
B= 0:16.08 us 5 eo ee ee 
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We conclude that a dislocation in a strained lattice will approximate 
to a smooth curve rather than to a square zig-zag, and that its effective 
tension when it deviates little from a straight line is approximately 
0-16 6%. This value may be compared with that given by (119) of Part II 
for a circular dislocation, in which a logarithmic term also appears, but 
the numerical factor 1/47 is replaced by 5/167. 

If % is given by =6/A instead of by (50), the second term in (49) is 
larger than the first, and the energy per unit length is 


— Ob4y uA 

167/12 = Te s 

where @ lies between 1 and 2. The increased length of the dislocation 
is now 6?/2/, and # is given by 

6b'44, A 
B= Bn ink? 
which lies between the value given by (133) and half of that value. The 
general conclusions are unaltered. 

We may obtain (48) by an alternative method. Two segments of the 
dislocation were originally of total length b/s and of pure edge type. 
The original energy per unit length was approximately 

bu A 
4ir(1—o) a 4nl’ 


W 


so that the original energy of these two segments was 
b* A 
aS 
A4n(l—oa)b  4rrF 
The segments are now inclined at an angle js to their original direction, 


while their Burgers vector is unaltered. They are now partly of screw 
type, and their energy per unit length is reduced to 


(od) ($28 ry) = (a) oh 


The length has increased from b/s to b sec #/, and the total energy is now 


ONT aoe 2)(1 Ay? 
(aay ee) (1—oys")(1+- 39"). 
Since o is necessarily less than 4, the increased length outweighs the 
reduction in the energy of unit length, and we obtain 

(1—20)b yp A 
Mos 8ir(1—o)A Marl’ 


which agrees with (48) if we take k=2(1—2c)/(1—o). For o=4, this 
coincides with the value k=1 already assumed. 

The value H=0-16 6% obtained in (54) is smaller by a factor of 
about 5 than the estimate of the line energy of a dislocation obtained by 
other methods (Frank 1950a). This result arises partly from the small 
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value of A, here assumed to be 100, a value probably appropriate to a 
precipitation-hardened alloy. For a soft material, the value A=1 0006 
would be more reasonable, which increases the value of H to 0-34 b'p. 
The rest of the discrepancy corresponds to the appearance of a numerical 
factor 1/47 in the logarithm in (53), whereas the corresponding factor 
in (119) of Part II is 2/e. Clearly the approximations used in deriving 
(54) could introduce errors greater than this. However, the difference 
‘may be real, for there is no obvious reason why the energy per unit 
length required to bend a straight dislocation into a sinusoidal curve of 
small amplitude and wavelength A should be precisely equal to the total 
energy per unit length of a circular dislocation of radius A, or to that of a 
straight dislocation lying along the axis of a cylinder of radius A. 


§ 5. DisLocaTIoNs IN SPECIAL CRYSTAL STRUCTURES 


So far, we have dealt with those properties of dislocations which are 
common to all crystals, both those which are lattices and those which 
have more than one atom in the unit cell. We now turn to those 
properties which distinguish one crystal structure from another, and 
examine such problems as ‘ What determines the glide system in any 
particular crystal? Why do metals with cubic close packing work 
harden more rapidly than those with hexagonal close packing? Why 
does the body-centred cube give deformation twins, whereas the 
face-centred cube does not?’. The answers to these questions usually 
depend on the idea of partial dislocations. Other new ideas, in particular 
those of dislocation reactions and dislocation nodes, can be illustrated 
from the properties of complete dislocations. Throughout the discussion 
we shall draw freely on the papers of Frank (1951 a) and Read and 
Shockley (1952). 


5.1. REACTIONS AND NODES 
Two dislocations with Burgers vectors b, and b, may combine to 
form a single dislocation with Burgers vector b,=b,+b,. This is a 
dislocation reaction. We may illustrate it by considering the dislocations 
[3a, 3a, 0] and [—}a, 0, $a], which both lie in the (111) glide plane of a 
face-centred cubic lattice of parameter a. Suppose these dislocations 
both lie along lines parallel to [211]. Then (cf. § 4.63 of Part IT) they 
attract one another, and coalesce. By considering the residue of 
displacement for a circuit surrounding both dislocations, which is 
unaltered by the motion and coalescence of the two dislocations, we see 
that the resulting dislocation has Burgers vector [0, $a, 4a]. This 
reaction is represented by the equation 
3@[110]+ 4a[101]=4a[011). ee 
This dislocation, which was made in the plane (111), also lies in the close 


packed plane (111). It can dissociate in this plane according to the 
scheme 


2a[011]=$a[110]+ 4a[ 101}. She sei 


—— = 
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The reactions (55) and (56) differ in that the former proceeds with release 
of energy, whereas the latter absorbs energy. If we neglect the factor 
(1—o) which relates the formulae for the contributions of the edge and 
screw components to the total energy of a dislocation, we have a simple 
rule for deciding whether a reaction of the type b, +b,—>b, releases energy. 
The initial energy is proportional to b,2-+0,2, the final energy to 
(b,+b,)’=6,°+b,+ 2b,b,. Thus energy is released if b,b, is negative, 
that is, if the angle between b, and b, is greater than a right angle. The 
neglect of the factor (l—c) is usually justified because we are already 
neglecting the influence of the crystal anisotropy on the energy of a 
dislocation. 

In general, a dislocation which is not the shortest possible in a lattice 
will dissociate. Thus, in the face-centred cubic lattice, the double 
dislocation a[110] dissociates into two dislocations 4a[110]. The disloca- 
tion reaction a[100}+4a[110]+4a[110] involves no change of energy in 
the present approximation. However, the dislocations such as 4a[110] 
can reduce the total energy by extending in the close packed planes 
(§ 5.2), while the dislocation a[100] cannot. In the body-centred cubic 
lattice, the reaction a[100}~a[111]+4a[111] increases the energy, since 
the angle between [111] and [111] is obtuse. The dislocation a[100], 
though not the shortest possible, is stable. Such dislocations may be 
removed by reactions of the type 


@lO0l-alOlOlstalLi1}--deflll].” . . ... (87) 


involving four dislocations. 

Two dislocations which lie in the same plane may combine to form 
a third even though their lines are not parallel, provided their lines and 
Burgers vectors all lie on the same plane. The dislocations b, b, b, then 
meet to form a node. Frank has suggested that it is then convenient 
to use a new sign convention, the Burgers circuit for each dislocation being 
taken in such a direction that it appears to be right handed when looking 
outwards from the node. The three dislocations are treated equally, 
and their Burgers vectors satisfy the symmetrical relation b,+b,+b,=0. 
More than three dislocations may meet at a node. If we remove the 
restriction that a reaction such as b,+b,—~—b, can occur by a gliding 
of the dislocations which moves the node, we no longer require that the 
three lines and the three Burgers vectors should lie in one plane. Thus 
the four dislocations in (57) can form a node, but the reaction (57) cannot 
occur by glide. Three dislocations free to move in the same glide plane 
can only form a stable node if the angles between their Burgers vectors 
are all obtuse, so that the triangle with sides b,, b, and b, is acute angled. 
This is because the equilibrium of the node requires that the tensions 
in the dislocations should form a triangle of forces with sides proportional 
to 6,2, 6.2 and 6,2. Hence b,?<6,’+6,", and b, b;<0. This is just 
the condition that the energy should fall when b, and b,; combine to 
give —b,. Similarly b, b;<0 and b, b,<0. 
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5,2, PARTIAL DISLOCATIONS 

If a crystal forming a lattice is cut along a half plane lying between 
two sheets of atoms, and the cut surfaces are displaced relative to each 
other by a lattice vector lying in the plane of the cut, we say that the 
crystal has slipped on this plane, and the boundary between the slipped 
and the unslipped regions is called a complete dislocation. If the Burgers 
vector of the relative displacement does not lie in the plane of the cut, 
one or more layers of atoms must be inserted between the cut surfaces or 
removed from them. If the displacement is still a lattice vector, the 
crystal is still perfect except near the edge of the region inserted or 
removed, which is again a complete dislocation. If the crystal structure 
has a basis, the same process may be carried out: the regions added or 
removed must be made of complete cells. 

It often happens in real crystals that a relative displacement of the 
two cut surfaces by an amount which is not a lattice vector, possibly 
combined with the addition or removal of layers of atoms, produces a 
mechanically stable configuration. 

The configurations assumed to be of low energy are usually those 
which are observed to occur in the composition planes of twinned 
crystals. The region of the cut is now a fault plane, and its edge is a 
partial dislocation. The properties of partial dislocations with their 
Burgers vectors in the fault plane are entirely different from those of 
other partial dislocations ; we discuss the two types in § 5.21 and § 5.22. 
The formal definition of the Burgers vector of a partial dislocation presents 
certain difficulties, particularly in the case in which a layer of atoms has 
been inserted which does not belong to the structure of the original 
crystal above or below the cut. It seems necessary for the Burgers circuit 
to begin and finish on the cut, in a region free from large elastic strains. 


5.21. Mobile Partial Dislocations 

If the Burgers vector of a partial dislocation lies in the fault plane, 
it can glide in this plane just as a complete dislocation glides in the 
slip plane. In the absence of any applied stress, there is a force on 
unit length of the partial dislocation tending to move it in the direction 
which will reduce the area of the fault. This force is equal to the 
surface energy y of the fault plane. For a stacking fault in copper, 
y=18 dyne cm (Fisher and Dunn 1952). If the partial dislocation 
moves out of the fault plane, it leaves behind it a region of misfit which 
in general has a very high energy, probably at least as large as the 
energy of a large angle grain boundary, which in copper is 535 dyne em-}, 

The most important partial dislocations of this kind are those 
occurring in the close packed metals. The cubic and hexagonal close 
packed structures both contain close packed sheets of atoms, and the 
relation between one sheet and the next is the same in both structures. 
They can only be distinguished by considering at least three planes. 
In the cubic lattice these lie (in the usual notation) A B GC, in the 
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hexagonal, A B A. If only forces between nearest neighbours were 
important, there would occur not only the regular structures 
ABCABCABCABC...and...ABABABABABAB...and 
the cubic twin... ABCABCBACBAGC..., but also the faulted 
hexagonal structure... ABABABCBCBCB... and such irregular 
structures as ... ABCBCBCBCABC..., which contain layers 
of cubic and hexagonal packing. Any arrangement in which the same 
letter does not occur twice consecutively is permissible. In practice, 
many metals change from cubic to hexagonal close packing with change 
of temperature, while cobalt in particular shows irregular piling. It 
therefore seems that the crystal structure of metals is determined largely 
by nearest. neighbour repulsions, and that the model of a pile of hard 
spheres provides a useful analogy. It is at once evident on handling 
such a model (Millington and Thompson 1924, Mathewson 1944) that the 
easiest movement is not in the glide direction [110] but in the twinning 
direction [121], where the plane is taken to be a cubic (111) plane. This 
leads from the cubic lattice... ABCABCABCABC... to the 


arrangement ...ABCABCBCABCA..., or from the hexagonal 
structure ...ABABABABABAB... to the arrangement . 
ABABABCACACA.... If this movement occurs over part of a 


plane, the boundary of the slipped region is an imperfect dislocation with 
Burgers vector ja[121]. Partial dislocations of this kind were first 
discussed by Heidenreich and Shockley (1948). Two such partial 
dislocations may combine to give a complete dislocation according to 


pee eehene PHN etal2ilj—tal lol, . . 0. % (88) 


The resolution of a dislocation 4a[110] into two dislocations 4a[121] 
and 4a[2I1], all lying in the (111) plane, is illustrated in Plate IX, 
where the plane of the figure is taken as (111) and the horizontal direction 
as [110]. The two partial dislocations are readily seen by looking along 
[101] and [011] directions; the existence of a dislocation $a[110] is most 
easily verified by counting steps round the circuit indicated. However, 
there is no exact analogy, for the glide plane in Plate IX is not the 
(111) plane of the figure, but the plane (111). If the extra rows of atoms 
along [101] and [011] are continued as extra planes (111) and (111) 
respectively, the line along which they meet is [110], which does not lie 
in the glide plane. The two partial dislocations of Plate LX cannot 
separate, but since the angle between them is acute, the complete 
dislocation in the real crystal would always dissociate if the irregular 
piling were not a state of higher energy than the regular. When the 
partials are separated by a distance r, the force of repulsion on unit 
length of each is of order b, bou/27r=a*p/24rr. If the energy of unit 
area of the surface of misfit is y, the force on unit length tending to 
bring the partials together again is y. Their equilibrium separation 
is thus given by 

r=a"/24ry. 
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Heidenreich and Shockley suggest that y may be estimated from the 
energy change EF per atom in the transition from hexagonal to cubic 
close packing according to a formula of the type y=8H/,/3a*, which 
gives r=/3a41/1927H, and should be of the right order of magnitude. 
They suggest that H may be about 100 cal mol-'=7 x 10-¥ erg per atom, 
while a3:=+5 x 10-1! erg per atom, giving a separation of r=20a. This 
suggests that cases may occur in practice for which the two partial 
dislocations are tightly bound together, and other cases in which their 
binding is so weak that quite small stresses will separate them. The 
approximate treatment of Frank and van der Merwe (1949) (§ 3.1) 
may be applied to this problem. Leibfried and Dietze (1951) have 
attempted to find solutions using an extension of Peierls’s method, 
allowing for the two components of the displacement in the glide plane. 
They find approximate solutions by a variational method, minimizing 
the energy. Results are given for the (100) plane, and also for half 
dislocations lying along various lines in the (111) plane. A solution is 
also given for a pair of partial dislocations which combine to form a 
complete dislocation in the (111) plane. Since their approximation 
neglects the energy y of a stacking fault, there should be no equilibrium 
until the partial dislocations have repelled each other to infinity. In fact, 
they find a solution in which the partial dislocations are separated by 
only a few interatomic spacings. The erroneous result arises because the 
form of solution they assume, while adequate for very small separations 
of the partial dislocations, cannot represent two partial dislocations 
widely separated; for wide separations it represents two dislocations 
with Burgers vectors nearly antiparallel. 

Suppose the dislocation $a[110] lies along the line [112], so that it is 
of pure edge type. The partial dislocations have equal edge components 
za[110], and opposite screw components ja[112] and ja[112]. A shear 
stress P acting across the glide plane in the glide direction [110] drives 
them both in this direction. A shear stress P acting in the direction 
[112] acts on the screw components, driving them apart with a force 
per unit length of Pa/2,/6. They will separate completely under a 
stress given by P=2,/6y/a. 

Dislocations in which the two partials are appreciably separated are 
referred to by Heidenreich and Shockley as extended dislocations. They 
do not have the same freedom of motion as complete dislocations, because 
the two partials have Burgers vectors which are not parallel, and there 
is no line along which they are simultaneously pure screws. Since only 
a pure screw can move out of the glide plane, there is an activation 
energy required in order to produce the reaction 


fa{121]+ ta[211]= $a[110]= af 121] + da[2T1], . yea 


which leads to a dislocation extended in the plane (111). 
Thompson (1952) has considered a dislocation 4a[110] which already 
lies partly in the plane (111) and partly in (111). The partial dislocation 
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ga[121] in (111) may continue as either of the partials ga[12T] or 4a[211] 
in (111). In the first case it igs linked to the join of the other pair of 
partials by a dislocation }a[001], and in the second case by a dislocation 
ga[110]. The pair of partial dislocations in (111) is separated by a 
stacking fault and the pair in (111) by a different stacking fault. The 
small dislocations 4a[001] and {a[110] lie along the junction of these two 
stacking faults. Thompson has also considered the behaviour of the 
stacking faults at a node corresponding to eqn. (55). 

If the partial dislocations travel independently through the lattice, 
leaving surfaces of misfit behind them, their motion is more constrained 
than that of a complete dislocation. A complete dislocation may be 
reflected at a free boundary with change of sign, or it may be followed 
across the glide plane by another similar dislocation. In either case the 
resulting displacement is 2b. If b is a partial dislocation, the displace- 
ment 2b does not lead to close packing, but to an unstable position of 
high energy. It is therefore necessary for the second partial dislocation 
to have its Burgers vector inclined to that of the first. The resultant of 
the two then lies along a <110> slip direction. Alternatively, the same 
displacement b may occur in successive glide planes of the cubic lattice, 
leading to the formation of a twinned lamella. In the hexagonal structure 
a displacement b which is produced between layers 1 and 2 by the passage 
of a partial dislocation cannot occur between layers 2 and 3. If the 
three possibie partial dislocations between the first two layers are 
b,, b, and bs, those possible between the next two layers are —b,, 
—b, and —b,;. The partial dislocation b, can travel only on alternate 
glide planes, and in doing so it creates a lamella of cubic crystal. 

Mathewson (1951) has considered the effect of partial dislocations on 
the law of resolved shear stress. He did not use the terminology of 
dislocation theory, but we present his-argument in this form. Figure 41 
represents a unit triangle in the stereographic projection. According to 
the analysis of Taylor and Elam (1923), a crystal deformed in tension 
will slip in a particular <110> direction on a particular {111} plane if 
the tensile axis lies within this triangle. If slip along [101] consists of 
successive movements along [211] and [112], the triangle ADE must be 
divided by a line BC into two regions, so that in ABC the resolved shear 
stress is greater in the [112] direction than in the [211] direction, while 
in BCED this relation is reversed. If a tension is applied with an axis 
in the area ABC, and gradually increased, the first effect is to move 
partial dislocations 3a[112] across their glide plane (111). Slip cannot 
continue unless the partial dislocations $a[211] then follow them across. 
the glide plane. If we make no allowance for the energy y of a stacking 
fault, this means that a tension applied in the triangle ABC will cause 
slip only when its component resolved into the (111) plane and the 
[271] direction reaches a critical value. It is usually assumed to be the 
component in the direction [101] which is constant ; if, in fact, the 
component in the [211] direction is constant, the [101] component varies 
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by about 4 of its value as the tensile axis moves through the triangle. 
The observed variation (Fenn, Hibbard and Lepper 1950) in 70-30« brass 
was only 1/27 of the mean value over a wide range of orientations. If we 
take into account the energy y of the stacking fault, the resolved shear 
stress required to heal a fault is less than that to create it. There must 
be some other locus B/C’ such that a tension along a direction lying 
in B'C’ which has a component resolved along [211] just large enough to 
create a stacking fault has also a (smaller) component resolved along [112] 
just large enough to heal this stacking fault. We may now extend 
Mathewson’s argument, and introduce the locus B”C” of directions of 
tension for which the making of a fault with a [112] partial dislocation 
and its healing with a [211] partial dislocation require the same tension. 
In AB’C’, the partial dislocations are dissociated by the stress, and it is 
the stress resolved in the ]211] direction which determines the strength. 


Fig. 41 
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A unit triangle of the stereographic projection. If the crystal is extended along 
an axis lying in the triangle ADE, slip occurs on one particular glide 
system; outside this triangle other systems are preferred. In ABC 
‘easy glide ’ along [112] occurs before ‘ easy glide’ along [211]. In AB"C” 
the applied stress required to create a stacking fault by this process is 
less than that required to heal an existing stacking fault by glide along [211]. 


Similarly, in B’C’ED, it is the [112] component which must reach a critical 
value. In B’B’C’C’, the partial which moves first is immediately 
followed by the other partial, healing the fault. The effective stress on 
the first partial is not only the applied stress, but also its interaction 
with the second partial, which pushes it further ahead. The partials are 
not dissociated, but travel together across the glide plane under the 
influence of the [101] component of the stress. Apparently this region 
B’B'C"C’ covers most or all of the unit triangle in the case of « brass. 
Since this is a material which readily twins and so presumably has a. 
rather low value of y, this will probably occur in most other face-centred | 
cubic metals. 
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5.22. Sessile Dislocations 

The idea of a sessile dislocation, which cannot glide in the lattice, was 
introduced by Frank (1949). Any partial dislocation with its Burgers 
vector inclined to the plane of the stacking fault is sessile; other 
dislocations, which can move by glide under moderate stresses, have been 
called glissile. We first discuss the type of sessile dislocation which 
Frank described. 

If Plate IX represents the (111) plane of a close packed cubic 
lattice, the dislocation shown may be regarded as having been produced 
by inserting two close packed half planes of atoms into the crystal, lying 
parallel to (111) and (111) respectively. The partial dislocation formed 
by inserting one of these half planes alone is perhaps the simplest 
imperfection of dislocation type that can exist in a close packed cubic 
lattice. If, for example, the crystal is considered to be made by piling 
(111) planes, a plane of type C may be inserted between any two 
consecutive planes A and B without spoiling the close packing. Similarly, 
a portion of a close packed plane may be removed. In the hexagonal 
packing, a portion of a faulted plane may be added, but the removal of 
a portion of a plane would-spoil the close packing. The Burgers vector 
of this partial dislocation is 4a[111], and it can only glide in a plane 
containing this direction. The dislocation only leaves a close packed 
lattice if it moves in a (111) plane, and this does not contain the Burgers 
vector. It is therefore locked in its position in the lattice. Partial 
dislocations of this kind can move by three processes: (a) the diffusion 
of individual atoms, (6) dissociation into mobile dislocations, (c) ineffective 
glides. Process (a) is slow. Process (b) can occur by a reaction such as 


4a{111}+4a[011]+4a[2T1], Re en 16351 (60) 


where all three dislocations lie along the line [110] which is the intersection 
of the plane (111) in which the first dislocation can lie and the third can 
glide with the plane (111) in which the second can glide. The energy 
increases from } to 2 in this process. On the other hand, if a wandering 
partial dislocation ja[211] meets the sessible dislocation, they can 
combine without change of energy to form a mobile complete dislocation, 
Ag an example of (c) we consider the result of removing part of a C sheet 
from a cubic lattice, giving ... ABCABCABA*BCABC... 
The plane marked with an asterisk can move between the B planes to a 
C position, its movement with respect to each being in a direction of easy 
glide. The resulting structure... ABCABCABCB*CABC... 
can move by a similar ineffective glide to... ABCABCABCAC* 
nub Cr... 2 2nd"s0 on. 

It seems certain that sessile dislocations will play an important part 
in the process of work hardening, because they form fixed sources of 
internal stress which interact with mobile dislocations. In addition to 
the dislocation 4a[111], the face-centred cubic lattice can contain the 
sessile dislocations 4a[001] and }a[110] which arise at the junctions of 
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stacking faults on different {111} planes. These last two dislocations 
do not exist in the hexagonal structure. The first can exist, but cannot 
be formed by the interaction of mobile dislocations, whereas in the cubic 
lattice the reverse reaction to (60) allows two mobile dislocations to 
combine, forming a sessile dislocation and releasing energy. If there is 
only one glide plane, any dislocations which come together by glide have 
their Burgers vectors in this plane. The results of their reaction also 
have Burgers vectors in the glide plane, and are not sessile. This explains 
the slow work hardening of hexagonal metals. 

All the partial dislocations of the close packed cubic lattice may be 
expressed as combinations of partial dislocations of the type ga[211]. 
We have already outlined the way in which Frank’s sessile dislocation 
can be formed from one of these and a complete dislocation. Cottrell 
(1952 c) has developed an idea of Lomer (1951) to show how a sessile 
partial dislocation 4a[110] may be formed from complete dislocations 
alone. Lomer considers a dislocation }$a[101] lying in (111) and a 
dislocation 4a[011] lying in (111). These meet along the line [110], 
combining to form a dislocation }a[110], with release of energy. The 
line and the Burgers vector define the slip plane as (001), which is not a 
usual slip plane for the face-centred lattice. The dislocation cannot 
extend in this plane, and therefore the resistance to its motion is likely 
to be larger than the resistance which a dislocation meets if it lies in {111 }. 
However, it is not sessile. Cottrell points out that it can dissociate 
according to the scheme 


4a[110}~+4a[112]+ ¢a[112]+ 4a[110], sf a: 


the energy falling from } to 7/18. The first two products move away in 
the planes (111) and (111) respectively. They leave stacking faults 
behind them, and are therefore not repelled to infinity, but only to 
distances of the order of 10a. The resulting group of three dislocations 
cannot glide. 

It is perhaps more accurate to imagine the initial dislocations as 
already extended in the planes (111) and (111). In this case the two 
partial dislocations ja{211] and ja[{121], one in each plane, combine to 
give ga[110], while the other two partial dislocations are unaltered by 
the reaction. 

The principal difficulty in relating this reaction to the experimental 
observations is that the sessile dislocations are formed on the line [110], 
which makes an angle of 60° with the slip directions [101] and [011], 
whereas the observed deformation bands are planes perpendicular to 
the active slip direction. 


5.3. THE CHOICE OF A GLIDE SYSTEM 


The passage of a dislocation of Burgers vector b translates part of the 
crystal a distance b with respect to the rest, and the resulting crystal 
is again perfect. It follows that b is a lattice vector joining equivalent 
atoms in the crystal. The energy of the dislocation is proportional to 62, 
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and we expect that the glide directions which will be active in practice 
will be those corresponding to the smallest permissible values of 2 
(Frank 1949 d). This is not because we expect a Maxwell—Boltzmann 
distribution of dislocations in thermal equilibrium, but because dislocations 
of large Burgers vector are mechanically unstable. Thus in a cubic 
lattice a dislocation 2a[100], with energy 4, dissociates into two 
dislocations a[100] widely separated, with total energy 2. The ease with 
which a dislocation can move through the lattice is probably also 
influenced by the value of 6. A given applied stress exerts on the 
dislocation a force proportional to 6. Some of the forces resisting its 
motion are also proportional to b, but others depend on higher powers of b, 
and these make the mobility a decreasing function of b. 

These considerations determine the glide direction. The choice of glide 
plane is not so clear, which agrees with the observation (Andrade and 
Chow 1940) that the glide direction is usually a fixed characteristic of 
the crystal class, while the glide plane may be a function of temperature. 
Chalmers and Martius (1951) have shown that the glide plane is always 
the most closely packed plane which contains the Burgers vector; if 
several planes are observed, these are the most closely packed, the next 
most closely, and so on. They explain this by assuming that if the 
spacing of the glide planes is a, while the Burgers vector is b, the energy 
is an increasing function of b/a. It follows from (6) and (7) that on 
Peierls’s approximation this argument is not correct. The energy is 
proportional to b?, but depends on a only through the logarithmic term 
inf. The change of energy for a moderate change in a is likely to be small 
compared with the change produced by the elastic anisotropy when the 
glide plane is rotated round the Burgers vector, and this latter effect 
has been neglected. Probably the most important effect is the deviation 
of the law of force between planes of atoms from the sinusoidal form. 
This allows the dislocation to lengthen and so reduce its energy and the 
resistance to its motion if it lies in certain planes. This is particularly 
effective in close packed planes, where the dislocation can extend into two 
partial dislocations. 

5.31. Special Lattices 

5.311. The Simple Cubic Lattice—In the simple cubic lattice of 
parameter a, all atoms are equivalent. The shortest possible lattice 
vectors are of the forms a[100] and a[110]. If the dislocation a[110] 
dissociates according to the scheme : 


a{110}+a[100]-+-a[010], Mae er eM ( 62) 


the strain energy is unaltered, and more delicate considerations would be 
required to decide whether the dislocations present after annealing are 
<110> or <100>. However, since the dissociation (62) requires no 
energy, a stress acting on a dislocation a[110] which produces a larger 
force on the component a[100] than on the component a[010] is able to 
move the component a[100] while the other is still fixed in the lattice. 
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The glide direction is thus necessarily <100>. There are three basic 
slip vectors a[100], a[010] and a[001] which are mutually perpendicular ; 
nodes and reactions such as (55) are impossible. Salts having the 
caesium chloride structure are crystallographically simple cubic. 
Thallium halides, which have this structure, do in fact show <100> as 
the glide direction, with {110} as the glide plane (Smakula and Klein 
1949). These salts cannot contain partial dislocations of low energy 
because of the electrostatic repulsion of like ions. No dislocation 
reactions of any kind are possible, and a dislocation moving in a particular 
glide plane can only continue to move in this plane. This is probably 
the explanation (Nabarro 1950) of the ‘ prismatic punching * observed 
in these salts: an indentation on one face of a crystal travels through the 
crystal and produces mounds on the other faces at points lying in 
<100> directions from the indentation. (Seitz (1950c) has given a 
different explanation. It could also be an example of kinking (Jillson 
1950).) 


_ 5.312. The Face-centred Cubic Lattice.—The shortest lattice vector is 
Za[110], and the glide direction is <110>. In the approximation of an 
isotropic elastic continuum there is no change of energy if the dislocation 
a{100] dissociates into 4af[110]+4a[110] or into 4a[101]+4a[101]. 
However, the last four dislocations lie in {111} planes in which they can 
extend and reduce their energy, whereas a[100] does not lie in this plane. 
Thus the glide plane {111} and glide direction <110> are determined. 


5.313. The Body-centred Cubic Lattice—The shortest lattice vectors 
are a[100] and }a[111], with energies 1 and #? respectively. The glide 
direction is thus <111>, and the glide plane is variable. A dislocation 
a[100] cannot dissociate spontaneously into $a[111]+4a[111], because 
the energy rises from 1 to 3. However, two such dislocations may be 
removed, with release of energy, by the reaction 


a{100]+-a[010}>4a[111]+4a[1lT]. . . . . (63) 


5.314. The Hexagonal Close Packed Structure-—In the close packed 
hexagonal structure (and the structure derived from it by changing the 
ratio c/a) all atoms in the same basal plane are equivalent, but atoms in 
the planes above and below are not equivalent to these. The normal 
glide directions are therefore the close packed directions in the basal 
plane, with hexagonal indices <1120>. Suppose the stress is so 
oriented that the shear stress acting across the glide plane is small, which 
means that the hexagonal axis lies close to a principal axis of stress. 
Glide cannot occur on the base plane. The shortest remaining lattice 
displacement between equivalent atoms is in the direction of the hexagonal 
axis, from one atom to an atom in the next close packed plane but one, 
with Burgers vector [000c]._ This is a principal axis of stress, and there 
is no component of shear stress tending to move this dislocation. Any 
other complete dislocation dissociates into dislocations of the types already 
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considered. A hexagonal crystal pulled in, or perpendicular to, the 
basal plane, cannot deform by slip. It either fractures or twins by the 
motion of partial dislocations. 


5.315. The Sodium Chloride Structure—This structure has been 
discussed by Seitz (1950). The Burgers vector for slip is the line joining 
two neighbouring ions of like sign, 4a[110], and this defines the slip 
directions <110>. 


5.4. TWINNING AND PHASE TRANSFORMATIONS 


It is obvious that suitable movements of partial dislocations can 
produce the shear displacements which occur in mechanical twinning 
and in some phase transformations in the solid state. Cottrell and 
Bilby (1951) have studied mechanical twinning in the body-centred 
and face-centred lattices from this standpoint, while Millard and Thompson 
(1952) and King (1952) have considered twinning in the hexagonal close 
packed structure. Frank (1952b) has discussed the martensitic 
transformation from face centred to body centred iron; in this case the 
mechanism of the homogeneous shear is not discussed, but the 
inhomogeneous deformations which also occur are most naturally 
expressed in terms of screw dislocations. Stacking faults in cobalt and 
the relation between the cubic and the hexagonal close packed structures 
have also been studied (Christian 1951). A treatment of twinning from 
a rather different standpoint is due to Lifshits and Obreimoy (1948). 
Before discussing the theories of Cottrell and Bilby, and of Lifshits and 
Obreimovy, we make some general observations on the dislocations 
associated with twinning. We consider only twins in crystal lattices 
in which all atoms are equivalent, and neglect twins of the second kind 
in which the twinning plane is irrational. 


5.41. Twin Lattices, Twinning and Pole Dislocations 

Consider any lattice of atoms generated by three vectors which are 
not coplanar. Suppose (fig. 42) that a and b are the vectors joining 
any atom O to any other two atoms A and B. If these two vectors are 
not parallel, and the lines OA and OB do not pass through any other 
atoms, the vectors a and b are basic lattice vectors. They generate a 
two dimensional lattice, and the three dimensional lattice can be resolved 
into a stack of similar and similarly situated two dimensional lattices. 
The third basic vector c is the line joining O to any atom in the next 
sheet ; we choose in particular that atom C whose projection normal 
to the ab plane is a point P lying within the two dimensional unit cell 
defined by OA and OB. If h is the normal distance between sheets, 
and p is the vector OP, the basic vectors are a, b and h-+-p. There is 
one other point P’ in the two dimensional unit cell which is symmetrically 
related to P. By regarding D as the origin it is easily seen that the 
vectors —a, —b and h+p’—a—b define a lattice which is the mirror 
image ot the first, the twin lattice. We may write this set of basic vectors 
more simply as a, band h—p. (The mirror image may be regarded either 
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as a mirror image in the plane OAB or as an image in the line perpendicular 
to this plane and passing through the mid-point of AB. For a lattice 
these are identical; for a structure with basis and no centre of symmetry 
these processes yield different twins.) The coordinates of any atom in the 
twin lattice differ from those of a point in the original lattice by a multiple 
of the twin vector t=p’—p=—2p. The partial dislocation with Burgers 
vector t is a twin dislocation. If it lies in an ab plane of the original 
lattice, it is mobile, and may be called a twinning dislocation. As it 
moves across the plane, it leaves behind it a sheet of stacking fault, which 
may be regarded as a layer of twin one atom thick. For certain special 
choices of the vectors a and b, this stacking fault has a low energy, so 
that the corresponding imperfections of the crystal are actually observed. 
A twinning vector t may, in a lattice of high symmetry, be common to 
two twinning planes, as occurs in the body-centred cubic structure. 


Fig. 42 


0 b B 


Twinning. O, A, B and D are lattice points, a and b lattice vectors in a plane. 
P is the normal projection of a point in a neighbouring plane. To the 
point P corresponds the symmetrical point P’, which is the normal 
projection of a point in the twin lattice. 


To produce a twinned layer m sheets thick, twinning dislocations 
must pass across each of the m planes separating one sheet of atoms from 
the sheet below. Cottrell and Bilby pointed out that this could occur by 
the action of a modified form of the Frank—Read source. In the ordinary 
Frank—Read source, the dislocation in the glide plane rotates about a 
fixed pole dislocation. The pole dislocation may be the moving 
dislocation itself which has turned into another line and is there anchored, 
or there may be a dislocation node. In either case, all the Burgers vectors 
lie in the glide plane of the moving dislocation. If the Burgers vector of 
the pole dislocation has a component n normal to this glide plane, the 
stack of glide planes is distorted into a helical surface of pitch n. Each 
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time the moving dislocation rotates, it finds itself displaced a distance n 
normal to the glide plane. If n=h, the normal distance between planes 
of atoms, and the moving dislocation is not complete, but a twinning 
dislocation, a thick layer of twin can be formed. 

We now consider the possibility of finding a suitable pole dislocation. 
In general, a dislocation of strength h cannot exist, for this would imply 
that a sheet of atoms originally lying distances h+-p from the atoms in 
the reference plane had been shifted to 2h+-p, while another sheet was 
interposed. The interposed sheet must lie either h+p or h—p from the 
atoms in the reference plane; the vector from an atom in the interposed 
sheet to an atom in the third sheet is then either h or h+2p. The vector 
h is not allowed, while h+2p is allowed only if 3p is a lattice vector, 
so that h+2p—h—p is a vector of the twin lattice. This occurs for 
the (111) plane of the cubic close packed lattice. If the dislocation h 
does exist, it is sessile, but cannot function as a pole for twinning on the 
plane to which its Burgers vector is normal, or on other planes of the 
same type, for it lies in its own plane and does not thread it, while the 
component of its Burgers vector perpendicular to the planes of its own 
type which it threads is necessarily less than h, their spacing. 
Dislocations of strength h+p or h—p are possible, since the atoms then 
lie 2h+ 2p or 2h from atoms in the reference plane, allowing the interposed 
atoms to lie at h+p or h-+p from the atoms in the reference plane. 
However, h-+p and h—p are lattice vectors of the original and the twin 
lattices respectively. The corresponding dislocations are therefore 
complete, and cannot be sessile. Their mobility may well interfere with 
the twinning process. Sessile dislocations of strength 2h can always 
exist, because their atoms lie 3h+p from the atoms in the reference 
plane, which allows two interposed sheets to lie at h+p and 2h, h—p 
and 2h, or h+p and 2h+2p. The dislocation 2h can be resolved into 
h-+-p and h—p. The former is mobile in the original lattice, but the 
dissociation will not occur if the angle between h+p and h—p is obtuse, 
that is, if h<p. The line of the dislocation lies in the glide plane of the 
twinning dislocation —2p, but if this plane is inclined at 60° to another 
plane of the same type, as happens in the body-centred cubic lattice, 
the component of its Burgers vector normal to this other plane is h, 
the interplanar spacing. 


5.42. Special Lattices 

5.421. The Body-centred Cubic Lattice.—Cottrell and Bilby consider a 
complete dislocation }a[111] which lies in (112) but not along [111]. 
This dissociates without change of energy according to the scheme 


fa[111l}-$a[112]+gaflll]... . . . . - (64) 


The first component is pure edge and sessile. It is of the type 2h, and its 
plane (112) is inclined at 60° to each of the equivalent planes (121) and 
(211). It threads each of these, and its Burgers vector has a component 
normal to each of h=a/,/6. The second component is a twinning 
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dislocation in (112), and can rotate into the direction [111], leaving a 
stacking fault in (112) behind it. In this orientation it is a pure screw. 
It is also a twinning dislocation in (121) and (211), and under a suitable 
applied stress it starts to move in one of these. The requirements for 
progressive twinning are now satisfied. 


5.422. The Face-centred Cubic Lattice —Here Cottrell and Bilby state 
that the dissociation analogous to (64) is 


4a[110}>-4af111]--hel 12] 


They show, by a detailed analysis, that progressive twinning cannot 
occur because the sessile dislocation 4a[111] lies in the glide plane of the 
twinning dislocation ja[112]. However, the analogy between (64) and 
(65) is incomplete, for 4a[111] is a sessile dislocation of type h, and the 
true analogy is with the sessile dislocation 3a[111], which is of type 2h. 
Since the cosine of the angle between {111} planes is 4, this dislocation 
gives a lift of only 2h to twinning dislocations in another {111} plane. 
A dislocation a[111], which would give the necessary lift, would be both 
mobile and unstable. It seems that the absence of mechanical twinning 
in the face centred cubic lattice corresponding to the (111) twins produced 
during recrystallization is due to the lack of any suitable dislocation 
mechanism. 


5.43. The Analysis of Lifshits and Obreimov 

In experiments on calcite, Garber (1947) found that indentation with 
a knife edge produced a narrow twin which progressed gradually into the 
crystal under an increasing force, and reverted to the original lattice 
when the force was removed. When greater loads were applied, the twin 
progressed further, and did not disappear on reducing the load. Finally, 
the twin ran right across the crystal when a critical load was applied. 
Lifshits and Obreimov have analysed the process using a simple atomic 
model in which the restoring forces on an atom are linear in its displace- 
ment until a critical displacement is reached, after which the restoring 
force is proportional to the displacement from a second, twinned, position 
of equilibrium. They show that the Green’s function for displacements 
may be expressed in terms of the lattice vibrations. The method of 
solution is to postulate a certain shape for the boundary of the twin, and 
verify that the relative displacements of atoms near the boundary still 
lie within the elastic range assumed. (Most of the calculations are, in fact, 
carried out for the case of a force applied to a line of atoms within the 
crystal.) ‘They find that three cases are possible, according to the value 
assumed for certain parameters of the material: (1) All twins are elastic, 
and reach a finite depth under a finite load. On removing the load, they 
revert to the original lattice. (2) Twins penetrating to less than some 
critical depth are elastic, but those which penetrate further remain as 
twins when the load is removed. (3) At some value of the applied load 
the twin propagates into the crystal without further increase in the load. 
There may or may not be an elastic region before this. 
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ADDENDUM 


“Properties of Helium Three at Low Temperatures’, J. G. Daunt, 
Advances in Physics, 1, 209, 1952. 


In §6(c) of the above article reference was unfortunately omitted to a paper 
by L. Goldstein and M. Goldstein (J. Chem. Phys., 18, 538, 1950) which 
considers theoretically the magnetic properties of liquid *He. In this paper 
the exchange energy originating in the interatomic mutual potential energy of 
the *He atoms (both Slater-Kirkwood and Margenau potentials are considered) 
is computed and is compared with the kinetic energy of the system. It is 
concluded, as reported in §6(c), that nuclear ferromagnetism in liquid *He is 
not a stable configuration. 
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